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1 Pin-jointed frameworks 


1.1 Structures of least volume of material 


The data for the present problem consist of (i) given forces applied at given 
points in space and specified by Cartesian components, (ii) points of support 
whose forces of constraint are together capable of reacting any given force and 
couple, and (iii) a pin-jointed framework of any given layout, whose nodes 
include the points of application of the given forces and the points of support, 
and whose members form a structure as opposed to a mechanism. The given 
forces may be thought of as embracing all components of force (including zero 
components) acting at all the nodes of the framework, excluding only component 
forces of constraint. This complete set of given force components is denoted by 
F,(=1, 2 ..n). The points of support must be sufficient to prevent overall 
motion of the structure, but may be redundant. End loads 7;( = 1, 2 . . m) in 
the m members of the framework and forces of constraint at the supports can 
always be found to balance any set of given forces [Fj] .¢ An independent set of 
nodal equilibrium equations for [7;] can be written 


. KjT;=F = G=1,2--n), (1.1) 
i=1 


where [Kj] is determined by the direction cosines of the members of the frame- 
work and is therefore known. The order of redundancy of the framework is 

m~n and should be greater than zero, so as to provide choice for the layout of 
an optimum structure. The matrix [K;;] will have rank 7, since (1.1) is always 
soluble for 7 of the 7;. 

The framework of (iii) can be given any practical degree of generality by an 
appropriate choice of its nodes and members. It is completely determined when 
the areas of cross-section of its members A ;(i = 1, 2 . . m) are given, as well as 
its material of construction. The magnitudes A; must satisfy 


A; 20 @=1,2..m), (1.2) 
and the ith member is present if A; > 0 and absent if A; = 0. 
+ This can always be achieved, even if it is not specified, by adding new supports and 
perhaps converting some given forces into forces of constraint. 


} The meaning of the symbols 7; and F; is illustrated for a particular structure in Fig. 1.1 
on p. 7. 
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1 Pin-jointed frameworks 


1.1 Structures of least volume of material 


The data for the present problem consist of (i) given forces applied at given 
points in space and specified by Cartesian components, (ii) points of support 
whose forces of constraint are together capable of reacting any given force and 
couple, and (iii) a pin-jointed framework of any given layout, whose nodes 
include the points of application of the given forces and the points of support, 
and whose members form a structure as opposed to a mechanism. The given 
forces may be thought of as embracing all components of force (including zero 
components) acting at all the nodes of the framework, excluding only component 
forces of constraint. This complete set of given force components is denoted by 
F,(=1, 2 ..n). The points of support must be sufficient to prevent overall 
motion of the structure, but may be redundant. End loads 7;( = 1, 2 . . m) in 
the m members of the framework and forces of constraint at the supports can 
always be found to balance any set of given forces [Fj] .¢ An independent set of 
nodal equilibrium equations for [7;] can be written 


. KjT;=F = G=1,2--n), (1.1) 
i=1 


where [Kj] is determined by the direction cosines of the members of the frame- 
work and is therefore known. The order of redundancy of the framework is 

m~n and should be greater than zero, so as to provide choice for the layout of 
an optimum structure. The matrix [K;;] will have rank 7, since (1.1) is always 
soluble for 7 of the 7;. 

The framework of (iii) can be given any practical degree of generality by an 
appropriate choice of its nodes and members. It is completely determined when 
the areas of cross-section of its members A ;(i = 1, 2 . . m) are given, as well as 
its material of construction. The magnitudes A; must satisfy 


A; 20 @=1,2..m), (1.2) 
and the ith member is present if A; > 0 and absent if A; = 0. 
+ This can always be achieved, even if it is not specified, by adding new supports and 
perhaps converting some given forces into forces of constraint. 


} The meaning of the symbols 7; and F; is illustrated for a particular structure in Fig. 1.1 
on p. 7. 


2 Pin-jointed frameworks 


If the structure defined by [Aj] is to transmit the forces [Fj] to the supports 
with safety, the stresses corresponding to [7;] must be restricted in magnitude. 
Let the ‘allowable stresses’ for a chosen material be oy in tension and gg in 
compression. It then follows that 


T; <o7Aj, —-T; <0¢A; (@@=1,2..m). (1.3) 


The relations (1.3) must be satisfied even if [7;] is found by an ‘elastic analysis’ 
of the framework. However, for the moment, the simpler procedures of ‘plastic 
analysis’ will be followed. They require only that a solution of (1.1) can be 
found which satisfies (1.3). It will be shown later that this limited approach does 
in fact lead to an optimum elastic design, which in the case of a single system of 
given forces coincides with that obtained by plastic or limit analysis. 

The optimum problem of the present section can now be formulated. It is to 
find that framework defined by [Aj] , which will safely transmit the given forces 
[4] to the given supports and make the volume V of material in the structure 
given by 


y= 3 Aili, (14) 
i=1 


where /; is the length of the ith member, a minimum. The conditions of safety 
require that non-negative [A;] be chosen so that (1.1, 3) have a solution [7;] . 
Mathematically the problem is that of minimizing V of (1.4), subject to the 
restrictions (1.1, 2, 3). 

The relations (1.1—4) can be made non-dimensional by taking [F;] propor- 
tional to F (a typical force), [7;] to F, a7 and og to o (mean of oy, dc), [Aj] 
to F/o, {u} to / (a typical length) and V to Fi/a. This means that the best 
material is that with the largest o. If a structure of least weight is required, then 
pV must be a minimum, where p is the specific weight, and so the best material 
is that with the greatest ‘specific stress’ o/p. 

The relations (1.3) may be replaced by 


T; +27; =0yA;, -T;+2T/=0¢A;, 1/20, T,>0 (=1,2..m); 
(1.5) 


where [7;] , [7;"] are ‘slack variables’, and the first two of (1 -5) solved to give 


A,=(Tj+T Jo, 7, = (077; - ocT;"io (i= 1,2..m), (1.6) 
where, 


0 = (0% + 0¢)/2. (1.7) 
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Substituting from (1.6) into (1.1) and (1.4) then gives a new formulation of the 
present problem in the form: 


min V= ¥ (7; +T;")i/o, (1.8) 
i=1 
subject to 
TOTO oh=42 or me (1.9) 
and 
> Kiy(orT; - o¢T,;")/o=F (/=1,2..n). (1.10) 
i=1 


Equations (1.8, 9, 10) define a problem of linear programming. 


1.2 The simplex method 
The problem of (1.8—10) is a special case of the standard problem of linear 
programming: 


P 
minz= > ¢;xi, (1.11) 
i=1 
subject to 
x20. G=1,2.-p); (1.12) 
and 
P 
> ayx:=b;  (G=1,2.-n), (1.13) 


where {c;}, [a3] and [b;] are given matrices. The linear function z of (1.11) is 
called the ‘objective function’. The relations (1.12) are termed the ‘non- 
negativity restrictions’ and the equations of (1.13) the ‘constraints’. Non-negative 
values of the variables [x;] are to be found which satisfy the constraints and make 
the objective function a minimum. 

It is assumed that [a;;] has rank n and so (1.13) can be solved, after a renum- 
bering of the variables perhaps, in the form : 


P 
Xj SNe ey xrees 0H 1; 2s=n)- (1.14) 
i=n+1 
This implies that 
xy=Xpy (G=1,2..n), x,=0 (i=nt1,.-.p) (1.15) 


is a solution of (1.13). It is termed a ‘basic solution’, meaning that it contains at 
most 1 non-zero components, The variables x;(j = 1, 2 . . 2) are termed “basic 
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variables’ and x,;(i= + 1, . . p) ‘non-basic variables’. Any solution like (1.14) 
divides [x;] into variables of these two kinds. 
Assume now that the constants xp;{j = 1,2 . . m) satisfy 


Xp PO G=12 ssn) (1.16) 


The solution (1.15) then satisfies (1.12) and is for this reason termed ‘feasible’. 
It is thus a ‘basic feasible solution’. It is always possible to find such solutions of 
(1.10), since any solution of (1.10) can be taken as giving for each of 7 values 
of i either 7; or —7;,". 

Equation (1.14) is equivalent to (1.13) and gives feasible solutions, with 
x;20(i=n+1,..p),solong asx; 20 (j= 1, 2 . . ) is satisfied. Substitution 
from (1.14) into (1.11) gives 


P 
z=zpt ¥ (-2)%, (1.17) 
i=n+1 
with 
n 
Zp= > 6 XB; (1.18) 
j=1 
and 
n 
Z=> Gy (i=n+1,..p). (1.19) 
j=1 


The quantity zg is the value of the objective function corresponding to the 
basic solution (1.15). The quantities z; (i= + 1, . . p) can be obtained once a 
solution (1.14) has been found and are thus defined for any basic solution. 

Sufficient conditions for (1.15) to give a minimum of z equal to zg are now 
seen by (1.17) to be 


Cj 22; (@i=n+1,..p), (1.20) 


since, if (1.20) is satisfied, no non-negative x; (i =n + 1, . . p) can decrease z in 
(1.17) below zg. a 
If (1.20) is false for at least one of i=n + 1, . .p, say i=K, that is if 


Ck <Zr, (1.21) 


then it is possible by giving x; positive values and setting x; = 0 (i # k) to 
decrease z below zg. The variables x; of (1.14) are then given by 


Xj =Xpj —VejXn 20 G=1,2..n), (1.22) 


and must satisfy the inequality imposed. This means that 
Xk SX pjlVeej, for all / with yj > 0. (1.23) 
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There is no restriction on x, if yp; <O (j= 1,2. . n) and in this case, by (1.17), 
Zz —© as X, > +°°, Such ‘unbounded solutions’ do not arise in structural 
applications; V in (1.8) cannot be negative. Assume therefore that there is a 
value of r taken from j = 1, 2 .. 2 such that 
Xr/Ver = min XBj/Yrj, for y,j > 0. (1.24) 


The assumption that x, = xg,/y;,, then satisfies all of (1.22) and the solution 


Xj =Xpj—VejXpr/Yer 9 (J =1,2..n; JF), 


Xp=Xpr/Vers 

25 
x=0 8 (fj=nt1,..p; 7#K, ae 
x, =0, 


in which the previously basic variable x, has been interchanged with the 
previously non-basic variable x,, is a new basic feasible solution. The corres- 
ponding value of z is by (1.17) 


Z=Zp t (Ce — Zp) XBr/Ver, (1.26) 


which thanks to (1.21) gives a smaller value of z than zg, unless xg, = 0. 

A basic solution with one or more of its basic variables equal to zero is said 
to be “degenerate”. If (1.15) is degenerate and Xp, = 0 for a value of j for which 
Ynj > 0, then by (1.24) xg, = 0 and no decrease of z is possible by inserting x, 
among the basic variables. However z is not increased by this insertion and a 
new set of basic variables x; (j= 1,2 ..; 7 #r), x, has been introduced. 

It is to be noticed that ambiguity in the choice of r in (1 -24) leads to 
degeneracy in the new basic solution of (1.25), since at least one of the 
equations in the first line of (1.25) will give a zero value. 

The process of interchanging basic and non-basic variables described above 
is called a ‘simplex transformation’ or ‘simplex iteration’. It can be continued 
as long as a k can be found, which is the index of a non-basic variable for which 
(1.21) is valid. It will eventually determine an ‘optimum solution’, for which z 
is a minimum, or an unbounded solution, because these are the only ways it can 
stop and the number of basic solutions like (1.15) are finite in number. A 
possible difficulty is ‘cycling’, which can occur if a basic solution repeats itself 
during the iterative process. This can only happen if there is degeneracy, because 
otherwise z is decreased at each iteration. It can be overcome by an appropriate 
selection of basic and non-basic variables for an interchange, in cases where there 
is more than one possibility. Fortunately cycling arises but rarely in practical 
calculations. Means for avoiding it are given in the standard works on linear 
programming (Dantzig 1963; Hadley 1965). 

Equations corresponding to (1.14), for the new basic variables of the 
solution of (1.25), can be obtained by solving the rth equation of (1.14) for Xt 
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and substituting in the remaining equations. This gives 


P 
Xj =(XBj Yai XBr/Ye— | Ow — Yj Virl Var) Xi + (Dri /Ver) *r 
gr 

e (j=1,2..0; 7 #7), (1.27) 

Xe =Xpr/Yur— { >  Oir/Yer) xi + *lnr| > 
i=n+1 
i#k 
which can be used like (1.14) to carry out a new iteration if required. 
The choice of k is governed by (1.21). If there is ambiguity, then (1.26) 

suggests that the most rapid approach to the minimum of z is obtained by taking 
that k, which makes 


(Ze —Ck)Xpr/Ver a Maximum. (1.28) 
In practical calculations the simpler criterion 
Ze —Ce iSamaximum (1.29) 


is usually used. 

The process of iteration described above is called the ‘simplex method’ of 
resolving problems of linear programming. It is due to G. B. Dantzig and is 
readily programmed for the digital computer. Many standard programmes are 
available, see for example Kienzi, Tzschach, and Zehnder (1968), and so there 
is no difficulty in solving quite complicated problems of the type defined by 
(1.8—10). The only limit is available storage. 

It has already been remarked that (1.10) may be solved in the form (1.14) 
with (1.16) satisfied. Any statically determinate sub-framework, preferably one 
which experience suggests is a “good design’, may be used to isolate m quantities 
T,, including the loads in this sub-framework, and then (1.1) may be solved for 
these variables in terms of the given forces [F;] and the remaining m —n 
quantities 7;. This defines a solution to (1.10) as well, since this is the same set 
of equations, and leads to a division of T;’, 7;" (f= 1, 2 . . m) into basic and 
non-basic variables, with an associated basic feasible solution, because a non- 
negative term in this solution can always be obtained by choosing one of T; or 
T," as the basic variable. It is to be remarked that if one of 7;', T;” is a basic 
variable, then the other must be non-basic and so, for all i, one of 7;,, 7,” is 
always zero in a basic solution. Both can be zero of course, if 7;' and T;' are 
non-basic variables. 

The simplex method leads to a definite minimum value for V, since, as has 
already been remarked, unbounded solutions are not possible because V > 0. 
The basic solution corresponding to Vinin defines an optimum framework and 
so demonstrates the existence of such a framework of least volume of material 
designed according to the methods of plastic design. This optimum framework 
is statically determinate and so is an elastic design as well as a plastic one. Now 
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since elastic designs have to satisfy the additional conditions of compatibility 
of strains, the minimum value of V for elastic designs cannot be less than that 
| for plastic designs determined here. However the optimum plastic design is an 
elastic design too and so its value of V is both < and > the minimum V for 
elastic designs. The two minima thus coincide and the least volume for elastic 
designs is equal to that for plastic designs. 


1.3 An example 


Consider the symmetrically loaded symmetrical framework of Fig. 1.1. The 
figure shows the most general possible system of given forces and a complete 


Fig. 1.1 


system of end loads in the members. An optimum design is sought for the 
special loading system 


F, =F, =F; =F, =0, Fs =F. (1.30) 
The equations of nodal equilibrium corresponding to (1.1) are 
T, + Tz//2 + Ta//5 = F, = 0, 
T3//S — Ta/\/5 — Ts = Fz = 0, 
-2T3//S — 2Ts/\/5 = F3 = 0, (1.31) 
-T,//2 — isla —0, 
2T3//5 + Ts = Fs =F, 


and have a solution 


T, =F-Ts/V5, Tz =—J/2F + 2V(2/5)Ts, Ts=-Ts, Ts =-2Ts/V/5, 
Ts=F-2Ts//5. (1.32) 
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Assume for simplicity that o¢ = 0¢ = 0. Equation (1.6) then gives 
Ais (Tet Te ont 2s ay CE15 208): (1.33) 
Equations (1.32) then yield, using (1.33), 
Ty =F+Ty' —(T4 —Ta')/5 = 0A, -Ty' 
Ty’ =/2F + Tz — 2V(2/5)(T4 - Ta’) = 042 - Th , 
T;' =T3 + Ts —Ts =0A3—T3, (1.34) 
Ts! = 2(Tq —Ta')/5 + Ts = 0A5 —Ts, 
Te =F -2(T4 -Ta V5 + Ts = 0A6—Te . 
Also, 
oAg=Ts + Ty’. (1.35) 
The choice of basic variables is clear in the case of T;, Tz’ and T¢, since these 
take positive values in the basic solution corresponding to (1.34). T3' and Ts! 
are chosen, since if members 3 and 5 do come into operation they will be in 
compression. 
The values of /; are given by 
4,= 11, V2, V5/2, V5/2, 1/2, 1). (1.36) 
Substituting from (1.34), (1.35), (1.36) into (1.4) or (1.8) then determines V 
for the complete structure as 
V=(2l/o)(4F + 2Ty' + 2/27; + V5T3 — Ta [s/5 + 6Tq'/\/5 + Ts + 2Te'). (1.37) 


The coefficient of Tj is negative and so the basic solution corresponding to 
(1.34) with V = 8/F/o does not give a minimum. The non-basic variable Ty must 
therefore be exchanged for a basic variable. The rule of (1.24) gives Ty ore 
for the exchange. Taking 7’ and solving the second of (1.34) for Ty gives 

Ta = V5F/2 + V(5/2)(T2 — Tr’)2 + Ta, (1.38) 


which when substituted in the remaining equations of (1.34) and in (1.37) gives 


Ty = F/2+T' — (Tz —T 22, 

T3' = V5F/2 + V(5/2)(T2 — Ta')/2 + Ts, 
Ts = F + (Tq —Ts' 2+ Ts, 

Te =—(Tq - T,' 2+ Te, 


(1.39) 


and 


V = (2I/0)(7F/2 + 2Ty) + 7Tq [2/2 + Ta! |2\/2 + VST3 +574’ + Ts + 2T¢'). 
(1.40) 
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The new basic solution gives an optimum, since V in (1.40) cannot be 
decreased below 7IF/o. The minimum V is thus 


Venitn = OO (1.41)+ 
with 
Ty’ = Ty = Ty' = T3 = Ty! = Ts =Te' = 0. (1.42) 
Substituting from (1.42) into (1.38, 39) and using (1.33) then gives 
T; = F(1/2, 0,-/5/2, V'5/2, -1, 0), (1.43) 
and 
A; = (F/a)(1/2, 0, V'5/2, V'5/2, 1, 0). (1.44) 


The optimum design is thus determined. The layout follows from Fig. 1.1, 
where now members 2 and 6 should be omitted. The node at the force 2F4 
should also be removed to avoid instability. The present theory does not require 
this, since it takes no account of such matters. 


1.4 Duality in linear programming 


The problem of (1.11—13) is now called the ‘primal problem’ and is assumed to 
have an optimum solution given by (1.15). Equations (1.13) are soluble for 

x, (= 1,2...) and this process may be made explicit by introducing co-factors 
Ay Ui, K=1,2..7), which satisfy 


n 1@=k) 
> ayAjn = i, k=1,2..n). 1.45 
24 Voge GETHIN (1.45) 
The constants of (1.14) can then be written 
Xpr = > dj Ajp (K=1,2..27), (1.46) 
j=1 
and 
n 
Vir = ¥ aijAjn @=nt1,..p; k=1,2..7). (1.47) 


j=l 


The minimum value zg of (1.18) can be expressed, using (1.46), as 


n n 
Zmin=Zp= > > dyAince (1.48) 
k=1j=1 


| This may be compared with the absolute optimum for the half-plane, which by (4.56) has 
Vinin ® 2nlF/o, 
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and the condition (1.20), using (1.19) and (1.47), can be written 


n n 
> > wang (Gi=n+1,..p). (1.49) 
R=t=2 


The ‘dual problem’ or the ‘dual’ of (1.11—13) can be written as: 
n 


maxw= > bu, (1.50) 
j= 


subject to, 
n 
> agupse; @=1,2..p). (1.51) 


The relation (1.51) can also be written as 
n 
DY aut y=c;, 1,70 @=1,2. .p) (1.52) 
j=1 
The dual problem depends on the same constants as the primal, but in a different 
way. The dual variables {uj} , which can have any sign, are m in number, and 
this by (1.13) is the number of constraints of the primal. The number of con- 
straints in the dual, which by (1.51) is p, is equal to the number of variables 
[x;] in the primal. 
Substitution from (1.13) into (1.50) and use of (1.52) and (1.11) gives 


n Pp n Pp P ‘ 
w= > dyujy= > Y xayuj = > (c;X; — x j0;)) =z —- Dd x, (1.53) 
j=1 i=l i=1 


i=1 j=1 
and so, since by (1.12), (1.52) x;, v; @= 1, 2 . . p) are all non-negative, 


n P 
w= ¥ bu; <Zmin S > 0X, = 2, (1.54) 
i=1 i=1 


where {uj} is a solution of the dual constraint (1.51) and [x;] is a feasible 
solution of the primal constraint (1.13). The relations (1.54) give a means of 
finding upper and lower bounds to Zin. They also show that w is bounded 
above. 

It can now be shown that the dual has an optimum solution given by 


n 
u= > Ajece (f= 1,2.-n). (1.55) 
k=1 


Substitution in (1.52) gives 


n n 
YX a*yAjnce + V4; =; (i=1,2..p), (1.56) 
j=1 k=1 
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and so by (1.45) and (1.49) 
y=0 (@=1,2..n); u,20 (=n+1,..p), (1.57) 


which shows that (1.55) satisfies the dual constraints (1.52) or (1.51). The value 
of w corresponding to (1.55) is, using (1.48), given by 
n n 
Wi SSO Ape 2 min» (1.58) 
k=1 


fad 
and so by (1.54) it follows that the dual has a solution and that 
Wmax = Zmin- a 59) 


The value of Z,,in Can thus be obtained by solving the dual problem. 
Equation (1.53) gives for the optimum solutions to the primal and the dual, 


Pp 
since (1.59) is valid for them,  x,v; = 0 or since each term of this sum is 


sat i=1 
positive, 


x20 (=1,2..p). (1.60) 


This relation is verified by (1.15) and (1.57). It shows that if an optimum 
solution of the primal has a component x; > 0, then the corresponding constraint 
of the dual is satisfied as an equality in its optimum solution. It shows further 
that the variable of the primal x;, which corresponds to a constraint of the dual 
satisfied in its optimum solution as a strict inequality, must vanish in the optimum 
solution. 

The conditions (1.60) are thus necessary for optimum solutions. They are 
also sufficient, since, by (1.53), they imply that the w and z, corresponding to 
solution of the dual constraints and a feasible solution of the primal constraints, 
ure equal to one another and so by (1.54, 59) equal to their common optimum 
values, This result can also be expressed as the following theorem. 
The necessary and sufficient condition that a feasible solution [x;] shall be 
optimum for the primal (1.11—13), which is known to have an optimum, is that 


there should be a solution {uj} to the dual constraints (1.51), which satisfies 
these constraints as an equality for all values of i for which x; > 0. (1.61) 


The dual problem was introduced in (1.50, 51) by definition. It can also be 
approached by introducing ‘Lagrangian multipliers’ {uj} in an attempt to solve 
the primal problem. The ‘Lagrangian’ Z corresponding to (1.11, 13) is 


Pp n P 
Z=)> mt > wy (41- DY ayxi) 
i=1 i=1 i=1 


n a, n 
=> buy+> (a ps ayy) Xi (1.62) 
jm (1 


j=1 
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and the necessary conditions that this function of non-negative [x;] should 
have a minimum are that (1.51) should be valid, since otherwise x; > % (any i 
with (1.51) false) implies Z > —°°, and that 


wuyj=c; (x; > 0). (1.63) 


W DAs 
a 


These are in fact the conditions of (1.61), which are thus expressed in a 
Lagrangian form. 

If Z is now thought of as function of non-negative [x;] and {uj}, then the 
condition that it is stationary with respect to {u;} is the constraint (1.13). If 
this last is imposed upon the problem ‘make Z stationary’, Z reduces to z and 
the primal problem is obtained, which is taken here to be a minimum problem. 
If on the other hand, the conditions of (1.51) and (1.63) are imposed, then Z 
reduces to w a function of {uj} and (1.51, 63) with [x;] eliminated, reduce to 
(1.51). This gives the dual problem, which is known to be a maximum problem. 
The primal and dual problems can thus be generated in a reciprocal way from 
the necessary conditions that the Lagrangian Z should be stationary with regard 
to {u;} and a minimum with respect to [x;]. 


1.5 The dual problem for structures of least volume of material 


The primal problem for structures of least volume of material is given by 
(1.8—10). Comparison with (1.11—13) gives 


[x;] = [isi Tn Titel wills 
foh-={ijetilas leeolehios 


(1.64) 
lay] = [opKyj/o,.. op Kp j/0, 0c Kyj/6, - --0¢K pj/0] , 
[5] = (Fi). 
The dual problem can thus be written using (1.50, 51) as 
n 
max W= } Fjuj/oe, (1.65) 


g=1 


subject to, 


n 
Kyuj <(celop)li, — SY Kizuj <(oeloc)l; (@i=1,2..m), (1.66) 
1 j=1 


n 
j= 
where the dual variables are u;/oe, with o given by (1.7) and € a positive 
infinitesimal. 

Let {u; }be interpreted as virtual displacement components at the nodes 


corresponding to the given forces [/j]. The objective function W is then the 
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virtual work of the given forces taken over the virtual displacements, divided by 
the constant oe. Equation (1.1) gives for any {u;} 


m n n 
y ( > Kyu) T= 3 By, (1.67) 
i=1 \i=1 g=u 

which is the principle of virtual work. The virtual strains {ei} in the members of 


the framework are thus given by 
n 
= > Kyuj/l; @=1,2..m), (1.68) 
j=1 


and so (1.66) can be written as 
€;<o¢e/op, —€;<€/o¢ (@=1,2..m), (1.69) 


which implies that the virtual strains must lie between —oe€/o¢ and aeé/ap. The 
dual problem thus seeks to maximize the virtual work W, subject to bounds 
placed on the values of virtual strains as given by (1.69). 

The theorem of (1.61) leads to the following conditions for optimum 
frameworks: 


The necessary and sufficient conditions that a pin-jointed framework, selected 
from a given framework, should have a minimum volume of material are that 

it should be capable of carrying its given forces with stresses Oy in its tension 
members and —dg in its compression members and should allow a virtual dis- 
placement of its nodes, which produces a strain of o€/0-p in its tension members 
and a strain of —o€/d¢ in its compression members and strains not outside this 
range in members of the given framework, which are not present in the 

optimum design. (1.70) 


The primal variables 7;', 7;", for any i, cannot both be positive, since this would 
imply contradictory equalities in (1.69). One or both of 7;’, T;", for each i, must 
thus be zero, If 7;' > 0, 7;" = 0, then by (1.6) 7; > 0 and 7;/A;= op. Also by 
(1,69) e; = o¢/oy. If T = 0, 7;" > 0, then by (1.6) 7; <0 and T;/A; = —o¢. Also 
by (1,69) e; =—-oe/o;. If T; = T;" = 0, then by (1.6) T; = 0 and A; = 0. Also 
(1.69) must be satisfied. 

Conditions of the type found in (1.70) were first given in Michell (1904) and 
are called Michell conditions in his honour. Michell considered sufficient condi- 
tions only and dealt with ‘continuous’ strain fields and structures. An account 
of his theory is to be found in Chapter 4. 


1.6 Direct solution of problems using the dual 


The dual problem of (1.65, 66) can be transformed to the standard form of 
(1.1113) by introducing non-negative slack variables in (1.66) and replacing 
{u,} by 

wtuy-u"; uj,uy 20  (f=1,2..n). (1.71) 
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Relations (1.65, 66) then yield 
n 
max W= > Fj (uj—u;)/oe, (1.72) 
j=1 
subject to 


n n 
Y Kyu —uj)+ i= (elon), — Y Kyu uj) + ai = (oeloc)i 
j=1 j=1 


Gaeeein)s 10:73) 
and 
ui,uj >0 (f=1,2..); py qi 20 G@=1,2..m). (1.74) 


Equations (1.73) may also be written 


n 
Y Kyu; —uj) + pi=(Gelog)li, py + qi = el; (1/op + log) = 1,2..m) 
at (1.75) 


which shows that {Pi} can be treated as ‘bounded variables’. This means that 
analysis by the simplex method can be somewhat shortened (Dantzig 1963, 
Chap. 18). 

Equations corresponding to (1.14) can be obtained from (1.73) in the form 


Dj = (Ge/op)l; — » Ki(uj — uj) 
4 (i=1,2..m), (1.76) 

qi = (Geog) + S Ky(uj—uj) 
i=1 


which give a basic feasible solution, from which the simplex method can begin. 

The final optimum solution or any basic solution for that matter, cannot 
contain both wu; and uj, for any j, among the basic variables, since (1.73) could 
not be solved for such a pair of variables among others, because the corresponding 
determinant vanishes. The apparent indeterminancy of uj, uj in (1.71) thus 
disappears for basic solutions. 

The optimum solution will determine the virtual displacements {uj} and also 
give those members which have p; = 0 or q; = 0, i.e. those which reach the strain 
limits o€/o7 and —e/0g. This gives a layout of tension and compression 
members from which the optimum structure can be found. In many cases more 
potential members will be supplied than are required to form an optimum 
structure. This surplus may be eliminated by imposing the static conditions 
(1.1), but sometimes a redundant structure will remain. This defines a class of 
optimum structures with the same volume of material, which is more general 
than could be obtained from a basic solution of the primal. 
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In the special case when 
Op =0g =O (1.77) 
the solution to the dual shows that, for the optimum framework, members in 


tension with stress o will have a virtual strain € and members in compression 
with stress —o will have a virtual strain —e. It follows that if 


e=9/E, (1.78) 


where £ is Young’s modulus for the material of the framework, then the virtual 
strains coincide with the real strains in the members and the conditions of 
compatibility of strain are satisfied. The optimum is thus an optimum elastic 
design, even though the optimum framework is a redundant structure. 


1.7 A further example 


Consider the problem defined by Fig. 1.2. This shows the geometry of a given 
layout for a framework, which is taken as symmetrical about the line through 


K-—--—--[--- >< ---l--- > 
uy, 5 act 


Fig. 1.2 


the nodes where the given forces F are acting. The enumeration of the members 
is shown as are the components of virtual displacements {u}. The virtual strains 
{ei} in the members are given by 


@) = uUy/l, €y =uU3/2l, €3 = (u, —uz)/21, €4 = (U3 — U2 fl, 


5 © (Ug Uy )/1, €6 = (Ug ~My — Uz )/21, €7 = (Ug — Us + U3)/21, (1.79) 
€g = (Ug ~ Us D/L. 
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Assuming for simplicity that op = og = 0, the constraints of (1.69) can be 
written 
€;<e, —€;<e @=1,2..8) (1.80) 


and so, by (1.79) introducing slack variablest { p;} and { q;}, yield 
u,t+p,=el, —u, +q,=4l, 
uz + pz =2el, —u3 + q2 = 2el, 
Uy — U2 + pz =2el, uy + uz +q3 =2el, 


Uz —U2 + Pg =el, U3 + U2 + qq = el, 


Ug—U, +ps=el, Ug tu, +qs =el, sedis 
Uy — U2 + pe =2el, Ug ty t+U2 + Ge = 2el, 
Ug —Us tu3 + pz = 2el, —Ug + Us — U3 + G7 = 2el, 
Ug —Ust+pe=el, Ue tus + qe = el, 
with 
p,20, gq;20 (=1,2..8). (1.82) 


The sixteen equations (1.81) can be solved for u; (j = 1, 2. . 6) and for ten of 

the slack variables. It is clear in all cases which sign the virtual strain of magnitude 
e will take, if the corresponding member is present in the optimum structure. 
This suggests that p1, P2, 73, Ps, Po and pg could be taken as non-basic variables. 
Following this suggestion leads to 


y= el'— pis Ga— 21 =p, 

uz =2el—p2, G2 =4el—pr, 

uz =3el— pi — qs, P3=4el—qs, 

Pa = 2el—p, + P2— 43, a =P1—P2 + Qs, 
ug = 2el—p,— Ps, qs = 2el—Ds, 

Us = 6el—2p1 —G3—Po, Io = 4el — Deo, 

us = Sel— 2p: —43— Pot Ps, 4a = 2el —Ps, 
Pz = 3el—pi +p2—q3+tPs—De+t Pa. 

47 =el+P1—P2+43—Pst+Po—Ps, 


(1.83) 


which defines a basic feasible solution. There is no need to introduce {uj} and 
{uj} from (1.71) here, since the determination of the maximum depends, thanks 
to (1.83), on the slack variables only.+ 


{ It is to be remarked that the quantities differ from those of (1.73) by numerical 
multipliers. 


+ It is always worthwhile to solve as in (1.'83) in simple problems. The result may not 
always lead to a feasible solution however and certain simplex iterations, involving the 
slack variables, may be required to achieve feasibility, before the process of maximization 
can begin, 
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For the present problem W of (1.65) is given by 
W = F(uz3 + ue)/o€ = F(8el — 2p; — P2 — 43 — Po) /0€, (1.84) 
where (1.83) has been used. Thanks to (1.82) the maximum follows without 
further calculation. It is 
Wax = 8Fi/o, (1.85) 
with 
Pi = P2 = 43 = Po = 9. (1.86) 
Equations (1.83, 86) now give 
u, =el, ug = 3el, uz = 2el, ug = 2el—Ps, Us = Sel + pg, Ue = 6el, (1.87) 


and 
Pp, =0, qi =2el; p2=0, qz=4el; p3=4el, q3=0; pa =2el, a4=0; 


ds =2el—ps: Po=0, do=4el; p7=3el+ps+Ps, q7=el—Ps—ps; } (1.88) 
da = 2el— Dg, 


which by (1.82) imply 
ps20, pg20, Ps tpg <el. (1.89) 


This solution is indeterminate. Basic solutions can be obtained by writing 
(i) Ps = Ps =0, (ii) ps = 0, pg = el and hence q7 =0,or (iii) ps = el, pg = 0 and 
hence q7 = 0. The resulting solutions are degenerate and give the following 
layouts: 
(i) members 1, 2, 5, 6, and 8 in tension, members 3 and 4 in compression, 
(ii) members 1, 2, 5, and 6 in tension, members 3, 4, and 7 in compression, 
and 
(iii) members 1, 2, 6, and 8 in tension, members 3, 4, and 7 in compression. 
Non-degenerate basic solutions can be obtained by writing (iv) ps = 0,0 <pg <el, 
(v) ps = 0,0 <ps <eland (vi) ps + pg = el, O<ps <el and hence q7 = 0. The 
corresponding layouts are: 
(iv) members 1, 2, 5, and 6 in tension, members 3 and 4 in compression, 
(v) members 1, 2, 6, and 8 in tension, members 3 and 4 in compression, and 
(vi) members 1, 2, and 6 in tension, members 3, 4, and 7 in compression. 


The final type of solution can be obtained by writing (vii) ps > 0, ps > 0, 
Ps + Py <el. This gives the layout: 
(vii) members 1, 2, and 6 in tension, members 3 and 4 in compression. 


Consideration of the equilibrium equations corresponding to u4 and us 
removes members 5, 7, and 8, together with their common node, from the 
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layouts (i) to (vi) and reduces them all to layout (vii). The principle of virtual 
work gives for the remaining members, using (1.79), 


Tyuy + T2u3//2 + T3(uy — U2) 2 + Ts (uz — U2) + To(ue — U1 — U2 )/V2 (1.90) 
= F(u3 + ue)/2, 
and so the remaining equilibrium equations are 
T, + T3//2 - To//2 = 0, 
-T3/\/2 — Ts — Ts//2 = 0, 


1.91 
T,//2 + Ts = F/2, ( ) 
Te/V2 = F/2, 
which can be solved to give 


T, =F+Ts, T, =F//2-V2Ts, T3= -Fh/2-/2Ts, Ts = F/V/2. (1.92) 


Finally, the specification of layout (vii) implies 
—-F/2 <T, <0. (1.93) 


The optimum framework is thus redundant, but the value of the redundant 
load T, must satisfy (1.93), so that the Michell conditions of (1.70) are satisfied. 
The extremes of (1.93) give T, = 0 and T; = —F/2 or T = 0. The simplest layouts 
thus have members 1, 2, and 6 in tension and either member 3 or member 4 in 
compression. The combined layout gives an optimum structure, so long as both 
members 3 and 4 carry compressive loads. The volume of material for all these 
optimum structures is given by (1.85). All the designs are both optimum elastic 
and optimum plastic frameworks. 

This example is interesting in that it shows that any basic solution (i) to (vi), 
which might be found by a purely numerical application of the simplex method,+ 
gives the redundant optimum structure. This would not be obtained by a 
similar analysis of (1.8—10). 


1.8 Approximation to absolute optima 


The methods of analysis considered in this chapter are confined by their nature 
to the selection of optimum frameworks from a framework of given layout. A 
search for absolute optima requires consideration of all possible layouts and 
hence a selection from an infinite set of possible frameworks. An approximation 
to this ideal can be achieved by covering the region of space in which the 
required structure is allowed to lie by, for example, a closely spaced rectangular 
grid of nodal points and allowing members to lie along all segments joining nodal 


+ For example by the use of a standard programme on the digital computer. 
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points. The nodal points must include the points of application of given forces 
or of discrete approximations to given continuously distributed forces and all 
points of support. 

Let (x;, ¥;, Z;) be the coordinates of the nodal points, /;; the length of the 


segment (i, j) given by 

lj ={@j — xi) + Oi yd? + | - 2}, (1.94) 
Aj the cross-sectional area of the member (i, /), Tj the end load carried by the 
member (i, /), (X;, Y;, Z;) the external force at node i, o = op = dg the allow- 
able stress and V the volume of material. The primal problem can be formulated 
as: 


min,V=" >. TigAy: (1.95) 
Pairs i, j 
if 

subject to 

YX TyXj- Xin Yj —Ie 4-2 lg + Xp Ye Z)=0 (alli), (1.96) 
isi 
and 

|Tj|\<oAg (all pairs i, j). (1.97) 


This formulation can now be converted to that of (1.1—1.4) by an appropriate 
change of notation and the omission of equations involving forces of constraint. 

Let (u;, v;, w;) be the virtual displacement at node i, € the limiting strain and 
W the virtual work divided by ce. The dual problem can be formulated as: 


max W= > (X,u; + Yin; + Zw,)/o€, (1.98) 


subject to 
(uy ~ uy — x;) + & — 20; —¥) + O4j —w) EG -z)i<e§ @lli/). (1.99) 


This formulation can now be converted to that of (1.65, 66) by an appropriate 
change of notation. 

Calculations of the size envisaged here can only be carried out using digital 
computers, Many standard programmes exist for the solution of the linear 
programming problems which arise here. The compilation by Kienzi et al. (1968) 
has already been mentioned. A programme which uses this dual formulation 
and which is specially adapted to the solution of framework optimization 
problems is given in Hemp and Chan (1966). Results obtained using this 
programme are given below. 

The problem of designing an optimum structure to transmit a force F to 
two fixed supports is illustrated in Fig. 1.3.4 The structure is assumed to be 
symmetrical about the horizontal line through the point of application of F. 


+ Taken from Hemp and Chan (1966). 


20 Pin-jointed frameworks 


The assumed grid of nodal points is shown and these are taken to have a 
completely general virtual motion consistent with symmetry. The solution to 
the dual problem of (1.98, 99) gives strains + € in the members shown in the 
upper-half of the diagram. The full lines have strain +e and are potential tension 
members, the dotted lines have strain —e and are potential compression 
members. Several of the members are seen to be unnecessary to carry the load 
F, but those that remain make up a structure with order of redundancy three. 


A statically determinate optimum structure is shown in the lower half of Fig. 1.3. 


The volume of material in all these optimum structures is given by 


Vinin = 4-83 Flo. (1.100) 


The same problem as that of Fig. 1.3 is illustrated in Fig. 1.4, but here a 
finer mesh of nodal points is used. The upper half of the diagram shows all the 


Fig. 1.4 
+ Taken from Hemp (1964). 
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potential members and the lower half a statically determinate optimum 
structure, The volume of material required at the optimum is: 


Vinin = 4-61 Fl/o. (1.101) 


It is interesting to compare (1.100, 101) with the exact solution obtained in 
section 4.9 below, which gives 


Venin = 4-34 . . Fife. (1.102) 


The best approximation is thus some six per cent above this exact value. 
Another example of a computer designed structure is shown in Fig. 1.5.+ 


Fh] fy : Fl; 


Fig. 1.5 


Here the problem is to design an optimum structure to transmit a load F, 
uniformly distributed over a span L, to two vertical supports at the end of the 
span. The distributed load is represented by fifteen concentrated loads as 
shown. The resulting structure of arch, hangers and tie has a volume given by 


Venin = 1-067 . . FL/o. (1.103) 


1.9 Alternative loading conditions 


Practical structures are usually designed to withstand more than one system of 
given forces. The present theory must be extended to deal with this more 
complex problem. 

A single system of given forces has been denoted by [F;] . Several systems 
can be written [F,;], where the subscript @ picks out a particular system by 
taking values 1, 2, .. s. New quantities [7;], [T2:] , [Tj] can be introduced 
in the same way to denote the corresponding end loads in members and slack 


| This was obtained by Dr. W. J. Supple while working at Oxford under contract to the 
Ministry of Technology (1968-70). 
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variables. Generalizations of (1.1) and (1.6) will be valid for the new variables, 
with of course [A;] unchanged, since there is only one structure to be 
designed. It is convenient in the present case to retain the variables [A], rather 
than to eliminate them as in (1.8—10), and to express the primal problem in the 
form: 


m 
min V= > Ajl;, (1.104) 
i=l 
subject to 
m 
¥ Ky(orTai — %cTai lO = Fog = (G= 1,2... a= 1,2. - 8), (1.105) 
i=1 
Ti 4T=0A; (G=1;2:.m; @=1;23.9); (1.106) 
and 
TNA S02 LG=1 2) sm ae Des), (1.107) 


This is a standard problem of linear programming of the same form as (1.11—13) 
and can be solved using the simplex method. It leads to larger calculations than 
(1.8—10) for a single system of forces. 

The dual can be formulated using (1.50, 51), but as an illustration of the 
alternative Lagrangian approach of (1.62), the Lagrangian V * corresponding to 
(1.104—106) will be written down and the dual deduced directly from it. 
Introducing Lagrangian multipliers [u,;] /o€ and [l;¥q:]/oe where [u,;] are 
virtual displacements then gives 


m s n m 
V*= > Adit YS YD Cailoe) | Fay - Y Kig OrTai — GcTai)/ 0} + 

i=1 a=1 j=1 i=1 

(1.108) 
= beng , ” 
+ > DY Citaile) {Tai + Tai — 8Ai}- 
a=1 i=1 

The conditions that V* should have a minimum for the variables [A], [Tz] 
and [7,;] are then 


s s 
(= y Tat) Ai=0, €> 3 Ya =1,2..m), (1.109) 
a=1 a=1 
ai - %r€ail9) Tai =0, Yai Fr Eai/O (@=1,2..m; a=1,2..5), 
(1.110) 
and 
(Yai + Fc €ai/6) Tai =, Yai>—Sc€aild (i= 1,2..m; a=1,2..5), 
(aly 
where {€;} are member virtual strains given by 


€ai= > Kyttajll; C12 emp a ©1572 25.8): (1.112) 
j=l 
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The equations in (1.109—111) are another way of writing relations of type (1.63), 
while the inequations express conditions, which prevent V * +—=, as one of the 
non-negative variables A;, T,;, Ta; tends to +-°. The dual problem, which is a 
maximum problem, is obtained by imposing (1.109—111) on V* to give the 
objective function W, which is seen to depend upon the variables [uj] alone, 

and by eliminating [A;], [Tz;] and [Tz;] from (1.109—111) to give the dual 
constraints. The dual problem can thus be written 


s n 
maxW= > > Fajllaj/oe, (1.113) 
a@=1 j=1 
subject to 
s 
€> > Yais Yai > 97 €ailO, Yai >—OcEai/ (i= 1,2..m;, a=1,2..5). 
ei (1.114) 
The function W of (1.113) is the sum of the virtual work of the various 
systems of given forces taken over their corresponding displacements. The 
quantities [y,;] can be eliminated from (1.114) to give 


s 
Y Créaile or —oc€gi/o)<e (i=1,2..m), (1.115) 
a=1 
where all the 2° alternative relations are implied. All these can be obtained from 
(1.114) by adding an appropriate selection of the second and third relations of 
(1.114) to the first. Conversely (1.114) can be obtained from (1.115) by 
writing? 

Yoi = Max (Op Egi/G, —Ge€qi/o) @=1,2..m; a@=1,2..5), (1.116) 
since, for each i, one of the 2° relations (1.115) will give the first of (1.114) and 
(1.116) satisfies the second and third of (1.114). The relations (1.115) can thus 
be taken as the constraints for the dual. They are 2*m in number and show how 
rapidly the size of the present problem grows with the number s of alternative 
systems of given forces. 

The relation (1.115) can also be written 


DY opleail/ot+ Y oelegillo<e  (=1,2..m), (1.117) 
€ai > 0 €ai <0 
or, if op = dg = 0, 
S | eas | <e (@i=1,2..m), (1.118) 
a=1 


which shows the constraints imposed on the virtual strains in a particularly clear 
form. 


t The choice of yay, for given é, is unique if (1.115) is satisfied as an equality; otherwise 
it is indeterminate. 
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The relations between the optimum solutions of the primal and the dual 
enable a number of deductions to be made about the optimum structure. The 
following cases can arise: 


(i) TZ; >0, Ty’; > 0, which by (1.106) implies A; > 0 and gives a stress 
between o7 and —d¢.7 Also since the corresponding dual constraints 
must be equalities (1.114) gives Yq; = O7 €ai/0 = —O¢ €gi/o and so 
€ai ~ Yai = 0. 

(ii) Tz; > 0, Tz; = 0, which implies A; > 0 and gives a stress o7.+ Also 
Yai = 7 Egil O > —O¢ Eqi/G and so Eg; > O and yq; =O. 

(iii) Tz; = 0, Tz > 0, which implies A; > 0 and gives a stress —0¢.7 Also 
Yai = —8¢ qi /O > O7 Eqi/o and so €g; <0 and Yq; > 0. 

(iv) To; = Te; = 0, which implies A; = 0. Also Yq; > 0 and —o7q;/0¢ < €qi 
<0Yqi/Or- 


s 

The case A; >Ogives } Yq; = € and so it follows that every member that 

a=1 
is present must be stressed to the limit by at least one system of given forces. 
A member that is stressed to the limit in tension will have a corresponding 
positive virtual strain, a member that is stressed to the limit in compression 
will have a corresponding negative virtual strain and a member that is not 
stressed to the limit will have a corresponding zero virtual strain. A member 
that is absent can have a strain of any sign. A member that is present will give 
an equality in (1.117) and so its virtual strains must satisfy 


> or leail/o + dy %elealio=e G=1,2..m; A;>0). 
Stes = ox Siena ot (1.119) 
For a member with A; = 0, the relation (1.117) must be satisfied by its virtual 
strains. 
The generalization of the Michell type theorem of (1.70) is as follows: 


The necessary and sufficient conditions that a pin-jointed framework, 
selected from a given framework, should have a minimum volume of material 
are (i) that it should be capable of safely carrying its alternative systems of 
forces, (ii) that each member should be stressed to the limit oy or —0¢ by at 
least one system of force, (iii) that, corresponding to each system of force, the 
framework allows a virtual displacement of its nodes, which produces non- 
negative strains in members loaded to stress O7, non-positive strains in members 
loaded to a stress —Gg and zero strains in members not stressed to these limits, 
(iv) that the strains of (iii) satisfy equation (1.119) and (v) that the strains in 
all members (including those not present in the optimum structure), correspond- 
ing to the displacements of (iii), satisfy the relation (1.117). (1.120) 


The general theory developed above is complete, but clearly likely to lead to 
quite large calculations, compared with those for design to a single system of 


+ This follows from Tq /Aj * (op Tea ~ %¢ Text (Tea + Tai) 
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given forces. In the special case of two systems of given forces (s = 2), with 
Oy = Gg = G, considerable simplification is possible. This is due to the fact that 
the problems can, in this special case, be reduced to two problems of optimum 
design with a single given system of forces for each. 

The relations (1.115) with s = 2 and o7 = og = 0 give 


€1; + €a3 SE, 4; — €ai SE, Eq; + €gi SE, —GqQi—Egi SE = (@= 1,2 .. m), 


(1.121) 
which suggest the change of notation 
U4; = uyj + Ug;, Uo; = uy; — U9; G=1,2..n), 
ak i la ct ns ad aa } (1.122) 
Eyz= €1; + €2i, Eoi=€4i—€2i = @=1,2..m). 
Equations (1.112, 122) give 
n 
Eqi= > Ki Uajll; (@i=1,2..m; @=1,2), (1.123) 
i=1 
while (1.121) can be written 
Eqi<é€, —Eqi<e€ (i=1,2..m; a=1,2), (1.124) 
and (1.113) as 
2 
maxW= ¥ Wy, (1.125) 
a=1 
where 
n 
Wa= > fog Uaioe  (@=1,2), (1.126) 
j=1 
and 


fay = (Fay + Foj)/2, fo; = (Fay — Fopl2 (7 =1,2..m). (1.127) 


The present problem thus splits into two, namely to maximize W, (a= 1, 2) of 
(1.126) subject to the constraints of (1.123, 124). The given forces for these 
problems are, by (1.127), equal to half the sum and to half the difference of the 
two systems of given forces. 

The theory of the present section, like that of section 1.1, has been developed 
on the basis of plastic design. The conditions for optimum structures that have 
been found do not yield optimum elastic designs, except in trivial cases, such as 
when each system of given forces is carried by separate non-interacting 
frameworks. 
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1.10 An example with two alternative loading conditions 


Consider the given layout of Fig. 1.2 and let two given systems of force be 
specified by 


Ryjy=0 (G=1,2..5), Fxg = 2F, 
a G ) Be } (1.128) 
Fy,=0 (f=1,2,4,5,6), Fo3 =2F. 
Equations (1.127, 128) give 
;=0 (=1,2,4,5), fis=he=F. 
fi G ). fis =fie } (1.129) 
foj=90 (=1,2,4,5), fos =F fog =F 


The system [f,;] defines the problem solved in section 1.7 and so by (1.85, 87, 
89), with an appropriate change of notation, a basic solution to this first 
problem can be written: 


(Wy )nax = 8Fi/o, 


Fy be 
Uy; = el (1, 3, 2,2, 5, 6). ames 


A solution to the second problem with forces [f2;] can be obtained by 
solving (1.81) in the form 
u, =el—Py, qi =2el—p, 
u3 =—2el + G2, P2 =4el—Qr, 
uz = —el + q2— 4a, Pa =2el— Qa, 
P3=Pi+q2—44, 43 =4el—pi —42 +44, 
ug = 2el— py — Ps, Qs = 2el—Ps, 
Ug = 2el—P1 + Q2—G4— Po, Wo =4el— Do. 
us =el— Pi + G2—Ga—Pot Ps, 93 = 2€l— Ds, 
Pz = 3€l—q4 + Ps— Pet Ps. 47 =€l + Gs— Ps + Po— Ps, 
which gives for W the expression 
W = F(ug — u3)/oe = F(4el — p, — G4 — Do )/a€. (1.132) 


The maximum is 4F1/o with p, = q4 = ps = 0 and so a basic solution to the 
second problem can be written, on changing to the notation of (1.130) and 
using (1.131), in the form: 


(1.131) 


(Wo) max = 4Fl/o, 
Ug; = el (1, —1, —2, 2, 1, 2). (1.133) 
The total volume of material required for the present problem of alternative 
loading is by (1.125, 130, 133) 


Wmax = 12Fl/o. (1.134) 
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The virtual displacements follow from (1.122, 130, 133) as 


uy; = el (1, 1, 0, 2, 3, 4), 
(1.135) 
Ug; = el (0, 2, 2, 0, 2, 2), 
and so by (1.79) the virtual strains can be written 
€3; =€(1, 0, 0,—1, 1, 1, -1/2, 1), 
li ( I ) (1.136) 


€o; = €(0, 1,—1, 0, 0, 0, 0, 0). 


Inspection of (1.136) shows that (1.118) is satisfied as an equality for all i, 
except i= 7. This means that A, = 0 and hence that 7,7 = 727 = 0. Equilibrium 
then shows that 7,5 = Tys = 0 and 75 = 72s = 0. It thus follows that 
As = Ag = Oand so the three members i = 5, 7, 8 can be omitted. 

Equation (1.116) shows that here Yq; = |€g;|(@ = 1, 2; # 7) and so 
Yai — €ai > 0, if Eg; <O and Yq; + €g; > 0 if €,; > 0. In the first case the second 
of (1.114) is satisfied as a strict inequality and the corresponding primal variable 
T,,; = 0. In the second case the same is true of the third of (1.114) and so 


T,,, = 0. Application of these results to (1.136) gives 
Ti, = Tia = Tie =0, 
os 3 16 a .137) 
To2 = To3 = 0. 


The equilibrium equations (1.91) adapted to the present cases give, writing 
T, = T; —T," , using (1.137), and changing to the present notation 
Thi + (Ti3 — T1'3)/V2 — Tie /V2 = 0, 
—(TMi3— Ti3)/V2 + Tis — Te //2= 0, 


1.138 
(Ti = Ti Wv2 = 14 =O c ) 
TislV2 =i 
and 
Tx, — Tz — T2332 - (Tre — Tx6)/V2 = 0, 
Tz3/\/2 — (Tza — Tz) ~ (Tae — Tre )/V2 = 0, 1.139) 


Tz2//2 + (Tia — Tz) = F, 
(Tze — Tx )/V2 = 0. 
Equations (1.106) give, using (1.137), 
0A, =Ti1 =T3, + Tr}, 
0A, = Ti2 + Ty = Tra, 
0A3 =TM%3 +73 = Tr, (1.140) 
0Ag = Ty's = Tr4 + Ty, 
0A 6 = Ti6 = Tr6 + Tre. 
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Equations (1.137—140), together with the non-negativity conditions (1.107), can 
be solved to give 
Ty, =2F- Tia, Tri =O, 
Tin =V2Ti5, T; = 0, 
T;3 = 0, Ty's =V2(F-T13). (1.141) 
Ti, =0, F/2<T4 <F, 
Tis = V2F, Ty¢ = 0, 
Tz, = 3F/2-T4, Tz: = F/2, 
Tx, = V2T;'5, T, = 0, 


Tz3 =0, Tr =V/2(F-T5). (1.142) 
Tra = F/2, Tr = Tra =F; 
Tie = F//2, Tre = F/V2, 


and 


A, =(2F-Tj)/o, A2=V2T3/0, A3=JV2(F-Tis)/o. 
Ag=Tyslo, Ag = V2F/o. 
The end loads in the members for the two loading cases follow from (1.141, 
142) as 
Ty =2F-Tia, Ty =V2T3, Ti3= -J/2(F - Ty). Ti4=-Ti4, Tis = - 
Ty =F-Th, Ba=V2Th, Ts =-V20—Tis), Tra =F-Tye, Tas =O 
(1.144) 


‘| (1.143) 


where, by (1.141), 
F/2 <4 <F. (1.145) 

The optimum solution obtained above is a redundant structure, except when 
T;4 = F and T,3 = T,3 = A3 = 0. It gives a range of designs depending on the 
parameter 7;4. In the first system of loading all members are stressed to the 
limits +o. In the second system members 2 and 3 are stressed to the limits, 
while members 1, 4 and 6 are understressed. It may be verified that (1.143) 
gives a V equal to W,,,, of (1.134) and that all the conditions of (1.120) are 
satisfied. 


1.11 Constrained minima of functions of non-negative variables 
Consider the problem 
min z = f(X1, X2 --Xni Vir V2 --Ym)> (1.146) 
subject to 
Oj 0X1, X2 --Xni Vir V2 --Ym) = 9 (j=1,2..p; p<m+n), (1.147) 


and 
Ve 20 (kK=1,2..m). (1.148) 
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Assume that values [x;] , [y,] are known that make z of (1.146) a minimum.+ 
It then follows, for all small variations [5x;] and [5y;], which by (1.147, 148) 
must satisfy 


pa Vix; ant 5 Giy, Se =0 (f=1,2..p), (1.149) 
and 
by, =0 (all k such that y, = 0), (1.150) 
that 
n m 
D ke, dxi+ D fy He FO. (1.151) 
i=1 k=1 


Only first order terms have been written in (1.149—151), since the arguments 
to be developed depend upon the signs of certain quantities, which are not 
affected, for small enough variations, by higher order terms. 

Introducing Lagrangian multipliers [A;] leads, using (1.149, 151) to 


n 


Pp m P 
v3 (f+ D) iO; ) Bx: + y (+> i Ory )OM 20, (1.152) 
j=1 k=1 j=1 


i=1 
which is valid for all [5x;], [5], which satisfy (1.149, 150). 
Assume now that (1.149) can be solved for 5x;(i =i), i2 - -i,), 

by, (k = Ki, ke ..k_), where s + t =p and the selection from [5y;,] is such that 
Yr > 0 (k =ky, kz . . ke). (1.153) 

This means that the matrix 

[Gix;5 iy, | has rank p (Gj=1,2..p; i=1,2..n; all k such that », > 0), 

(1.154) 

which is a ‘regularity condition’ for the present problem. The remaining 

variables in (1.149) are denoted by 5x;(i=i,41,--in), 5¥e(K =Keo1.--km) 

and these can be regarded as independent variables, taking any values, subject 


only to (1.150). 
The multipliers are now chosen to satisfy the equations 


fer + Fy Gx, =0 (i=i;, 2. <i); 

- (1.155) 
Syn +2 Ny Giy, =O — (K=Ki, ka. - er), 

“1 


| The functions f and [¢] are assumed to be continuously differentiable in a neighbourhood 
of these values. A variable as a subscript is used to denote a derivative. 
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which is always possible, since the determinant of (1.155) is the same deter- 
minant as that assumed to be non-vanishing in the solution of (1.149). 
Substituting from (1.155) into (1.152) then gives 


in Pp Rm P 
s (%: +E Gin, )ox: eS (Sn +> % tion) 9% >0, 
i=i,41 j=l k=Resy ie (1.156) 
where the variations 5x;(i = i,41,-- in), OV_ (K = Krai, -- kn) are independent 
variables. 

Now (1.156) implies 


(tit EW O0,)8 20 O=feeay in) 
j=1 
P (1.157) 
(n+ A %)8m 20 Keer.) 
{=1 


since it is possible to set all but one of the variations in (1.156) equal to zero 
and leave the remaining one arbitrary. It then follows that 


P 
fet EO Gej=0 Cheer. -tns 
s=1 


5 (1.158) 
Sut DM by, =9 (k =Kp+1,--Km and k such that », > 0), 
i=1 


P 
fin t EWN Giyp2O (R= Keaa.-- hm and k such that y, = 0), 
io (1.159) 


since if one of the equations in (1.158) is not valid, then, by giving the 
corresponding 5x; or 5y;, a value of appropriate sign in (1.157), a contradiction 
can be obtained and if (1.159) is false, then the corresponding inequation of 
(1.157) can be contradicted by giving the 5y,,, which occurs in it, a positive 
value, which (1.150) allows. 

A combination of (1.155, 158, 159) gives as necessary conditions for the 
minimum of z of (1.146), subject to (1.147, 148), the relations 


fz, s 4 dy, =0 G@=1,2..n), (1.160) 
j=1 
and 


P P 
(f+ > Ay O1yy) Me =O. font DN Fin, 20 (K=1,2..m). 
hee oe (1.161) 
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The relations (1.147, 160, 161) give p + m + m equations for the p +1 +m 
variables [Aj], [x;], and [),]. The solution which gives minimum z, if it exists, 
is to be found among the solution to these equations and must also satisfy the 
inequations of (1.148, 161). It is not possible to say in advance which values of 
k give y, = 0 and so the various alternatives have to be examined. It may be that 
(1.154) cannot be verified in advance and must be checked after a solution has 
been obtained. The process of solving (1.160, 161) for A; usually provides a 
proof of (1.154). Irregular or degenerate solutions for which (1.154) is not 
satisfied are possible. These can be avoided by making suitable small changes 
in the constraints of (1.147, 148). 

The conditions of (1.160, 161) can also be expressed by the statement that 
they are the conditions that the Lagrangian 


Z=f+ FG; (1.162) 
j=1 


regarded as a function of [x;], [y,], with [Aj] as parameters, should have a 
minimum, subject to (1.148). 


1.12 Structures of maximum stiffness 


Consider, as in section 1.1, a framework of given layout, which is defined by 
the areas of cross-section [A;] of its members. Replace (1.2) by the condition 
that no A; can be less than a positive limit Ag. This can be written 


A;=AotpP;, pj20 (i=1,2..m). (1.163) 


The relation (1.163) is introduced so as to avoid degeneracy in the minimum 
problem to be considered. The limit Ao > 0 will be imposed at the end of the 
analysis. 

Let the framework be loaded by forces [F;] and let the resulting displace- 
ments, corresponding to an elastic deformation, be denoted by {y} to distinguish 
them from the virtual displacements {u;} of section 1.5, with which they are of 
course kinematically identical. Equation (1.68) gives for the strains {e;}in the 
members 


E=> Kyr/li (@=1,2.. m). (1.164) 
i=1 
The end loads [7] follow from 
T; = EAjé; (=1,2..m), (1.165) 


where £ is Young’s modulus for the material of construction. Equations (1.1, 
164, 165) yield, on eliminating { e;}and [7;], the displacement form of the 
equilibrium equations 


> > FA KyKanl=F  (=1,2..n). (1.166) 
i=l kw 


32 Pin-jointed frameworks 


Since all members of the original structural layout of Section 1.1 are present, it 
follows that equations (1.166) always yield a unique solution for {uv}. 

The problem of determining structures of maximum stiffness will be inter- 
preted here as the problem of finding that structure of given volume of material 
V, which, under the action of the given forces, stores the least amount of strain 
energy U or in other words undergoes the least generalized displacement. 
Mathematically this can be stated as 


n 
minU= > Fy/2, (1.167) 
j=1 


subject to 


> Ah=V, (1.168) 


and to relations (1.163, 166). The Lagrangian U * for this problem can be 
written by (1.162) as 


n m m 
u*t=> Fuml2+2( AsV) + > ku; (Ao + p;-A) + 
k=1 i=1 i=1 


+¥ @/2) (4- 5 ie FAK Kn2ll), (1.169) 
j=1 i=1 k=1 


where A, {u;}, {vj} are multipliers, and conditions for U *to be a minimum, 
deduced by applying (1.160, 161) or directly, can be formulated as 


. > EA Ka Kyvjlls=Fr (k=1,2..n), (1.170) 


i=1 j=1 
A= > > EK jK in ¥j 0/21? + uj (i=1,2..m), (1.171) 
j=1 k=1 
and 
uyp,=0, p20 (@=1,2.. m). (1.172) 


Equations (1.170) for {y; }are identical with (1.166) for{ a} and since this 
last has a unique solution, it follows that (1.170) can be solved and that the 
solution is 


y= (=1,2-.n). (1.173) 


Equations (1.171) can now be solved for {ui} in terms of \ and an equation for 
d can be found by writing yu; = 0 for any i for which p; > 0 or, by (1.163), 
A, > Ag. Since Ao > 0, there must be at least one such i, because otherwise 
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(1.168) would be contradicted. It thus follows that a set of equations for the 

present problem like (1.155) is soluble or that the regularity criterion of 

(1.154) is satisfied. The derivation of (1.170—172) is thus justified. 
Substituting from (1.173) in (1.171) and then making use of (1.164), gives 


N= Ee}/2 + py; (@@=1,2..m). (1.174) 
Equation (1.172) shows that \ > 0 and so writing 
A=07/2E (1.175) 
and interpreting (1.172), using (1.163), it follows that 
€,=+0/E (A;>Ao) 
le] <O/E eet (@i=1,2..m), (1.176) 
or taking the limit Ao > 0, 
€;=+0/E (A;>0) 
le <O/E aes (@@=1,2..m). (1.177) 
Equation (1.165) gives 
T/A;=+o0 (=1,2..m;A;>0). (1.178) 


The structure of maximum stiffness thus satisfies the equilibrium equations 
(1.1), the conditions (1.177) for an allowable stress 0 and also allows its own 
displacement [v;] , which, by (1.177), gives a strain @/E in its tension members, a 
strain —0/E in its compression members and strains between these limits in 
members of the original layout, which are absent from this optimum structure. 
It follows by (1.70) that the present optimum structure of maximum stiffness 
is identical with the structure of least volume designed to carry the same loads 
[F)] with an allowable stress 0.4 The previous methods of resolving this last 
problem are thus available to solve the problem of maximum stiffness. 

It only remains to determine J. This follows from (1.168). If the Vinin for 
the equivalent structure of least volume is 


Vin = BFC, (1.179) 


where F is a typical force, / a typical length and § a number, like the 7 of (1.41), 
determined by the optimum analysis, then 


a= BFI/V, (1.180) 
where V is the given volume of (1.168).+ 


+t This result was first given by H. L. Cox, see Cox (1965). 

} Safety requires 3 < 0, which can be satisfied by reducing /, if necessary, since this does 
not affect stiffness. If F is actually required to be carried by the structure, V must be 
not less than sFl/o. 
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The generalized displacement 5 corresponding to F is given by 


F3=> Ry, (1.181) 


and so by (1.167), since the strain energy density is 0°/2E, 


5 =0° V/EF. (1.182) 
The ‘stiffness’ F/5 then follows from (1.180, 182) as 
F/5 = EV/6? 17, (1.183) 


and shows that the best material is that with the largest F. If the weight of the 
structure pV is given, then the best material is that with the largest E/p or the 
largest ‘specific modulus’. 

The problem of optimum elastic design in which min V is sought for given 
forces [F;] and a given displacement 5, defined by a formula like (1.181), is 
easily shown to lead to the same problem of optimum design, with allowable 
stress J, obtained above. The stress 3 is now given by (1.182), where V is given 
by (1.179) and so 


a= 5/61 (1.184) 


The prescribed 5 must not be too great or otherwise @ will exceed o and the 
structure will be unsafe. The general problem of optimum elastic design in 
which limitations are placed both upon stress and displacement presents difficul- 
ties because of the non-linearity of (1.165). Equations (1.164, 5) as well as one 
or more displacement restrictions have to be added to the equations of 

section 1.1 and this generates a problem of non-linear programming. Reference 
may be made to Pope and Schmit (1971) for details of numerical methods for 
tackling this kind of problem. 


[2 | Beams 


2.1 Beams of least volume of material 


Consider the problem of designing a beam of variable section to carry given 
normal loads over a span from x = 0 to x = /, with given boundary conditions 
at the ends. Let the distributed loads be w(x) per unit length and the concen- 
trated loads be w; at x = x;(i= 1,2 . .m), together with loads at the free ends 
(if any). Denote the shear force by S(x) and the bending moment by M(x). 
The conditions of equilibrium require that 


dS/dx +w=0, dM/dx +S =0, ea 
(AS)e=2;¢wi=0 (f= 1,2..n), 
and 
M=0(pinned ends); M=0, S given (free ends). (2.2) 
Equations (2.1) can be integrated in the form 
S = Mpo/l—M,/1-—dm/dx, 3) 


M =M,(1 —x/1) +M,x/1 +m, 


where Mo, M, are constants of integration (end moments) and m(x) is the 
bending moment for pinned ends. Fig. 2.1 shows a graph of a typical m(x) and 
a straight line with ordinates —Mo(1 — x/l) — M,x/I. The difference in the 
ordinates gives M(x). 

Let the beam be designed according to the methods of plastic design. Denote 
the required limiting plastic bending moment by M, (x) and let it have a lower 
positive bound M,o. It then follows that 


M(x) + p(x) =M,(x), p(x) 20, 
M(x) + q(x) = My (x), q(x)>0, (24) 
Myo + r(x) =M,(x), r(x) >0, 


where p, q, and r are non-negative ‘slack functions’. 

The relation between the area of cross-section of the beam and M, depends 
upon assumptions about the nature of the section. For a rectangular section of 
given width the area varies as M,”. For similar sections of given shape the area 
varies as Mp’ - For a beam of ‘sandwich construction’, with given depth and 
width but variable skin thickness, the area varies as M,,. This is also true for an 


idealized I-section of given depth and flange width but variable flange thickness. 
it will be assumed here that the area of the cross-section varies as Mp, where « 
is a positive number not greater than one. The condition for minimum volume 
of material can then be written 


t 
min v-{ KM ‘dx, (2.5)t 
0 


where k is a known constant depending upon the geometry of the section and 
the yield stress of the material. The ‘variables’ for this minimization problem 
are the functions S and M or the constants My and M,, together with the 
functions p, q, r, and M,. The ‘constraints’ are given by (2.1, 2, 4) or by 
(2.3, 4). This is a problem of the calculus of variations,t albeit a special one 
that can be readily reduced to an ordinary minimum problem. 

A solution to (2.5) will be sought on the assumption that M,, is a continuous 
function. Since the bending moment M(x) is continuous it follows by (2.4) that 
p, q, and r are also continuous functions. 


+ The case of ‘cost functions’ like (2.5) with a > | is treated in Marcal and Prager (1964) 
by an ingenious but special method. The present treatment may be applied to this case 
as well. 

t A good modern text is Pars (1962). The standard works do not however give much : 
attention to problems with non-negative functions and so this subject will be dealt with 
in some detail here. 


a 


Beams 377 


If p, q, r > O at a point and hence in an interval, it follows from (2.4) that 
M, can be reduced in magnitude in that interval. Equation (2.5) then shows 
that V can be reduced. It follows that one at least of p, q, and r must be zero 
at every point. If either p or q is zero, then (2.4) gives |M| = M, > Myo and if 
r=0 the conclusion is |M] <M, = Myo. The alternatives for M, are thus |M| 
or Myo and the greater of the two must of course be taken. The positions at 
which M = Mp are shown in Fig. 2.1 and an interval (a) in which M, =M, 
two intervals (b) in which M, =—M and two intervals (c) in which Mz, =Myo 
are also indicated. 

Equation (2.5) can now be expressed using (2.3) in the form 


min vip x (oc  x/D) +Myx[t +m dx + 
(a) 


(2.6) 
+ [ {aoc = sit) —Myxit~m)* dx + { KMS, dx, 
(b) e) 
where regions (a), (b), and (c) are defined by 
(a)p=0,q>0,r>0; (b)p>0,q=0,r>0; (c)p>0,q>0,r=0. 
(2.7) 


It is to be remarked that (2.4) forbids p = q = 0, but allows p =r = 0 and 

q =r =O in regions where m(x) is a linear function of x. These degenerate cases 
must be excluded, if (x) has such a linear part. This can be accomplished by 
adding a ‘perturbation’ to the loading, which consists of a uniform infinitesimal 
load 7 spread over the whole span. The conclusions of the present investigation 
can then be applied and when this has been done, the perturbation can be 
removed by allowing 7 > 0. It is clear that in practice the perturbation has 
only to be carried out in imagination. 

The problem of (2.6) is a minimum problem with variables My and M,. The 
variation of the boundaries of (a), (b), and (c) with 5M and 5M, produces no 
first order change in V, since M, is continuous at these boundaries and so the 
gains and losses of the various integrals in (2.6) cancel out. The conditions for 
the minimum of V in (2.6) can thus be written, remembering (2.2), as 


i) ME —xidae= [ My” *(1—x/l)dx or Mo given, 
Ka) ) (2.8) 
[ Mp Mxitax= [Mg x/dae ME a 
(a) » 


where 
M, =Min(a), M,=—Min(b) and M, = Mpo in (c). (2.9) 


Equations (2.8, 9) give two simultaneous equations for My and M, and complete 
the solution of the present problem. A solution can be obtained by trial and 
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error variation of My and M,, i.e. the position of the straight line in Fig. 2.1. 
For a cantilever or a pin-ended beam Mp and M, are given by (2.2) and (2.8) 
is not required. For a clamped-pinned beam M, is zero and only the first of 
(2.8) is needed. Both relations of (2.8) are required for a doubly-clamped 
beam. 

The meaning of (2.8) becomes particularly simple when a = 1. In this case 


f (l—x) dx = (l-x)dx or Mp given, 
fa) (b) 


| xar= [ xdx or M, given, 
(a) » 


and the equations state that the moments of the regions (a) and (b), taken about 
the ends x =/ and x = 0 respectively, are equal to one another. For the doubly- 
clamped case the lengths of (a) and (b) are equal and their centroids coincident. 


(2.10) 


2.2 An example 

Consider the problem of a doubly-clamped beam loaded by a uniform normal 
force —w per unit length acting on the intervals 0 <x <a and!—a<x <I. 

This gives a constant mm(x) for a <x </—qa and thus requires the introduction 
of a perturbation consisting of a uniform load —7 per unit length over the whole 
span. The unperturbed bending moment m(x) is given by 


m(x) = wax — wx?/2 (0 <x <a), 


m(x) = wa?/2 (a <x <//2). cM 


and is symmetric about x = 1/2. 

Assume for simplicity that M,o9 > 0 and that a= 1. This means for the 
perturbed loading that there are no regions (c) and that (2.10) gives the points 
where (a) and (b) meet as x = 1/4, 31/4. These values must satisfy M(x) = 0 and 
so by (2.3), since Mp = M, by symmetry, 

Mg + m(I/4) = 0. (2.12) 
Equations (2.11, 12) yield, taking n > 0, 
Mg = —wl(8a —1)/32 (a>1/4), 
Mo =—wa?/2 (a <1/4). 

The total bending moment M(x) = Mo + m(x) and M, = |M| for the present 

case. Equations (2.11, 13) thus give: 


(2.13) 


Case a > 1/4 
M, = wl(8a —1)/32 — wax + wx?/2 (0<x <//4), 
M, = ~wl(8a —1)/32 + wax — wx?/2 (/4 <x Sa), (2.14) 


M, =w(4a-1)/32 (aSx <1/2), 


Case a <1/4 
M, =w(a-x)?/2 (0<x <a), 


M,=0 (a<x<i/2). (2.15) 


Equations (2.14, 15) define the optimum design. In the case a <//4 the beam 
divides into two separate cantilevers. 
The volume Vyin is given by 


Venin = Wk(1? — 12la? + 481a? — 32a°)/96 (a > 1/4), 


Vinin = wka?/3 (a <l/4). (2.16) 


2.3 Application of the calculus of variations 


The conclusions of Section 2.1 may be obtained by a direct use of the calculus 
of variations. This alternative approach, although much less simple than that of 
Section 2.1, does throw new light on the problem. It also provides the oppor- 
tunity to rehearse the techniques of this calculus, making use of a simple problem 
whose solution is already known and which is capable of direct graphical 
representation (Fig. 2.1). 

It is necessary first of all, as it was in Section 2.1, to remove the degenerate 
cases p =r =O and g =r = 0 by adding a uniform infinitesimal loading n. It is 
also necessary, for the special case when Mo = 0, to remove the even more 
degenerate case p = q =r =0. This can be done by introducing a positive Myo 
and letting Mo 0 after the analysis is complete (cf. Section 2.2). 

The minimum solution must consist, after (2.7), of regions (a), (b) and (c). 
lt will be assumed that the span of the beam can be divided into a finite number 
of intervals for which one or other of these regimes is valid. 

Consider now ‘variations’ 5S(x), 5M(x), 5p(x), 5q(x), 5r(x) and 6M, (x) 
from an optimum solution S(x), M(x), p(x), q(x), (x) and M,(x), which is 
assumed to exist. These must satisfy, by (2.1, 2, 4), 


d5S/dx=0, d5M/dx+5S=0, 5S continuous, (2.17) 

5M = 0 (pinned or free ends), 5S = 0 (free ends), (2.18) 

56M +5p=5M,, —5M+65q=65M,, 5r=5M,, (2.19) 
and 

6p2>0(p=0), 5g>0(q=0), 5r>0(r=0). (2.20) 


They must also imply, by (2.5), that 
t 
ov=[ akMs * 5M, dx >0. (2.21) 
0 


The next step is to solve (2.17, 19) and thus to divide the variations into 
‘dependent’ and ‘independent’ variations. The dependent variations must be such 


that no additional restriction is placed upon their values and they are thus able 
to take any values the independent variations impose upon them. The independent 
variations must include all those which are subject to restrictions like (2.18, 20). 
A problem in which this division of the variations can be carried out is a 
‘regular’ problem, otherwise it is ‘irregular’. The reason for the exclusion of 
degenerate cases, in which two or three of p, q, and r are zero, is that (2.20) 
requires that two or three of the corresponding variations must be taken as 
independent. It is easily seen that this is not possible. 

Equations (2.17) can be solved to give 


5S =5Mo/l—5M,/l, 5M =5Mo(1 —x/l) + 5M,x/l, (2.22) 


where 5M and 5M, are variations of the constants Mo and M, of (2.3). These 
last may be restricted by (2.18) and should be taken as independent. Equation 
(2.22) shows that this is possible since 5S and 5M have no restrictions placed 
upon them, in the open interval (0, /) at any rate. Equations (2.19) may be 
solved for 5M, and for two of dp, 5q, and 5r in terms of the remaining 
variation of a slack function and, after (2.22), in terms of 5Mo and 5M,. This 
is satisfactory for each of the regions (a), (b), and (c), since by (2.7) only one 
of p, gq, and r vanishes and so by (2.20) only one of the corresponding 
variations need be independent. The scheme for the solution of (2.19, 22) is 
thus 


Region Dependent variations Independent variations 
(a) 5S, 5M, 5M,, 5q, 5r 5Mo, 5M, 5p 
(b) 5S, 5M, 5M,, Sp, dr 5Mo, 5M, 54¢ 
(c) 5S, 5M, 5My, dp, 54 5M, 5M,, 5r 


(2.23) 


It is unnecessary to write the actual solutions. 

Lagrangian ‘multipliers’ must now be introduced just as in (1.152). It is 
convenient to introduce the constant k as a factor, to introduce a positive 
infinitesimal € as a divisor and to write the multipliers corresponding to the 
first two of (2.17) and to (2.19) as ku(x)/e, KO (x)/e, kx '(x)/e, kx" (x)/e and 
kko(x)/e. Equations (2.17, 19, 21) then give 


i 
{ {eams* 5M, + vd5S/dx + 0(d5M/dx + 5S) + x'(5M + 5p —5M,) + 


P 2.24 
+ x"(COM + 5q—BM,) + xo(6r~8Mp) | dx >0. a9 
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Integrating by parts, assuming v and @ continuous, and collecting like terms 
gives 
[ {(cams* —K'—K"—Kko)5M, + (0 — du/dx) 5S + (k’ —K" — d0/dx)5M + 
: (2.25) 
+x'Spt+x"5q+ Kor} dx + [v6S + 96M], > 0, 


which is valid for all variations satistying (2.17, 19, 21). 
The Lagrangian multipliers must now be chosen so as to remove from (2.25) 
the dependent variations of (2.23). This requires 


k' +k" +Kko=eoM, +, (2.26) 
dujdx=0, d0/dx=x'—K"=k, (2.27) 

and 
k'=0(p>0), x"=0(q@>0), ko=0(r>0), (2.28) 


where a new function k has been introduced. Equations (2.26, 28) can always be 
solved for x’, k” and Ko in each region (a), (b), and (c). Equations (2.27) can then 
be integrated over the whole span to yield continuous functions v and 4. The 
required choice for the Langrangian multipliers is thus possible. 

Imposing (2.26—28) on (2.25) then gives, using (2.22) for the integrated 
terms, 


[ eoMS*5pdx+ | caMy *5qdx+ eaMe* br dx + 
» ‘() te) (2.29) 
+ [{) — v0) }/1 — 0(0)] Mo — [{v@) — vO) }/!— A) 5M; > 0, 


where all the variations are independent and can be given arbitrary values, 
subject to (2.18, 20). This means that each term in (2.29) must be non-negative, 
since the other terms can be set equal to zero. The first, second, and third terms 
of (2.29) are non-negative, thanks to (2.20). The last two terms can be zero if 
Mo or M, are given by (2.2), but otherwise they can always be made negative, 
unless the coefficients of 5M and 5M, are zero. This gives 


{v(I) - v(0)}/1=0(0) or Mo ag (2.30) 
{o(1)—vO)}/1= 01) or M, given. 
Equation (2.27) gives 
{ K(1 ~x/l) dx = -0(0) + v(0) — o(O)}/, 
(231) 


i 
[ k(a/l) dx = 0(1) ~{ 01) -v0)} 
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and (2.26, 28) yield 
k=eaMg~*in(a), x=-eaM¢*in(b), x =Oin(c). (2.32) 


Substitution from (2.32) into (2.31) and from thence into (2.30) then gives 
(2.8) once more and the present investigation agrees with that of Section 2.1. 

The Lagrangian multipliers can be given a kinematic interpretation. If vis a 
virtual normal displacement, then by (2.27) @ is the corresponding virtual 
rotation and x the corresponding virtual curvature of the beam. Equations (2.30) 
then give, on removing a rigid body motion by writing v(0) = v(/) = 0, 

6(0)=0 or Mo=0, O(1)=0 or M,=0, (2.33) 

which are the usual kinematic boundary conditions for the doubly-clamped and 
clamped-pinned cases. 

In the special case when M9 > 0, the solutions obtained above can be 
treated as optimum elastic designs. Regions (c) are now absent and in regions (a) 
and (b) the curvature « is given by (2.32), while the corresponding bending 
moments are M = + M,. It follows that 


M/k =M2 “Jea. (2.34) 

For sections, which have second moments of area / given by 

1=kyM2-*, (2.35) 
where k, is a constant, it is only necessary to write 

e€=1/Eak,, (2.36) 
where E is Young’s modulus, to convert (2.34) into the usual bending moment- 
curvature relation for the elastic deformation of beams and hence to allow u(x) 
to be identified with the real displacement. The relation (2.35) is valid for 
rectangular beams of constant width (a = 1/2), similar sections of varying size 
(a = 2/3) and for sandwich beams and idealized I-sections (a = 1). For the elastic 


case M,, should be taken as the bending moment at which plastic flow begins, 
with k of (2.5) defined to suit. 


2.4 Use of the Lagrangian 


The results of Section 2.3 can be obtained, formally at any rate, by constructing 
a Lagrangian in a manner analogous to that used to form (1.62), and then writing 
the conditions that this Lagrangian should have a minimum. 

The Lagrangian V* for the problem of (2.1, 2,4, 5) can be written 


evra fea + v(as/ax + w) + (am/ax +S) +x'(M +p-M,)+ 
F (2.37) 
+"(-M + q ~My) + Ko(Mpo +r-M,)\dx+ 3 v, [AS enn +i}, 
jel 
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where v(x), 0(x), k(x), k" (x), Ko(x) and v;(i= 1, 2 . . 7) are Lagrangian 
multipliers and € the positive infinitesimal introduced before. Conditions are 
to be sought for min ue subject to the boundary conditions (2.2)+ and the 
non-negativity conditions for p, g, and r. 

Integration of (2.37) by parts gives 


1 
ev* -*[ {eg + x" + 0)My + (6 —du/dx)S + (x — x" — d0/dx)M + 


+x'p+x"qtKkortowt KoMyo} dx+k > [fe - v(x) (AS) =x; + vw] + 
i=1 


+k[vS + 6M]°, (2.38) 
and so the condition for V* to be a minimum can be written 


[ ((coms* —K'-K"—Ko)5Mp + (8 — dv/dx)5S + (k'— x” — d0/dx)5M + 


n 
+k'5p+x"5qt+kodr}dx+ > {v;—v(%;) } (A5S) =<; + [v5S + 95M], >0. 
a4 (2.39) 
The variations in (2.39) can be taken as arbitrary apart from dp, 5q, and 5r, 
which must satisfy (2.20) and 5S, 5M at the ends, which by (2.2) will be zero, 
when not arbitrary. This means that each term in (2.39) involving a variation 
as a multiplier in its integrand must be non-negative and indeed the integrands 
of such terms must themselves be non-negative for all x. It follows that the 
factors multiplying completely arbitrary variations must vanish and those which 
multiply non-negative variations, like 5p for p = 0, must themselves be non- 
negative. Application of these rules to (2.39) gives 


k' +k" +ko=eaMy ', 

dv/dx =0, d@/dx =x'—K", 

«'=O0(p>0), x’ >0(p=0), 

k"=0(q>0), x" 20(q=0), (2.40) 
Ko =O(r>0), x. 20(r=0), 

vu, = v(x;) (@=1,2..n), 

v=0 or Sgiven; @=0 or M given (at ends x =0,/). 


1 These can be included when well defined. If M, = 0, a term —0,M, can be added to (2.37), 
where 0, is a constant multiplier. 


For the regions defined by (2.7) the equations (2.4U) give the same results as 
those of Section 2.3. The first three of (2.40) agree with (2.26, 27) and the 
conditions for p, g, and r > 0 in (2.40) and (2.28) are the same. It follows that 
the values of x’, x”, and ko agree in the regions (a), (b), and (c) and that (2.32) 
follows from (2.40). The boundary conditions (2.33) also follow from the last 


line of (2.40). Equations (2.8) can thus be derived just as they were in Section 2.3. 


The necessary conditions for the minimum of (2.5), subject to (2.1, 2, 4), can 
thus be obtained by writing the necessary conditions for the minimum of the 
Lagrangian (2.37), subject to (2.2) and p, g, and r > 0. This conclusion has been 
established on the assumption that the minimum solution is non-degenerate. 
However the results of (2.40) are valid for degenerate cases as well, when 
properly interpreted, just as were the results of Sections 2.1, 3. 


2.5 Dual form of the linear problem 


The Lagrangian of (2.37) becomes a linear functional when @ = | and for this 
special case a theory of duality can be developed, very similar to that for 
linear programming in Section 1.4. 

The dual problem corresponding to the primal problem of (2.1, 2, 4, 5) 
can be obtained from the Lagrangian (2.37) by imposing the conditions (2.40) 
as constraints.} The resulting problem will be shown to be a maximum problem 
and can be obtained by substituting from (2.40) into (2.38) and by eliminating 
p,q, 7, S, and M from (2.40). This gives the dual problem in the form: 


max W = (k/e) ({ (wu + MyoKo) dx + 3 wivtx)} ‘ (2.41) 
i=1 


subject to, 
k' +k" + Ko =e, 
x’, K", Ko>0, 
dv/dx=6, dO/dx =x'-K", 
v=Oatasupport, @=0ata clamp. 


(2.42) 


Let x’, x”, Ko, v and @ be a ‘feasible solution’ of the dual, i.e. any solution 


t See the discussion of (1.62) above. 
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of (2.42). Let S, M, p, q, r, and M, be a feasible solution of the primal, i.e. any 
solution of (2.1, 2, 4). For these 


W =(k/e) (f (-vdS/dx + M,oko) dx + = wivta)} , by (2.41) and (2.1), 


i=1 ’ 


=(k/e) [ (8S + Myoko) dx, by (2.1), (2.2), and (2.42), 

=(k/e) [ (-8dM/dx + M,oko) dx, by (2.1), 

=(k/e) [ {moe -K")+Myoko Jax, by (2.2) and (2.42), 
wef or, — p)x' + (Mp —q)k" + (My —r)ko las, by (2.4), 


1 
= V-(k/e) { (pk' + qk" + rK9) dx, by (2.42) and (2.5), 


< V, by (2.4) and (2.42). (2.43) 
It follows that, 
W< Vain <V, (2.44) 


which shows that W is bounded above and that upper and lower bounds can be 
found for Vain by using feasible solutions to the primal and dual respectively. 

If it is now assumed, as in Section 2.3, that the primal problem has an 
optimum solution, then feasible solutions to the primal and dual problems can 
be found, which satisfy (2.40), the necessary conditions for a minimum of the 
primal. These give by (2.40, 43) W = V and so by (2.44) W = Vain. This W is 
also Wx, a8 follows from (2.44), and so the dual has an optimum solution 
with 

| el oe (2.45) 


The volume of material in the optimum beam can thus be found by solving the 
dual problem. The present discussion can be compared with that of Section 1.4. 

Use of the kinematic interpretation of the multipliers allows, when M,o > 0 
at any rate, the interpretation of W as the virtual work of the given loads taken 
over the virtual displacements v. Equation (2.42) imposes the restriction that 
the absolute magnitude of the virtual curvature |x’—«”| may not exceed e. 

Necessary and sufficient conditions for an optimum design can be formulated 
for beams, just as in (1.70) for frameworks. The present duality theory gives the 
following result: 


A distribution of M, gives an optimum design, if M, =| M|for|M|> Myo and 
My ™ Myo for |M| <Myo, where M is a possible statically correct bending moment 


distribution, and if there is a virtual displacement v, which satisfies the kinematic 

boundary conditions and has a virtual curvature € when M = M,, a virtual 

curvature —e when M = —M, and a virtual curvature of zero when| M|<Myo- 
(2.46) 


This is a Michell theorem for beams like that of (1.70) for frameworks. 


2.6 Beams of maximum stiffness 

The present problem, like that of Section 1.12, will be interpreted as meaning 
the determination of that structure of given volume, which stores, for given 
loads, the least amount of strain energy U. For a beam, loaded as in 

Section 2.1, this means 


min 2U -[ wodx + ¥ w;v(x;)+ [se]; (2.47) 
i=1 


where v(x) is now the real displacement of the beam and the last term has been 
included to allow for the case of loaded free end. 

The second moment of area / is assumed to be related to the area of the 
cross-section A by 


A=kol® (2.48) 


where Ko is a given constant and § a positive number not greater than one. For a 
rectangular section of given width 6 = 1/3, for similar sections 6 = 1/2 and for a 
sandwich beam or an idealized I-section 8 = 1. The condition of given volume V 
can thus be written 


[ P ax = V/ko. (2.49) 
0 


The equations governing the deformation of the beam are 
dS/dx + w =0, dM/dx +S=0, (AS),-.;+w;=0 ((=1,2..n), 
M=Elkx, x=d0/dx, 6 =dv/dx, (2.50) 
S given, M=0 (free end); v=M=0(pinned end); v=6=0 
(clamped end). 
A positive lower limit J is placed on J by the relation 
iotr=I, r20. (2.51) 


The problem of maximum stiffness for a beam is thus to find the minimum U 
of (2.47), subject to (2.49, 50, 51). The design is determined by /(x) and the 
minimum problem involves in addition the functions S, M, x, 0, v, andr. The 
methods used in Section 2.3 may be applied to the equations as they stand, 
but it is easier to reduce the size of the problem first. 


Beams 41 
Use of (2.50) or of a standard formula enables (2.47) to be replaced by 
1 
min 2u=[ (M?/ET) dx, (2:52) 
where, again by (2.50) or by (2.1, 2), M is expressed in terms of Mg and M, by 
(2.3). The values of Mo and M, for a given /(x) follow from (2.50), using (2.31), 
or more usefully, using Castigliano’s theorem, from 
3U/dMy=0 or Mo given, dU/aM,=0 or M, given. (2.53) 
Let [(x) be varied to (x) + 5/(x). This will induce changes in My and M,, 
unless one or both of them are given. Such changes themselves will not affect U 
in (2.52), thanks to (2.53). The condition that (2.52) gives a minimum can thus 
be expressed by 
[ (-M?/EI*) &Idx > 0, (2.54) 
or, using (2.50), by 


1 
{ x?5/dx <0. (2.55) 


Equation (2.55) will be valid, by (2.49) and (2.51), for variations which satisfy 
3 "sax =0, 51=6r, 6r>0 (r=0). (2.56) 
0 


The span of the beam can be divided into two regions defined by (a) r>0 
and (b) r = 0. Equations (2.55, 56), with 5/ eliminated, can then be written 


[ orde+ | x?5rdx <0, 

(a) (b) 

when (2.57) 
[ Psrax + | 8 5rdx=0 and 5r>0in(b). 
(a) (b) 


Taking 5r = 0 in (b) and since r is not restricted as to sign in (a), it follows that 


i) x?8rdx=0, when | Pr dx =0. (2.58) 
« 4a) 


48 Beams 
This impliest 

x? =CP-) in(a), (2.59) 
where C is a positive constant. Substituting from (2.59) in (2.57) gives 


i) (cig —x«?)6rdx>0 (2.60) 
(b) 


and so, since by (2.57) 5r > 0 in (b), 
x? <ci8-}, in(b). (2.61) 
Equations (2.59, 61) combine to give, recalling (2.51), 
eScri Th); 
<a G=h). 


The solution to the present problem is thus given by (2.50, 62), with the 
constant C determined by (2.49). 
In the case when J, > 0, equations (2.48, 62) give 


Elk? |2A = EC/2ko, (2.63) 


which says that the average density of strain energy at a section of the beam is 
the same for all sections. This rule determines the design of a beam of maximum 
stiffness. Exactly the same rule can be established for the elastic design of beams 
of minimum volume of material considered in Section 2.3. The area of section 
was written there as KM, and / was given by (2.35). The bending moment M was 
always +M, and so 


(2.62) 


M?/2EIA = 1/2Ekyk, (2.64) 


which again gives constant average strain energy at a section. The two kinds of 
optimum design are thus identical in form. The volume can be made identical 
(equal to V of (2.49)) by an adjustment of k in (2.5). The overall displacement 
follows from the total strain energy U, which by (2.64) is given by 


U= V/2Ek,k, (2.65) 


where an adjustment must be made to k, , corresponding to that made to k.§ 


+ This is a standard lemma of the calculus of variations. Equations (2.58) imply 
f (x? [19-1 — C)P—15r dx = 0, with P~15r = «2/171 — C and the constant C deter- 
(a) 


mined by [  (x?//2-! — C) dx = 0. It follows that r («3/181 _ C) dx = 0 and hence 
fa) fa) 
(2.59). 
+ & can be varied by varying the assumed yield stress. 
§ If k is multiplied by p, then k, k must be muitipued by p”/“, 


2.7 Pin-ended strut of least volume of material 


Consider the problem of designing a pin-ended strut, with the least volume of 
material, to carry an end load P over a span /. Let the relation between the least 
second moment of area / of the section and the section area A be given by 
(2.48). The objective can thus be expressed by 


min v-[ kol® dx. (2.66) 


The design of the strut is defined by /(x) and, corresponding to any such 
design, there will be a least critical load or Euler buckling load Pg with its mode 
of buckling defined by a displacement v(x). The following equations will be 
satisfied: 


EI d*v/dx? + Pgv = 0, (2.67) 
v(0) = v(/) = 0, (2.68) 
where the solution with the smallest Pg is to be taken. The function v(x) can be 
fixed completely, apart from sign, by imposing the condition 
i] 
{ (dv/dx)? dx = 1. (2.69) 
Equations (2.67—69) give the Rayleigh formula 
Pe =k EI(d?v/ dx? dx. (2.70) 


The safety of the strut is ensured by 
P+p=Pg, p2>o. (2.71) 


However if p > 0, a design, with J replaced by PI/Pg, which by (2.67—69) will 
have the same mode of buckling v(x) and a critical load P, will give (2.66) a 
smaller value of V. It is therefore possible, for the present purpose, to replace 
(2.71) by 

P=Px. (2.72) 


A variation 5/ from the minimum solution must satisfy, by (2.66, 72) 
2 1>"51dx>0, when 5Pz =0. (2.73) 
A change of / to / + 5/ will produce a corresponding change in buckling mode. 


Such a change will however only induce a second order change in 5P_ (Rayleigh’s 
principle), since Pg is the minimum of the integral of (2.70) subject to (2.68, 69). 


It results that the condition of (2.73) may be written 
1 
5Pp = | E$I(d? v/dx?)? dx =0. (2.74) 
0 
Assuming that for a minimum? 
I(x)>0 (O0<x<i), (2.75) 
it follows, since 5/ may have its sign reversed without becoming inadmisible, 
that 
ree 7 
[ -*§rdx=0, when i E51(d?v/dx?)? dx = 0. (2.76) 
0 0 


This is of the same form as (2.58) and so yields, as in (2.59), what is in fact the 
same relationship, 


(dv/dx?)?=CP®-* — O<x<I). (2.77) 


The solution to the present problem is thus governed by (2.67—69, 72, 77). 
Equations (2.67, 77) give 


T=Cyv7t6*2). (2.78) 


where C, is a positive constant and v is taken with a positive sign.+ Equations 
(2.67, 72, 78) give 


d?u/dx? + (P/EC,)v®~ P/F *Y =, (2.79) 


The integral of (2.79), which satisfies (2.68), can be written in the parametric 
form 


x =1F(0, B)/F(x, 8), v= {I/F(m, B)} 8 * *{PB/EC,(B + 1)} 8 * ?/? (sin 8) E* VB, 
(2.80) 
where 


6 
F(@,8)= | Gin@yt 40 (2.81) 
0 
The constant C, may be chosen to satisfy (2.69), but there is no need to do 
this, since it does not affect the design. Equations (2.78, 80) give 
1={I/F(x, 8)}? {PB/E(B + 1)} (sin@)?, (2.82) 
and hence by (2.66) 
Vexin = kof l/F(n, 6); ** (PB/EG+ 1} Ffm, B2B+1)}. (2.83) 


7 A zero of /(x) not at an end will introduce a pinned joint at the same place. It also gives, by 
(2.67), v = 0 which cannot give minimum Pr. 
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The following results derived from (2.81) are required for applications: 


F(@,1)=1-—cos®@, F(z, 1)=2, 
F(@, 1/2) =/2—(1/4) sin 20, F(a, 1/2)= 7/2, 


F(@, 1/3) = 2/3 — cos @ + (1/3) cos*@, F(m, 1/3) =4/3, C8 
F(x, 1/4) = 32/8, F(x, 1/5) = 16/15. 

Equations (2.83, 84) give for the minimum volumes: 

V =(1/12)ko(P/E)I? (8 = 1, sandwich beam of given overall dimensions), 


V =(V3/2m)ko(P/E)'!? 1? (8 = 1/2, similar sections of varying size), 
V = (37/3/5)ko(P/E)'/3 15> (6 = 1/3, rectangular section of constant width), 
(2.85) 


which may be compared with the result for a uniform strut, given by 
V = ko (1/n#) (P/E) 178 **. (2.86) 


The volume ratios for 8 = 1, 1/2, 1/3 respectively are x? /12 ~ 0-82, 3/2 ~ 0-87 
and (37)7/3/5 ~ 0-89. These give a measure for the saving of material obtained 
by tapering a strut.7 


2.8 Structures built from beams 


Examples of structures built only from beams are those stiff-jointed plane 
frameworks, whose stability derives from their stiff joints, and grillages made 
up of beams with low torsional stiffness. If, for simplicity, each component 
beam is assumed to have a uniform section, the optimum design of such 
structures reduces itself to a problem of linear programming, when a cost 
function is used, which is linear in the unknown limiting plastic moments. The 
two kinds of structures present identical mathematical problems and the treat- 
ment given here is applicable to both. 

Consider a structure of given layout loaded by given concentrated forces 
F,(1=1,2..n). Each portion of beam between nodes has a uniform section 
and the ‘design’ of the structure is determined by the limiting plastic moments 
M,j(i = 1, 2 . . m) of these beam elements, which are m in number. The ‘cost’ 
of the beams is assumed to be jointly proportional to their M,;¢ and their 
lengths /;(i = 1, 2 .. m). The objective can thus be formulated as 


m 
min V = > Myil;- (2.87) 
i=l 
+ This is by no means a new problem. The result for 6 = 1/2 was given by Clausen at 
St. Petersburg in 1851. A recent paper which places lower bounds on / is Trahair and 
Booker (1970). 
+ Costs proportional to a power of M,,; are considered in Prager (1956). 
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The method of plastic or limit design will be used. This requires the 
consideration of all possible plastic hinges, those at the ends of each member 
and those at internal points where forces are applied. Since it is not known in 
advance which M,; at a node is the least, hinges must be given to each member 


at that node. The structure becomes a mechanism, when the hinges are introduced. 


Let the number of degrees of freedom be fand let ¢, (k= 1,2 . . f) be a set of 
Lagrangian coordinates, which define completely the ‘modes of collapse’ of the 
mechanism. Denote by 6;(j = 1, 2 . . A) the rotations of the hinges and by 
u,(1= 1,2 ..n) the displacements corresponding to the forces [Fj] . The geometry 
of the structure determines the relations 
f f 
6;= D Tee G=1,2-.h), w= Y Unde W=1,2-.n), (2.88) 
k=1 k=1 


where [7j;,] and [Uj] are known matrices. 


If M;(j = 1, 2 . .h) are the bending moments at the positions of the hinges, 
which correspond to {6;}, then the principle of virtual work gives 


h n 
> Mo;= > Fim. (2.89) 
j=1 


1=1 
Substitution from (2.88) into (2.89) yields 


h f n f 
EM LY Toe = TM YL Unde, (2.90) 
l=1 k=1 


fob k=1 


which is true for arbitrary {¢,}. The equations of equilibrium for the moments 
[M,] are thus 


h n 
YMTn= > AU (k= 1,2... (291) 
j=1 Ii=1 


The safety of the structure requires, in addition to satisfying (2.91), that the 
moments [Mj] are limited by 


Jmuil< 5 tiMps = G=1,2--H), (2.92) 


where 4; takes the value unity if the jth hinge lies on the ith member, but is 
otherwise zero. The limiting plastic moments are assumed to be bounded below 
or that 


Mp >My >0 (i= 1,2... m). (2.93) 


The problem of (2.87, 91—93) is a problem of linear programming.+ Introduc- 
sng slack variables p; and q;(j = 1, 2 . . h) in (2.92) and surplus variables 
r;(i=1, 2. . m) in (2.93) puts the constraints into the equation form 


m m 
M; + Pj)= L ujiMpi, Mj +4; = LY wyiMpi, where pj, q; >0(j=1,2..h), 
i=1 i=1 


M,;=Myo +17;, where r;>0(§=1,2.. m). 
(2.94) 


The conditions for minimum V of (2.87) and the dual problem are most 
readily obtained by constructing a Lagrangian as in (1.62). Introducing 
Lagrangian multipliers $,(k = 1,2 . .f), 0;,0;(/=1,2..h) and a;(i=1,2..m), 
as well as a positive infinitesimal €, yields, on combining (2.87, 91, 94), 


m ft n h 
VS [eM + (Moo +7: Myi)|+ & ou(S Fin MT) + 
i=1 k=1 t=1 i=1 / 


h m m 
+> (0; (™ +pj- 2 HM) +6; (™“ +a, HeMyi). (2.95) 
j=1 i=1 i=1 


The conditions for min V *, subject to the non-negativity conditions of (2.94), 
can now be written as 


h 
>2 (8; + 07) uy; +. = El; (i=1,2..m), 


7=2 


0;=0;-0;= > Ted (f=1,2..h), 
k=1 (2.96) 
6;=0(p;>0), 6;>0(p;=0) (f=1,2..h), 
6;=0(q;>0), 0;>0(q;=0) (j=1,2.-.h), 
a;=0(7;>0), a;,>0(7;=0)  (=1,2..m), 


which by the general theory of Section 1.4 give necessary conditions for min V 
of (2.87) and together with (2.91, 94) give sufficient conditions for this same 
minimum. 

The Lagrangian multipliers can be given a kinematic interpretation. If {on} 
are taken as generalized coordinates as in (2.88) defining a virtual deformation 
(collapse mechanism) of the structure, then {6} , as defined in (2.96), gives by 
(2.88) the corresponding hinge rotations. The quantities {¢,} and {0;} determined 
by (2.96), do in fact define the collapse mechanism for the failure of the 
optimum structure. This may be seen by considering the possible forms that the 


| Methods for the solution of this problem are given by Heyman (1951), Livesley (1956) 
and Heyman and Prager (1958). 
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optimum solution can take for the jth hinge on the ith member (4; = 1). If p; 
and q; are both zero, then (2.94) gives a contradiction. The possible cases are 
thus: 


(a) pj =0, q; >0, r;>0 with 0;=6;>0, 67=0, M;=M,;>Mpo, 

(b) p; >0, q;=0, r;>0 with 0;=-07 <0, 8;=0, Mj =—M,;<—Myo, 

(c) pj >0, q;>0, r,>0 with 0; = 0; = 0; =0,|M;|<M,:, Mpi >Mpo. 
(2.97) 


where (2.94, 96) have been used. The hinge rotations which occur in the virtual 
deformation are thus associated with bending moments of the same sign at their 
limiting plastic values. 

It is clear from (2.97) that | 6; | = 6; + 8; and so the first of (2.96) may be 
written as 


Pa |0;|+a,=eh, wherea;>O  (i=1,2..m), (2.98) 
Hinges on 
ith member 


or recalling (2.94) as 


> | 94| =el; (alli with M,; > Myo), 
Hinges on 4 
ith member 


(2.99) 
2 | 9 <el; (alli with M,; = Myo), 


Hinges on 
ith member 
which are conditions for the optimum structure, expressed in terms of its collapse 
mechanism, first given in Foulkes (1954). 
The dual form of the present problem can be obtained by imposing (2.96) on 
(2.95). This gives, introducing {u;} by (2.88) 


n m 
max W=(l/e)! > Fu; +Mpo > a}. (2.100) 
'=1 i=1 

which is, apart from the term in M,o, the virtual work of the given forces taken 
over the virtual deformation. The constraints for (2.100) are (2.88, 98). Any 
solution to the constraints gives a W in (2.100) which is a lower bound to Vain- 
The solution to the dual problem gives Wax = Vnin and also determines the 
active hinges, which by means of (2.91, 97) lead to the optimum design. 


2.9 Example of a portal frame 


Consider the problem of Fig. 2.2(a) with M,o = 0 for simplicity. The possible 
plastic hinges / = 1, 2. . 7 are shown and when they are present the resulting 


uy Pz Ur 
%; 3 Ps 
(b) (c) 
Fig. 2.2 


mechanism has four degrees of freedom. The deformations corresponding to 
each of four Lagrangian coordinates ¢,(k = 1, 2, 3, 4) are shown in (b) and (c). 
Equations (2.88) now take the form 

4, =d1, 92 =$1 — $3, 93 = 2 — $3, 94 = 262, O5 = G2 + Oa, 


(2.101) 
66=61-—¢4, 97=61, and uy =/9,, uz =162/2. 


The equations of equilibrium (2.91) can be written down, using the principle 
of virtual work and (2.101), for each of the deformations ¢(k = 1, 2, 3,4). 
This gives 


M, +M,+M.,+M,=FI, M3 + 2M, + Ms; = FI, -M,-M3=0, M; —-M,=0. 
(2.102) 


It is to be expected that M,,.M,, and M, will be positive and at their 
limiting values. Equations (2.102) will allow the same for M; and Mg, but not 
for M, and M,. Equilibrium will allow six bending moments to take the three 
limiting values for the members and will determine these limiting values and the 
seventh bending moment. Take as a trial 


M,=My,, M3=Mg=Ms=M,2, Me =M7z=Mgs. (2.103) 
Equations (2.102, 103) then give 
M,, = 3F1/4, Mp2 =M,3 = FI/4, M, =-Fl/4, (2.104) 


which determines a safe design since | M3] <M, ,. Equations (2.97) give 82 =0 
and 0, > 0(i= 1,3. .7). The first of (2.99) and (2.101) give 


, =el, 46. —$3 +d, =el, 26,-¢4=el and ¢$,=¢3, (2.105) 
with a solution, 


o,=3=¢4=6€l o2 =El/4. (2.106) 
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Equations (2.101, 106) give 8; <0, whereas (2.103, 104) give M3 > 0, contrary 
to (2.97). The trial of (2.103) has not led to an optimum solution. 

The first trial of (2.103) led to a contradiction at the third hinge. This 
suggests relaxing conditions there and writing 


M, =M,=M,;1, My=Ms=Mp2., Me =M7=M,3. (2.107) 
Equations (2.102, 107) give 
M,, = FI/8, Mp2 =Mp3 =3Fl/8, M3 =-FI/8, (2.108) 


with | M,| < Mp2. Equations (2.97) give @3 = 0 and the rest non-negative. 
Equations (2.99, 101) give 


26,-63=€l, 362+¢4=el, 26:-¢4=e1 and ¢2.=$3, (2.109) 
with a solution, 
$, = Sel/8, 2 = 63 = bq = El/4. (2.110) 


Equations (2.101, 110) now give 4; > 0 (all i, i # 3) and 43 = 0, which by (2.97) 
is consistent with (2.107, 108). The second trial gives an optimum solution. 
Further confirmation can be obtained by calculating Vin from (2.87, 108) and 
Wnax from (2.100, 101, 110). The result is 


Vina = Weiax = TEI2/8. (2.111) 


3 | Circular sandwich plates 


3.1 Plates of least volume of face material 


Consider the problem of designing a circular sandwich plate of radius R to 
balance a system of given normal forces or to transmit such forces to supports 
at its edge. The sandwich will be assumed to have uniform thickness h, with 
faces of variable thickness f(r), depending only upon the polar coordinate r 
(Fig. 3.1). The filling will be assumed to be able to carry such forces as may be 
required of it. 


L t(r) u(r) 


Fig. 3.1 


The given forces are assumed to be axially symmetrical. They are defined by 
the shear force per unit length (stress resultant) Q(r), which gives the total force 
normal to the plate, within a circle with centre O and radius r, as 27rQ. The 
function rQ(r) has discontinuities at circles of concentrated forces and tends to 
a finite limite as r > 0, if there is a concentrated force at the centre. 

Let the bending moments per unit length (stress couples) be M, and My in a 
polar coordinate system with origin at the centre of the plate. The conditions of 
equilibrium can then be written as 


d(rM,)/dr — Mg = rQ (3.1) 
and 
M,(r) continuous; M,(R) = 0 (edge free or simply supported). (3.2) 
The safety of the plate is assured if 


|M,|, |Mo|, |M,—Mg| <ohe, (3.3) 


58 Circular sandwich plates 


where a is an allowable stress and the Tresca criterion is assumed to apply to 
the material of the faces. The relations (3.3) may be written as equations, on 
introducing slack functions p,, q; (i= 1, 2, 3), in the form 
M,. + p, = oht, —M,. + q, = oht, 
Mp + p2= oht, —Mo + q2= oht, (3.4) 
M,— Mo + p3=oht, —M,+ Mg + q3= oht, 
where 
Pi. G1, Pa, G2, Ps, 4320. (3.5) 


Equations (3.3, 4) are given a geometrical representation in Fig. 3.2. The point 
(M,, Mg) must lie within or on the hexagon for safety. The lines upon which 
the slack functions vanish are shown. 


q,=0 (-oht) 


Fig, 3.2 
The objective of the present problem is to minimize the material in the face 
plates. This can be written as 
nk 
min v=4n| tr dr. (3.6) 
0 
The mathematical problem is thus to minimize V of (3.6) subject to the 
constraints (3.1, 2, 4, 5). 
The problem as formulated is in fact ‘irregular’ and allows degenerate 
solutions. To avoid the difficulties that this implies it is necessary to introduce 
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infinitesimal perturbations of the loading and of the constraints. The addition 
of a uniformly distributed load p(p > 0) over the whole plate adds pr?/2 tor 
and enables the assumption to be made that rQ cannot vanish over a finite 
interval. The second perturbation introduces a factor 1 + 7 (n > 0) on the right 
hand side of the second line of (3.4) to replace this by 


Mo + po = oht(1 +n), —Mp + q2 = oht(1 + 7). (3.7) 


Both p and 7 are retained as finite until the analysis of the conditions for a 
minimum has been completed and are then set equal to zero. 

Equations (3.4, 7) or Fig. 3.2 show that it is impossible for three or more 
slack functions to vanish simultaneously, unless ¢ = 0 and all the slack functions 
vanish. This implies M, = Mg = 0 and so by (3.1) rQ = 0. It follows that three 
or more slack functions cannot vanish on an interval. It also follows that ¢ 
cannot vanish on an interval. 

There are six ways in which pairs of slack functions can vanish, without the 
remainder vanishing. These correspond to the vertices of the hexagon of Fig. 3.2. 
It is convenient to label points or regions of the plate, to which the conditions 
at these vertices apply, as follows: 


(a) Pi =p2=0, (b) p2=43=0, (c)qi=q3=9, (d)qi =42= 


0, 
3.8 
(e) P3 =42=0, (f) pi =p3 =0. ib 


It will be assumed, similarly to the case of the beam in Section 2.1, that the 
interval (0, R) can be divided into a finite number of sub-intervals, such that the 
minimum solution of the present problem satisfies in each of them one of the 
conditions (a) to (f) of (3.8) or the condition that one or no slack function 
vanishes in that interval. 

Variations from the minimum solution must, by (3.1, 2, 4, 5, 7), satisfy 


d(r6M,)/dr —6Mp = 0, (3.9) 
6M, continuous, 5M,(R) = 0 (edge free or simply supported), (3.10) 


5M, + 6p, =ahdt, —5M, + 6q, = ohdt, 
5Mo + p2 = ohdt(1 +7), —6Mg + 5q2 = ohdt(1 +7), G.I) 
5M, —5Mg + 5p3=ahdt, —5M, + 5M + 5q3 = ohdt, 


5p, 20 (pj, =0), bq, 20(q,=0) = (i= 1,2, 3). (3.12) 
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It can now be shown that (3.9, 11) can be solved in ways that enable the follow- 
ing division into dependent and independent variations to be made: 


Region Dependent variations Independent variations 
(a) 5M,, 5Mg, St, 641, 642, dps, d43 51, 6p2 
(b) 5M,, 5Mo, St, pi, 841, 542, 5p3 5p2, 543 
(c) 5M,, 5Mo, 5t, 5p,, 5p2, 542, ps 841, 543 
(d) 5M,, Mg, 5t, pi, D2, dps, 543 841, d42 
(e) 5M,, 5Mg, St, dpi, 641, dp2, dq3 5p, 542 
() 5M,, 5Mo, 5t, 841, 6P2, 642, 843 5px, bps 


(3.13) 


As an example consider (b). The conditions p2 = q3 = 0 give, by (3.12), 

6p, 20, 6q3 2O and these are the only variations that need to be placed in 
the independent position, since in fact none of the other slack functions or the 
function ¢ vanish in (b). Equations (3.11) may be solved to give 


6p, =—5M, + ohbt, 6q1 = 6M, + ohSt, 
56My =-5p2+ohdt(1 +n), 5q2=5Mg + ohdt(1 +n), (3.14) 
6p3=-6M,+5Mg + ohdt, 5M,=ohndt— dp. + §q3, 
and so, substituting from (3.14) into (3.9), yield 
oh(nr d5t/dr — 6t) = rd 5p./dr — d(r5q3)/dr. (3.15) 


Equation (3.15) integrates in the form 


ét= (c +f [{rddp2/dr— d(r6q3)/dr} John!" **] ar) rm” “@ie) 


where C is a constant of integration. Equations (3.14, 16) show that 6M,, 5Mg, 
51, 5py, 5q1, 5¢2 and 6p; can be expressed in terms of 5p. and 5q3. 

The whole interval (0, R)+ can be split into regions (a) to (f), since for the 
present purposes a region in which a single slack function vanishes may be 
included in one of two of the regions (a) to (f) and a region in which no slack 
function vanishes may be included in any of them. This implies a broadening 
of the definition of, say, (a) in (3.8) to include regions in which two, one or 
none of p, and p> are zero, while the rest of the slack functions are non-zero 
there. It makes no difference to the conclusions of (3.13). With this interpreta- 
tion, each interval, into which (0, R) is divided by the minimum solution, has a 
solution like (3.14, 16) for the variations, and in particular for 5M,, which 


+ This interval should be taken as ‘open’ at the end r @ 0, since solutions with logarithmic 
infinities at r = 0 do occur, All the integrals involved can be shown to be convergent for 
the actual solutions. These difficulties can be avoided by replacing the lower limit of 
(3.6) by Ry and taking the limit RX, ~ 0 in the solutions obtained. 
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contains an arbitrary constant. These constants may be chosen to satisfy the 

conditions of (3.10), with one to spare when the edge is clamped. It is to be 

remarked that if 7 is set equal to zero in (3.15) one of these constants is lost. 
The condition for a minimum of V in (3.6) gives 


‘R 
{ dtrdr>0, (3.17) 
0 


which on introducing Lagrangian multipliers rA;, ru;(i = 1, 2, 3) and ¢, as well 
as the usual positive infinitesimal €, can be written using (3.9, 11) as 


R 
| [coher St + rd, (5M, + dp, — ohdt) + ru,(—6M, + 6q, — ohdt) + 
0 

+ Po {5Mo + 5p2 — ohSt(1 + n)} + rtp {-5Mp + 8q2 —ohSt(1 + n)}+ (3-18) 


+7d3(5M, — 5M + 5p3 — cht) + ru3 (5M, + 5Mg + 5q3 — ohSt) + 
+ o{d(r6M,)/dr —5Mg}] dr>0. 
Integration of (3.18) by parts, assuming ¢ to be continuous, gives 
R 
{ [an ters = My — (Az +2 (1 +0) — Ag — a} Ot + 


0 


+ (Ar —aa + As Hs — dg/4r)5M, + Xo He Aa tHs—G/r)5My + —«F-19) 


3 
+d (Xp; +549] rdr + [¢r5M,]¢ >0. 
i=1 


The Lagrangian multipliers must now be chosen so as to remove the 
dependent variations from (3.19). This requires 
Ar + ea + (Ag + M2) +) +A3 + M3 = €, 
Ai — Ha + A3 — M3 = dd/dr, (3.20) 
NoMa —A3 +H = O/r, 


Ng = Hy = Ma = M3 = Oin (a), Ay = A3 = M1 = M2 =O in (b), Ay = Ag =A3 = Ha = Oin (c), 
Ay = Ag =A3 = Hs = O in (d), Ay = Az = Wr = Ma =O in (€), Az = Mi = Me = 3 = Oin (f), 
(3.21) 
and 
WR) =0 (edger =R clamped); (0) = 0 (edger =R free or simply supported). 
(3.22) 
Equations (3.20, 21) may be solved in each of the regions (a) to (f). For 


example in region (e), equation (3.21) requires that Ay = Az = 41 = M3 = O and so 
from (3.20) 


Ag = dd/dr, Hy = —dd/dr ~ p/r, ndg/dr + (1 + n)¢/r =e. (3.23) 
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The last of (3.23) integrates to give 
$ =—er/(1 + 2n) + K/rO 7, (3.24) 


where K is a constant of integration. The solutions for all the regions contain 
such constants which may be chosen so that ¢ is continuous, as assumed above, 
and so that (3.22) is satisfied. It is to be remarked that if 7 = 0 the third of 
(3.23) ceases to be a differential equation and a constant K is lost. The problem 
is then irregular and the present argument breaks down. 

Applying (3.20, 21, 22) to (3.19) gives, using (3.10), 


| (Ardp1 + Adp2 )rdr + | (2 5p2 + U35q3)rdr +{ (41541 + M35qa)rdr + 
(a) (b) ke) 


t \ (u15q1 + M252 )rdr + \ (A35p3 + M2 5q2 )rdr + a Q4udpi + A36p3)rdr— 
(d) e) 
-(or8M,),_, 20. (3.25) 


It has been assumed in (3.25) that r6M, is finite as r > 0. This can be shown 

to be true for the solutions of (3.9, 11) in all the regions (a) to (f). The last 

term in (3.25) is, after (3.22), only required for the clamped edge case and for 
this, as was remarked above, the solution for 6M, contains an arbitrary constant. 
This means that the variation 5M,, as r > 0, can be taken as an independent 
variable in (3.25). The variations of the slack functions in (3.25) are also 
independent, as is shown by (3.13). Equation (3.25) can thus be dealt with 
term by term in the usual way, to give, remembering (3.12) and (3.21), 


A, = 0, where p; >0 and \; 20, where p; = 0 
(F=15203)" °6'20) 
4; =0, where gq; >O and yu; 20, where g; =0 


and 
(0) = 0, (3.27) 


valid now for all edge conditions. 

The necessary conditions for the minimum of (3.6) subject to (3.1, 2, 4, 5) 
have thus been shown to be the existence of multipliers \;, u;(i = 1, 2, 3) and ¢ 
which satisfy (3.20, 22, 26, 27), with now n written equal to zero. 


3.2 The dual problem 


The Lagrangian for the problem of (3.1, 2, 4, 5, 6) can be written, using the 
same multipliers as in (3.18), as 


R 
V* =(4n/che) | [ohert + rd, (M, + p1 — ht) + ruy(—M, + q, — oht) + 
0 


+1Q2(Mg + p2 — ht) + rua(—Mg + q2 — oht) + rA3(M, — Mg + p3 — oht) + 
+ ru3(-M, +My +43 — cht) + ¢{d(rM,)/dr — My -rQ} | dr, (3.28) 
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or on integrating by parts, assuming ¢ to be continuous, 


/ pR 
V*= (rlone( | {ance Ay Ha Ae — M2 Ags )t t+ 
0 
+ (Ar ~ Hr + A3 — Ms — dd/dr)M, + (Az — M2 — 3 + 3 — O/r)Mg + 
3 
+ DY Opi t wii) - oo} rdr + [orM,] : (3.29) 
i=1 
The conditions for the minimum of V* subject to (3.2, 5) and t > 0 are 
Art+MitAr tM tAztuz3=e (¢>0), 
AitMitAgth2tAz3tu3<e (¢=0), 
ArH +3 Ms = dd/dr, Az — M2 —A3 + M3 = @/r, (3.30) 
Ap, =0, A420, wq;=0, nw, 20 @=1, 2, 3), 
@R)=0 (clamped edge), 
where it has been assumed that rM, > 0, as r > 0. These agree with those 
obtained above except the condition for ¢ = 0, which was excluded in 
Section 3.1 and the absence of a condition on (0) from (3.30), except the 
implied finiteness. It will be seen that in actual solutions multipliers can be 
found which satisfy the second of (3.30) as an equality ¢ and that (3.30) does 
in fact imply $(0) = 0. 


The objective for the dual problem, which will be shown to be a maximum 
problem, is obtained by imposing (3.30) on (3.29). This gives 


R 
max W =—(41/ohe) | oQr dr, (3.31) 
0 


subject to, 
Aa ty + Ao t+ ho +A3 + M3 Se, 
Ar — Hi + Ag — M3 = do/dr, 
de — M2 —A3 + Us = O/r, (3.32) 
A420, 4,20 @ = 1,2, 3), 
@(R) = 0 (edge clamped), 


where ¢ and p;, q;(i = 1, 2, 3) have been eliminated from (3.30). The dual 
problem can be given a kinematic interpretation by taking ¢ as the 0-wise 
virtual component of rotation of the middle surface of the plate. The curvatures 


+ @ continuous is the ‘first corner condition’ for min V* (Pars 1962, Section 2.5). 
+ The values of these multipliers are not changed by taking the limit p > 0. 
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are then d@/dr and ¢/r. The objective W is easily seen to be proportional to the 
virtual work of the external forces and is to be maximized subject to restric- 
tions on the curvatures imposed by (3.32). These restrictions can be shown to 
be expressible as | d¢/dr |+| 4/r |+ | dd/dr + $/r |< 2e. 

The inequality between the W of (3.31), evaluated for a solution of the dual 
constraints (3.32), and the V of (3.6), evaluated for a solution of what can now 
be called the primal constraints (3.1, 2, 4, 5), can be established in a similar 
way to (2.44). 


W =~-(4n/ohe) ie ¢Ordr, (by 3.31), 
nag 
= (4nlone) | {Mp ~ (M,)/ar} dr, by (3.1), 
0 
= (4n/ohe) i {M,d@/dr + Mg¢/r} rdr, by (3.2) and last of (3.32), 
0 


R 
= (4n/che) | {M,Qu — t1 + X3— us) +t Mg Qa — ba - As +u3)}rdr, 
st by (3.32), 


R 3 3 
= (4nlane) | {ont Y Atud- Y pi +uiao} rdr, by (3.4), 
0 i=1 i=1 


R 
<4n( trdr, by (3.5, 32), 
0 


= V, by (3.6). (3.33) 


The usual deductions can now be made. Equation (3.33) places bounds on 
Vain. The dual problem is a maximum problem. If V has a minimum, then 
conditions (3.30) are satisfied by the solution and these define a dual solution 
with W = V, which must therefore give both the maximum to W and the 
minimum to V, The necessary and sufficient conditions that a solution of the 
primal constraints (3.1, 2, 4, 5) gives min V of (3.6) is that multipliers can be 
found which satisfy (3.30). 

The virtual deformation defined by ¢ may be thought of as a ‘collapse 
mechanism’ as in plastic design theory. The argument used to derive (3.33) then 
shows that the dissipation per unit area of the plate M,d@/dr + Mg ¢/r = ohet for 
the minimum volume solution given by (3.30). It follows that the dissipation 
per unit volume is a constant ohe/2. This is a special case of a general theorem 
given by Drucker and Shield (1956) and derived using the principles of plastic 
design theory. All the sufficient conditions of the Michell type like (1.70), (2.46), 
(3.30) and others like them, which have been deduced here using duality theory, 
are likewise special cases of this general theorem. The application of this theorem 
to the problems of this chapter is to be found in Onat, Schumann, and Shield 
(1957). More general loading cases are treated in Prager and Shield (1959). 
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3.3 Nature of the minimum solution 


The minimum solution cannot have a finite interval in which all the slack 
functions p;, q;(i= 1, 2, 3) are non-zero, since this would make, by (3.26), all 
the multipliers \;, u;(i= 1, 2, 3) vanish, in contradiction to the first of (3.20). 

The minimum solution cannot have a finite interval in which a single slack 
function vanishes, while the rest are non-zero. This means by (3.26) that all but 
one of the ),;, u;(i= 1, 2, 3) must be zero. Substitution in (3.20) then leads in 
all cases to a contradiction. For example, if A3 is the only non-zero multiplier, 
(3.20) gives 

A3 = d¢/dr = —¢/r =e, (3.34) 
which yields inconsistent formulae for ¢. 

The remaining possibilities are the six states (a) to (f) defined by (3.8) and the 
degenerate state in which all the slack functions and ¢ vanish. All of (a) to (f) are 
possible on an interval with the exception of (b) and (e). The multipliers for the 
state (e) are given by (3.23, 24), where the limit 7 > 0 can now be taken. The 
result is ¢ = —er and \3 =—e, which contradicts (3.26). A similar contradiction 
can be obtained for (b). 

The possible solutions valid on finite intervals thus reduce to (a), (c), (d), (f) 
and the degenerate case. The first four split into two pairs which can be obtained 
from one another by a reversal of the signs of M,, Mg and rQ. The three 
essentially independent regimes can thus be defined by 


(A) P1 =P2 = 0, rest > 0; (B) qi =q3 = 0, rest > 0; (C) p; = 4; = 0 (i= 1, 2, 3). 
(3.35) 


Equations (3.1, 4, 5) give, for intervals (a, @), (b, B) and (c, ¢) which have 
regimes (A), (B), and (C) respectively, the results: 


(A) M, = Mo = oht = M,(@) ch Qdr>0 @<r<a); 


(B) M, =—oht = (b/r)M,(b) + (1/r) \\ rQdr<0, Mg =0 (b<r<hb), (3,36) 
b 


(C) M,=Mo =t=0 (c<r<0). 


Equations (3.20, 26), with 7 = 0, give for the same intervals: 

(A) ¢=er/2+A/r, 4, =6/2—-A/r’, Ap = €/2 + A/r’, rest zero (a <r <a), 
with the constant A restricted by —ea?/2 <A <ea?/2; 

(B) ¢=—er+B, pw, =2€-B/r, wy =—e+ B/r, rest zero (bh <r <b), 
with the constant B restricted by eb <B < 2b; 

(C) As for (A) or as for (B) (e <r <2). 


(3,37) 
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The solutions given for (C) are not of course the only ones possible, but are in 
fact the only ones required. 

The regime (B) cannot apply to an interval ending at r = 0, since for this 
b = O and by (3.37) 5 = 0. If the regime (A) applies to an interval ending at 
r= 0, then a = 0 and by (3.37) A = 0, which removes the singularity at r = 0. 
The only regimes valid near r = 0 are (A) and (C) and for both these (0) = 0. 
Finally it is to be remarked that no single regime can give the solution to a 
clamped edge problem. The only real candidate is (A) and with A = 0, (3.37) 
gives $(R) = €R/2 #0. 


3.4 Examples 


In the case where the edge r = R is simply supported or free, solutions using (A) 
only are possible. Writing a = 0,@ = R in (3.36) gives, recalling (3.2), 


R 
M, = My = oht =~ { Qdr>0 (0<r<R). (3.38) 
r 


The conditions of (3.37) impose no restriction in this case. The solution (3.38) 
R 

is valid so long as i} OQ dr is negative (0<r<R). If it is positive then the sign 
r 


of the external loads can be changed. If it changes sign, two or more different 
regimes must be introduced (see clamped edge solutions below). 
The case of uniform loading Z per unit area covering (Ro <r <R;_) gives 


Q=0 (<r<Ro), 
Q=-Z(P-Ro*)/2r (Ro<r<R;), (3.39) 
Q=-Z(Ry?-Ro?)/2r (Ry <r<R). 
Substituting from (3.39) into (3.38) yields 
M, = Mg = oht = (Z/2) {(R,? —Ro?)/2 — Ro? log (Ri/Ro) 
+ (Ri? —Ro”)log(R/Ri)}  (O<r<Ro), 

M, = Mg = ht = (Z/2) {(R1* — 1° )/2 —Ro? log (Ri/P) (3.40) 
+ (Ry? ~Ro”) log (R/Ri)}_ (Ro <r SR), 

M, = Mg = oht =(Z/2)(R1? —Ro?) log (R/r) (Ri <r <R). 


The volume is most readily obtained from (3.31), using 6 = er/2 from (3.37), 
and (3,39). The result is 


V =(F/4oh)(2R? —Ro’? —R,”), where F=1Z(R,? — Ro’). (3.41) 


The special cases of uniform loading over the whole plate (Ro =0,R, =R) with 
V = FR?/4oh, concentrated loading at r=Ro(R, =Ro) with V = F(R? —Ro*)/20h 
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and concentrated loading at r= 0(R, =Ro > 0) with V = FR?/20h are readily 
derived. For this last case 


M, = Mg = oht = (F/27) log (R/r) =(0<r<R), (3.42) 


which gives an unbounded ¢ as r > 0. Note however that 7M, > 0 as assumed 
above in the derivation of the boundary conditions of (3.30). 

The case of a free plate uniformly loaded and supported at the centre by a 
force F is readily solved as above to give 


Q=—-F/)2nr + Fr/2nR?, 
M, = Mg = oht = (F/2n) {log (R/r) ~ (1 —P/R?)/2}, (3.43) 
V = FR? /4oh. 


In the case when the edge r = R is clamped solutions made up using (A) and 
(B) are possible. (A) is valid around r = 0 and gives, by (3.37), @ = er/2. (B) is 
valid near r = R and so (3.37) and (3.22) give ¢ = e(R —r) and b >R/2. 
Continuity of ¢ is possible only if b = 2R/3, which is >R/2 as required. 
Continuity of M, at r= b =@ requires by (3.36) that ¢ and M, should vanish at 
this point.t Equations (3.36) then give 


2R/3 
M, = Mo = ont =~ | Qdr>0 (0<r<2R/3), 


r 


(3.44) 


iF 


M,=-oht = (0/7) | rOdr<0, My=0 (2R/3<r<R). 


2R/3 


2R/3 
This solution is valid so long as | Qdr<0 (0<r<2R/3) and 


r 


is 
| rQdr<0O (2R/3<r<R). If the integrals are positive the sign of the 
2R/3 


external loads can be changed, but if they have opposite signs or have zeros, 
more complicated solutions have to be constructed. 


2R/3 
The special case for which | Qdr=0, which occurs when there is no 


loading inside r = 2R/3, can however be dealt with using (C) and (B) only. The 
regime (C) extends from r = 0 tor =R, > 2R/3, which is the radius at which 
loading begins. The multipliers are chosen to be the same as for the previous 
solution. This means (A) of (3.37) is valid 0 <r < 2R/3 and (B) of (3.37) for 
2R/3 <r<R. The multipliers for (C) 0 <r <R, are thus chosen as for (A) in 
part of the range and as for (B) in the rest. This gives a valid solution just as it 
did previously. The continuity of M, gives M,(R ,) = 0 in this present case. 


+ Since M, is continuous, ¢(b — 0) = =1(b + 0), with both limits > 0. 
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Equation (3.36) then gives 


M,=Mg=t=0 (0<r<R), 

: (3.45) 
M, = ~ohe = (1/7) { rQdr<0, Mg=0 (Ri <r<R), 

R, 


) 
which is valid when Q = 0(0 <r <R,), with Ry > 2R/3, and | rQ dr is 
negative. Ry 

The special case of uniform loading for which Q = —Zr/2 is solved by (3.44). 
The results are: 


M,=Mg = oht =(Z/4)(4R7/9-7?) = (0<r<2R/3), 
M, =—oht = —(Z/6r)(r? — 8R?/27), Mg =O =(2R/3 <r <R), (3.46) 
V = 19FR?/1620h, where F = 1R°Z. 
The case of concentrated loading f per unit length at r = R,, for which 
Q =—(fR,/r) 1 (7 — Rj) is solved by (3.44) or (3.45) according as Ry < 2R/3 or 
22R/3. The results are: 
Case Ry <2R/3 
M, = Mg = oht = fR; log(2R/3Ri) + (0<r<R;), 
M, = Mg = cht = fR, log(2R/3r) (Ry <r <2R/3), 
M, =—oht = —(fR,/r)(r — 2R/3), Mp =0 = (2R/3<r<R), 
V =(FR?/30h)(1 — 3R;7/2R”), where F = 2nRyf. 
Case R, > 2R/3 
M,=M,=oht=0 (0<r<R), 
M, =—oht =—-(fR,/r)(r-R1), Mg =0 (Ri <r<R), (3.48) 
V =(FR?/oh)(1 —R,/R)*, where F = 2nRyjf. 


(3.47) 


The case of a concentrated load F at the centre r = 0 follows from (3.47), by 
letting R, > 0, with F remaining finite. The result is 


M,. = Mg = oht = (F/2n) log (2R/3r) (0<r<2R/3), 
M, = —oht = -(F/2mr)(r —2R/3), Mg =O =(2R/3<r<R), } (3.49) 
V = FR?/3oh. 
The examples given above are all based upon plastic design. However those 


for simply supported or free edges, which are based upon solution (A), with the 
constant A = 0, are optimum elastic designs as well. The stresses in the faces are, 


by (3.38), given by M,/ht = Mo/ht = 0. The curvatures dd/dr and ¢/r are by (3.37) 


with A = 0, given by e/2 in both cases. The corresponding components of strain 
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in the faces are he/4, The design is elastic, with virtual strains equal to real 
strains, if € = 40(1 — v)/Eh, where E is Young’s modulus and v is Poisson’s ratio. 
These elastic designs are uniformly stressed and hence have a uniform density 
of strain energy. This means that they are also designs of maximum stiffness for 
given volume of material, since it can be shown that the criterion for this kind of 
optimum design is that the strain energy is uniformly distributed throughout 
the structure. The proof of this result for a sandwich plate is effectively the 
same as that for a plate loaded in its own plane. Reference may therefore be 
made to Section 5.9 below, where the problem of maximum stiffness for this 
latter structural form is dealt with. 


4 Michell's structural continua 


4.1 Michell’s sufficient conditions 


In Chapter | the optimum frameworks were selected from layouts containing a 
finite number of nodes and members. Complete generality requires that any 
point in space can be a node and any segment can be occupied by a member. 
However any attempt to erect a general theory, analogous to that of Section 1.1, 
using conditions of equilibrium and restrictions on stress levels, is clearly not a 
practical possibility. Generalization of the dual problem of (1.65, 66) does 
however lead to new developments, which were first expounded in the 
important paper Michell (1904), 

The problem of (1.65, 66) can be formulated in completely general terms 
as follows: 


Find that virtual displacement of space, subject to the kinematic conditions 
imposed on the frameworks and to the restriction that the overall strain for 

all segments lies between —0€/0¢ and ge€/dy, which is such that the virtual 

work of the given external forces is a maximum, (4.1) 


This is a problem of the calculus of variations and, as in all such problems, 
progress can only be made if the mathematical character of the unknown 
functions is restricted (Pars 1962, Section 1.3). Let is be assumed that the virtual 
strain of (4.1) is piece-wise continuous. This means that the Cartesian compo- 
nents of virtual displacement are continuously differentiable functionst of the 
Cartesian coordinates, except at a finite number of surfaces of discontinuity. 
This is the usual assumption made in the derivation of a theory of strain and 
leads in the well known way to the existence of principal strains and principal 
directions. The statement of (4.1) can now be altered to read: 


Find that virtual displacement of space, subject to the kinematic conditions 
imposed on the frameworks and to the restriction that the principal strains 
everywhere lie between —0€/0g and 0€/0p, which is such that the virtual 

work of the given external forces is a maximum. (4.2) 


This is valid since the condition of (4.1) implies that all direct strain components 
at all points of space lie between the limits —o€/a¢ and oe/a-p and since, con- 
versely, the condition of (4.2) implies that the overall strain for any segment lies 
between these same limits. 


+ Up to the third order at least, if a theory of compatibility of strain is to be developed. 
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If (4.2) has a solution, then the minimum volume required for a framework, 
which carries the given forces, may be expected to be given by the maximum 
virtual work divided by oe. The layout of the corresponding structure may be 
expected to be given by those segments, which have strains equal to —a€/a¢ or 
o¢/oy and define possible compression or tension members. This gives definite 
locations for members when the lines of principal strain are straight lines. In 
general however such lines are curved and so possible members can only lie along 
those infinitesimal segments into which the lines of principal strain may be 
imagined to be divided. The attempt to generalize the theory of Chapter 1, to 
allow for all possible layouts, has thus led to a theory of frameworks with 
curvilinear members, which, as equilibrium of individual members shows, must 
form a continuum. 

The transition from (4.1) to (4.2) requires the introduction of special assump- 
tions and so, even if (4.1) generates all optimum frameworks, the formulation of 
(4.2) may well only generate a sub-class. Sufficient conditions, corresponding 
to (4.2), that a framework should be an optimum, are easily established. There is, 
however, no assurance that such a structure will exist, for all possible systems of 
given forces. A large number of optimum structures have been found, neverthe- 
less, by the use of these conditions and so their study is an important part of 
the subject of optimum frameworks. 

Michell’s sufficient conditions, for a framework to have the least volume of 
material, can be formulated as follows: 


A pin-jointed framework has the least volume of material, if it can carry its 
given forces, with stresses in its tension members equal to Op and stresses in its 
compression members equal to —0g and if a virtual deformation of a region of 
space, in which the competing frameworks must lie, satisfies the kinematic 
conditions imposed on the framework and gives strains of o€/oy in its tension 
members, strains of —0€/0¢ in its compression members and has no direct 

strain lying outside these limits. (4.3) 


The formulation of (4.3) restricts the frameworks to a definite region of space. 
This is normally the case in practice, but the whole of space may be used if 
‘absolute optima’ are required. The theorem of (4.3) is a generalization of (1.70) 
at least as far as sufficient conditions are concerned. Structures which satisfy 

the conditions of (4.3) or similar conditions may be called ‘Michell structures’. 

A proof of (4.3) may be given for any kind of framework with discrete and/or 
continuous distributions of members. The class of all competing frameworks, 
which can safely carry the given forces and which includes the Michell frame- 
work, is subjected to the virtual deformation of (4.3). The virtual work W is the 
same for all the frameworks, since the forces and deformations are the same. 
The principle of virtual work gives 


W=d Tyl, (4.4) 


where 7 is a typical end load in a member, y the corresponding virtual strain 
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and / the length of the member. The ‘sum’ > sums over all members of a frame- 
work and must be replaced by appropriate integrals for continuous distributions 
of members. 

For the Michell framework 


Ty = (o7A )(oe/a7~)=aeA (tension members), 


| as 


Ty =(-0¢A)(—oe/0¢)=0€A (compression members), 


where A is the area of cross-section of the corresponding member. Equations 
(4.4, 5) give 
W=0e> Al=oeVy, (4.6) 


where V,, is the volume of the Michell framework. 
For any of the frameworks under consideration 


Ty <T(oe/oq)<oeA (tension members), 


Ty = (-T)(-y) S(-T)(ce/o¢) <0€A — (compression members), en) 
since —0cA <T < oA and —0e/0g <y <ce/op. Equations (4.4, 7) give 
W<oed Al=aeV, (4.8) 
where V is the volume of any framework which carries the given forces. 
Equations (4.6, 8) complete the proof by giving 
Vana SV. (4.9) 


A special case of (4.3) arises when a framework, whose given supports, if any, 
are not redundant, carries the given forces with all its members in tension with 
stress op or in compression with stress —¢g. Such a framework is an optimum, 
since it allows a uniform dilatation of space with linear strain € or —e, as the 
case may be. This result was deduced by Maxwell from a more general theorem 
on frameworks which are stressed to their limits (Maxwell, 1890). This theorem 
can be stated in the form 

Op Vey — Oc Ve = Wole, (4.10) 


where Vp is the volume of tension members, V¢ the volume of compression 
members and Wo the virtual work of the external forces for a uniform dilata- 
tion with linear strain e. The proof of (4.10) follows immediately from (4.4). 
If Vo (or Vp) is absent from (4.10), it follows that all frameworks, which carry 
given forces with their members fully stressed in tension (or compression), have 
equal volumes of material. Many interesting examples of single-stress structures 


and of the application of Maxwell’s theorem (4.10) are to be found in Cox (1965). 


A number of general conclusions can be drawn about the layouts of Michell 
frameworks and their associated strain fields from the properties of principal 
strains and directions. If a single tension (compression) member passes through 
a point then its direction has a principal strain o€/0p (—o€/0¢) and the principal 


Michell’s structural continua 73 


strains at right angles are limited in magnitude. If a pair of distinct tension 
(compression) members passes through a point then the direct strains in all 
directions in their plane are oe/o-p (—o€/o¢) and any number of coplanar tension 
(compression) members are possible; the principal strain at right angles to these 
members is limited in magnitude. If a pair of tension and compression members 
meet at a point, they must be orthogonal and have strains ge/oq and —oe/og. 
No other members can be coplanar with them and the principal strain at right 
angles must be limited in magnitude. If three non-coplanar tension (compression) 
members meet at a point, the direct strain in any direction at this point must be 
o€/0~(—Ge/0¢). If two tension (compression) members and one compression 
(tension) member meet at a point, then the compression (tension) member must 
be at right angles to the other two, the strain in any direction in the plane of the 
tension (compression) members must be ge/a-~ (—o€/0¢) and the strain in the 
direction of the compression (tension) member must be —ae/o¢ (ae/o7). 

Michell layouts in three-dimensions have a very special character. They have 
been little studied however. Most progress has been made in two-dimensions, 
with which this chapter is mainly concerned. 


4.2 Two-dimensional strain fields 


The case where the two principal strains are both oe/oq or —oe/d¢ gives no 
guidance on layout at all. All members are in tension or alternatively in 
compression and can be arranged in any way that is suitable for carrying the 
given forces. As was seen above, the volume of material required is independent 
of the layout. 

The case where one principal strain is o€/o7 and the other —ae/ag gives 
orthogonal layouts of members, which provide most of the known Michell 
frameworks. The analysis of this strain field is the immediate task of this 
section. 

A system of strain with principal strains o¢/o and —0€/0g will have a 
definite orthogonal system of lines of principal strain. These may be used to 
define a curvilinear coordinate system with coordinates (a, 8) (Fig. 4.1). This 
system will be taken as right-handed} with strains oe/oy in the direction of the 
a-lines and strains —oe/og in the direction of the B-lines. The line element ds 
will be given by an expression of the form 


ds? = A? da? + B? d6?, (4.11) 


where A and B are positive functions of a and B. If ¢ is the angle between the 
axis Ox of a Cartesian coordinate system and the tangent to an a-line in the 
direction of a increasing, then 


0x/Ada=cos >, dy/Ada=sin gd, dx/BdB=—sind, dy/BdB=cos@. (4.12) 


+ A rotation from the positive direction on an a-line to the positive direction on a p-line 
agrees in sense with that from Ox to Oy, 
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Por Bs) 
Fig. 4.1 
Integration of (4.12) yields 
(a, B) 
xtiy=xotino+| e'? (A da + iB dB), (4.13) 
(0, 0) 


where (xo, Yo) is the point at which a = B = 0. The condition of single-valuedness 
of (4.13) gives 


3¢/da=—8A/BdB, 39/08 = 0B/Ada, (4.14) 


which integrates in the form 


$= bo +" e {c 9A/B9B) do + (8B/A da) dp t, (4.15) 


(0, 0) 
where ¢o is the value of o when a= B = 0. The single-valuedness of (4.15) gives 
(8/d0)(8B/A da) + (0/08) (dA/BaB) = 0, (4.16) 
which all such functions A, B must satisfy. 
Let the displacement components along the coordinate lines in the directions 


of a and B increasing be denoted by u and v, The direct strains in the a and B 
directions are oe/ap and ~oe/d¢ respectively and the associated shear strain is 
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zero. Denoting the rotation by w, the standard formulae in curvilinear 
coordinates} give 

du/A da — vd9/A da = ce/op, dv/BOB + ud¢/BdG = —oe/a¢, 
dv/A da + du/BOp + udd/A da — vdd/BAB = 0, (4.17) 
0v/A da — du/BOB + ud¢/A da + vdG/B dB = 2w. 
Solution of (4.17) for the derivatives of u and v yields 
df{el? (u + iv)} = e'? {o€(A da/og — iB dB/o¢) + iw(A de + iB dB)}, 
(4.18) 
which integrates to give 
i(b,-9) -, (a, B) | 
utiv=e!%o~ (yy + ing) + et? | e'? {ae(A defor —iBdB/oc) + 


(0, 0) (4.19) 
+ iu(A da + iB dB)}, 


where uo, Vo are the displacements at a = 6 = 0. The single-valuedness of (4.19) 
requires that 
dw/da = oe(1/op + 1/0¢)0G/da, 9v2/0B =—ae(1/op + 1/0¢) 06/08, (4.20) 


which integrates to give 
(a, B) 
= wn oe(Nlon + 1/0c) { ~  {(@6/2a) da~(09/08)46}, (4.21) 
(0, 0) 


where Go is the value of «, when a = 6 = 0. The single-valuedness of (4.21) gives 
974/dad6 = 0, (4.22) 
which characterizes the layout of Michell frameworks. The result of (4.22) can 
also be written, using (4.14), as 
(0/0) (8B/A de) = (0/08)(8A/BAB) = 0, (4.23) 


which are equations for A and B. Equation (4.22) can be integrated over a 
curvilinear rectangle bounded by a = ap, @ = a,, 6 = Bo and B =f, to give 


$(a, Bo) — $(o, Bo) = (a1, B1) — (G0, Br). (4.24) 


This shows that the angle turned through by the tangent to an a-line (8 = Bo), as 
it moves between two fixed B-lines (a@ = @, a ), is equal to the angle obtained 
by using the a-line (8 = 8; ) instead or indeed any other a-line. The angles 
involved here are marked on Fig. 4.1. The same result is found if the roles of 

the a and £-lines are interchanged. This geometrical property is the characteristic 
of sets of orthogonal curves known as Hencky nets. Slip lines in two-dimensional 


+ See, e.g, Love (1927), Chap. 1, equations 36 and 38, 
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perfect plastic flow form such nets. A knowledge of slip-line fields is thus of 
great value for the discovery of Michell layouts and the techniques for calculating 
the form of slip lines can be used in this new application. Reference may be made 
to Hill (1950) Chap. VI and to Geiringer (1937) for an account of these 
techniques.+ 

The boundary conditions at a fixed line of support, along which u + iv = 0, 
may also be noted. Equation (4.18) gives 


w=to¢eV(oq0g), Ada=+V(o7/o¢)BdB. (4.25) 


The first of (4.25) determines w on the fixed line, while the second gives the 
angle between the a@ or §-lines and the tangent to that line. 

States of strain in which only one of the principal strains is at its limit are 
also possible and give layouts consisting of a single set of non-intersecting 
members. These must all be straight lines, unless suitable distributed forces act 
in the region concerned. Consideration will therefore be limited here to a 
system of straight lines, which are assumed to have strains oe/o7 and therefore 
to be alines. These lines will envelop an evolute and the §-lines will be involutes 
(Fig. 4.2). The coordinate 8 can be taken equal to ¢ and the a-coordinate can be 
taken as the distance along an a-line measured from a fixed involute. It follows 
that 


A=1, B=a+F(§), (4.26) 


where F(6) is an arbitrary function depending on the form of the evolute, which 
will have the equation a + F(8) = 0. Finally 


o=8, xtiy=xotivetae® +i f° e'® F(B)d8, (4.27) 


where the second result follows from (4.13, 26). 

The strain in the @-lines is o€/oq. Let the strain in the 6-lines be € . Then 
equations like (4.17) can be written, taking into account the special forms of 
(4.26, 27). This gives 

Bu/da=eloy, (8v/4B + u)/{a + F(B)} = er, 
dv/da + (Au/dB —v)/ {a + F(B)} = 0, (4.28) 
8v/8a — (u/dB — v)/ {a + F(B)} = 2. 


The first and third of (4.28) integrate to give 


u=ealoy + GB), v= GB) + {a+ F(6)} H(B), (4.29) 


+ This analogy with plasticity theory was noted by Prager (1958) and by Hemp (1958). 
Prager also gives a powerful method for the graphical construction of Hencky nets. 
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where G, H are arbitrary functions. Equation (4.28) then gives for w and €, 
the results: 


w =H), 
€2 = H’(B) + {ea+ G"(B) + F'(8) HB) + G()} | {a + F()}. 
This system of strain is valid for a Michell structural layout if 
—a€lOg S€2 <a€/oy. (4.31) 


(4.30) 


4.3 Conditions of equilibrium 


The structure corresponding to a strain field with principal strains oe/o7, 
—oe/d¢ will have, say, members along a-lines with ‘equivalent thickness’ or 
cross-sectional area per unit width f,, and members along 6-lines with 
equivalent thickness tf. The end load carried in the a-direction by an element 
with coordinate difference df will be opt,Bd§8 and that carried in the 
8-direction by an element with coordinate difference da will be —ogtA da. It 
is convenient to introduce 


T, =o7Bt,, T, =—ocAtz, (4.32) 


which are end loads per unit coordinate difference in the a and B directions 
respectively. 

The differential equations of equilibrium can be written down using standard 
results for curvilinear coordinates.+ These give 


aT, /da—T, 39/08 =0, AT,/08 + T,d¢/da=0, (4.33) 


in the absence of forces in the body of the structure. 

Boundary conditions at edges, which may carry ‘concentrated’ members of 
finite cross-sectional area, are easily established. For an edge lying along a 
B-line, with the structure on the side of a decreasing, which has a concentrated 
member carrying a compressive load P and a distribution of external forces of 
magnitudes F,, Fg per unit length in the a and directions respectively, it 
follows that for equilibrium 


T, = P09/08 + FB, aP/06 = FB. (4.34) 


The corresponding result for an edge along an a-line, with a member carrying a 
tension T, is 
T, = T0¢/8a+ FzA, 9T/8a=—F,A. (4.35) 


The values of T and P follow from the second of equations (4.34, 35), together 
with initial values of end load induced by a given concentrated force, which 
requires concentrated members for its transmission. 


+ See, e.g. Love (1927), Section 331. However (4.33) are easily established directly. 
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The sizes of the distributed members follow from (4.32). The ‘total 
equivalent thickness’ f is given by 


t=t, +t, =T,/o7~B-—Th/ocA. (4.36) 


The structure corresponding to a straight line system with strain oe/o7 and 
with the second principal strain not at either limit, will have t. = 0. Equation 
(4.33) gives T, as a function of 8 only, as it should, and so in this case T, 
follows from a boundary condition of the form of (4.34). Equations (4.26, 32) 
give 

t1 = Ti@)/oz {a + FH). (4.37) 


4.4 Mathematical considerations 


The simplest Hencky net consists of two sets of parallel straight lines which are 
orthogonal to each other. For these a, 8 can be taken as Cartesian coordinates 
with A = B = 1. The angle ¢ is constant and (4.13) gives the usual relation 
between two distinct Cartesian coordinate systems. Equation (4.21) gives 

«@ = constant and (4.19) gives 


U= Ug — Woh + Gealoy, V=V9 + Woa- GEP/dg, (4.38) 
which is a rigid body motion combined with strains o€/o7, —o€/o¢ in the 
coordinate directions. 

If one coordinate line, say an oline, is straight, then, by (4.22) or (4.24), it 
follows that all the o-lines are straight. These a-lines will in general envelop an 
evolute as in Fig. 4.2. Equations (4.26, 27) are then applicable. Equation (4.21) 


Involute 
Evolute 
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gives 
W = Wo — eB (1/o7 + 1/a¢), (4.39) 
and equation (4.19) yields 


u + iv = (Uo + vp )e“® + [oe/agq + i {29 — o€B(1/o7 + 1/ac)}Ja+ 
+e L {-ioe/og — wo + o€B(1/ay + 1/o¢)} Fe dB. (4.40) 
Equations (4.33) integrate in the form 
T, = To) iat T,(a)da, Ts = T3(Q). (4.41) 


where the functions Tjo(8) and 7>(a) follow from (4.34, 35). 
In the case where both coordinate lines are curved it is possible to define the 
coordinates so that 


$=¢9—aa+bf, a,b=+lor-l. (4.42) 


This requires « = —a (¢ — ¢9) on 6 = O and 8 = b(¢ — Go) on a = O. It is thus 
necessary to define a as the angle through which the tangent to 6 = 0 turns, as 
it moves from a = 0 in the direction of a increasing, and to define 8 as the 
corresponding angle on a = 0 (Fig. 4.3). The choice of the signs of a and b 
depends upon the direction of rotation of these tangents. In Fig. 4.3,a@=b=1, 


but three other combinations of sign will be obtained if the curvature of one 
or both of the lines 8 = 0, a = 0 is reversed. Equation (4.42) is a consequence 
of these definitions and the rule of (4.24). It is clear that this coordinate 
system will only be valid between points of inflexion on § = 0 and a = 0. 
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Equations (4.14) now give 
04/88 =aB, dB/da=bA, (4.43) 
and so A, B satisfy 
07A/dadB=cA, 37B/dadB=cB, c=ab. (4.44) 


Equations (4.44) can be integrated by Riemann’s method. A particular solution 
is given by 


A=I,(t), t=2V {c(a—4(G-n)}, (4.45) 
where J is the modified Bessel function, which satisfies 
d?Jo/dt? + d/,/tdt —Ip =0, Jo(0)=1, (4.46) 
and £, 7 are parameters. It follows that 


(8/8e) (A BIo/ 98) — (8/28) (Io A/Ae) = A (37Io/BaaB — clo) — 


—I9(8?A/aaaB-cA)=0, 47) 


and so by Green’s theorem that 
$ (Uo 3A/da) da + (A dIy/38)d8 } = 0. (4.48) 
Application of (4.48) to the circuit DO, CPD of Fig. 4.4, where CP is a 


B P (cn) 


Fig. 4.4 « 
B-line, DP an a-line and P the point (£, 7), gives, since (Jo )a=¢ = (lo )g=n = 1, 


AG.2) =Ap + { {(lo9A/8a) da + (AaFo/2B) 46} (4.49). 
Do,c 
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Equation (4.49) gives A at P, when A and 0A/da are known on DO,C. If DO,C 


is replaced by DOC, where OC is an a-line and OD a §-line, then (4.49) yields, 
after an integration by parts, 


A(E, n) = A olo {2V(cén)} + { (Io 3A/ax) dee + { (1994/08) dB. (4.50) 
oc oD 


Equation (4.50) gives A at P when 0A/da is known on OC and 34/08 = aB is 
known on OD. 

A formula for B can now be obtained from the first of (4.43). Equation (4.13) 
can be used to derive x and y. Equation (4.21) gives, using (4.42), 


@ = Wo — o€(1/a7 + 1/ag)(aa + bB), (4.51) 
and (4.19) can be integrated to give u and v. 
Equations (4.33), (4.42) give 
8T,/08 =aT,, 9T,/da=bTp, (4.52) 


which can be integrated to give T, and T, in the same way as A and B follow 
from (4.43). 


4.5 Straight lines and circles 

A particular case of the Hencky net, which consists of straight lines and 
involutes, arises when the evolute contracts to a point. This gives a set of 
concurrent straight lines and an orthogonal set of concentric circles. Such a 
figure clearly satisfies (4.24). It is convenient to use polar coordinates and to 
write 


a=r, B=8. (4.53) 
Equation (4.27) and (4.26), with F(§) = 0, give 
$=0, A=1, B=r. (4.54) 


Equations (4.39, 40), with wo = uo = Up = 0, give, on changing the notation for 
the displacements, 

@w=—oe(l/op+1/oc)@, u,=cerfoz, ug =—oe(1/o7 + 1/o¢)r6, (4.55) 
where it is to be remarked that w and ug are multivalued and so the virtual 
deformation of (4.55) cannot be applied to the whole plane. 

This layout of circles and radii can be used to solve the problem shown in 
Fig. 4.5.+ Here an optimum structure is sought to balance the forces shown and 


to lie inside one of the half-planes defined by the line through their points of 
application. The solution consists of a concentrated member of area F/ag at 


+ All the problems of this section were solved in Michell (1904). 
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r=land a complete fan of radii —7/2 <@ < 2/2. Equation (4.34), applied to 
r=, gives P= F and (T,),=; = P= F. The second of (4.33) and (4.35) gives 
T, =O and hence, by (4.32), t2 = 0. The first of (4.33) gives T, independent 
of a and therefore equal to F. Equation (4.32) then gives t; = F/apr. The fact 
that ¢, and f, are non-negative and that the member at r =/ is in compression 
shows that the conditions of (4.3) are satisfied. The total volume of material 
in the structure is 


Venin = TIF(1 Joy + 1/G¢), (4.56) 


which may also be found by applying (4.6). The displacement at each of the 
forces F is, by (4.55), equal to oe(1/a7 + 1/o¢)in/2 in the direction of these 
forces. The point of application of the force 2F has zero displacement in (4.55). 
It follows that W/ae is equal to Vain in (4.56). 

A Michell strain field, which is valid for the whole plane, can be constructed 
by joining two opposed systems of circles and radii of half-angle 99 < 1/4, with 
strains of opposite sign, to two systems of orthogonal straight lines and two 
regions of hydrostatic compression (Fig. 4.6). The displacements in the upper 
fan are given by (4.55); those in the lower fan can be obtained by an interchange 
of strain limits. The displacements in the region of orthogonal straight lines, 
referred to axes O(x, y) in Fig. 4.6, are given by 

u, = 0€x/oy + o€(1/oq + 1/oc)Oo¥, Uy =—GE(1/o_ + 1/d¢) 80x — aey/oc. 
(4.57) 
These follow from (4.38) with up = v9 = 0 and wo adjusted so as to give con- 
tinuity of displacements at 8 = 49 or y = 0. The displacements in the regions of 


hydrostatic compression, which are referred to polar coordinates (p, ¥) in 
Fig. 4.6 are given by 


Up =—GEeplog, uy =—oe(l/oz + 1/oe)Bop. (4.58) 
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The first equation of (4.58) gives a uniform dilatation of linear strain —o€/o¢ 
centred on O, while the second gives a rotation about O. Continuity with the 
displacements of (4.57) is satisfied at y = 0 or x = 0. The displacements of 
(4.58) match the displacement in the lower fan at 8 = —@9, which are given by 
(4.55) with the signs reversed and oy interchanged with oc. 

One interesting feature of the present strain field is the existence of lines 
whose points have no motion at right angles to the line of symmetry @ = 0. 
These are given by either 


(a) uz sin 89 t+uy cos89=O0 or (b) up cos(y +09) —uy sin (W + 69) =0. 
(4.59) 


Substitution from (4.57, 58) in (4.59) gives 
(a) x/y = {(o7 + 0¢)9o sin Bo — G7 cos Ao}/ {G7 + Fc) 0 cos Bo — de sin a} (4.60) 
(b) tan (Y + 90) = o7/(o7 + 9c), 

which define straight lines through O. If their direction is determined by x, the 
angle made with the normal to @ = 0, then since 


(a) x= 9o- tan” (x/y), (b) x= +o, (4.61) 
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equations (4.60) give 

(a) tan x = {oy cos? 99 — o¢ sin* Bo} (oy + Oc) {0 —(1/2) sin 26g.) 

(b) tan x = o7/(Gq + Gc) 80. J 
Relations (a) are valid for 7/4 > 69 > y, where, as is easily seen from (4.60), 

tan y= Oq/(G7 + O¢)7- (4.63) 
Relations (b) are valid for y > 49 > 0. For the special case when o7 = dg the 
value of y is 37-4°. 

The strain field of Fig. 4.6 may be used to construct solutions to the problems 
illustrated in Fig. 4.7. A force 2F acts along the perpendicular bisector of the 
segment of length 2/ joining two fixed points. The point of action can be at any 
distance h from this segment. The form of the optimum structure depends upon 
the angle x = tan”! (h/I), which must be identical with the x in (4.62), since the 
line (shown dotted in Fig. 4.7), whose points have no displacement normal to 
the line of symmetry of the structures, must pass through the fixed supports. 


(4.62) 


as 
(b) ale 
ee h 
2F ak 
(c) 
/ \ 
/ \ 
(a) Q=n/4, (b) n/4>%>7, (Cly>% 
Fig. 4.7 
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The value of 69 follows from (4.62) and the form of the structure depends 
upon whether @p is greater than or equal to ¥ or less than y. Examples of the 
first case are shown in Fig. 4.7 (a) and (b); the second case is shown in (c). 
Inspection of the structures of Fig. 4.7 shows that the signs of the stresses 
in the various members agree with the signs of the virtual strain field of 
Fig. 4.6. The condition of zero displacement at the supports parallel to the line 
of symmetry of the structures is not satisfied by (4.55, 57, 58), but the addition 
of a suitable translation parallel to that line of symmetry will bring these 
supports to rest. The conditions of (4.3) for an optimum structure are or can 
be satisfied. 
The volume of material in the structures is most easily obtained using (4.6), 
at least when 09 > y. The coordinates of a point of support are given by 


x =Isec x sin(@9—x), y =I sec x cos (@9 — x). (4.64) 


Equations (4.57, 64) then give the displacements of this point, in particular, 
the displacement parallel to 9 = 0 in Fig. 4.6, which is u, cos 9 — uy sin 49. 
This is the displacement of the force 2F, when the supports are taken to be at 
rest in the virtual deformation. It follows that 


Vinin = 2F (u, cos 09 — Uy sin O9)/a€, (4.65) 


where u,, Uy are to be taken from (4.57) and to be evaluated for values of 
x, y given by (4.64). The result, after some transformation, using the first of 
(4.62), is 


Venin = 2FI{(Oy + 5c)? 90° — F¢ Fe} /OyGe(Gr + Fc) {Bo — (1/2) sin 2004 66) 
(89>y or h/l<tany), ~~ 


which must be taken in conjunction with the first of (4.62), with tan x = h/I. 
The result for the case @9 <7 is easily obtained by direct calculation, which gives 


Venin = (2Fl/o¢) (l/h + h/l) (89 <y or h/l> tan y). (4.67) 

The figure which is obtained from the Hencky net of Fig. 4.6 by taking 
69 = 2/4 is of great importance in its own right. It is reproduced in Fig. 4.8. 
It can be used to generate optimum structures for systems of parallel forces, 
such as are illustrated in Fig. 4.9. These structures are all determinate and it is 
clear by inspection that the signs of the end loads and stresses agree with those 
of the virtual strains of Fig. 4.8. The volume of material required for cases (a) 
and (b) can be calculated by observing that the x, y coordinates of the point of 
application of a force F are 


x=(I-AyJ2,  y=(1+h)V/2, (4.68) 


and that the displacement of this point resolved in the direction of F is 
(ux —uy)/\/2. The volume is thus 


Veoin = V2Flux ~Uy oe, (4.69) 
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or 
SEER 


where u, and uy are given by (4.57) with 69 = 7/4 and x and y by (4.68). The 
resulting formula for Vin is thus 
Venin = (Fl/oq0¢){(Gr + O¢)(1 + 1/2) + (7 — o¢)h/l}. (4.70) 

The special case of (a) where A = 0 should be compared with (4.56). The 
structure of Fig. 4.9 (a) has a volume (1 + 7/2)/z times that of Fig. 4.5 and is 
thus some 18 per cent lighter, although it carries the same forces. The reason 
lies in the fact that the strain field upon which Fig. 4.5 is based is, unlike that 
of Fig. 4.8, not valid in the whole plane. Fig. 4.5 does not give an absolute 
optimum, but only that for the plane cut along a ray through the origin, so 
that the displacements are single valued. 

The volume of material for case (c) of Fig. 4.9 follows by observing that the 
coordinates of the point of application of a F are 


r=V(h? +P), 6=tan"'(I/h), (4.71) 


and that the displacement of this point is given by —u, cos 6 + ug sin 8, where 
u, and ug are given by (4.55), with their signs reversed and oy and ¢ inter- 
changed. The resulting formula for Vin is thus 


Vesin = (2FYOq0¢){(Oq + Oc) tan“! (Wh) + ogh/I}. (4.72) 


A three-dimensional strain field can be generated from Fig. 4.6 by a 
rotation about the line of symmetry @ = 0. The original two-dimensional strain 
field is imposed upon all the co-axial planes so obtained, to produce an axially 
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(b) 
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(c) 
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Fig. 4.9 
symmetric field. The only ‘new’ strain component induced by this process is a 
hoop strain €y given by 
€y = (u, sin 0 +ug cos @)/rsin6, (u, +uy)/(x+y) and 
{up cos (89 + ¥)—uy sin (9 + W)}/p cos (09 + Y), 
in the various regions of Fig. 4.6. Substitution from (4.55, 57, 58) gives 


(4.73) 


€y = 0€/oq — 0e(1/og + 1/o¢)8 cot @ and —ve/dg + oe(1/o7 + 1/o¢)8 cot 4, 
(4.74) 
€p = Ge[x {1/oq — (1/oq + 1/0¢)80}-¥ {log — (/oy + 1/o¢)Oo} (x +¥), 
(4.75) 
and 


€y = ~0€/d¢ + Ge(l/oz + 1/0) Oo tan (4 + W). (4.76) 
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It then follows, since 0 <4 <4 <n/4, x, y >0 and 0) <0) + WV <2/2- 4%, 
that 


—0€/0¢ S€y < 0€/o7, (4.77) 


for all the expressions (4.74, 75, 76). This three-dimensional strain field has its 
strain limited as is required by (4.3) and so can be used to generate optimum 
structures. The structures of Figs. 4.7 and 4.9 can be rotated about their lines 
of syrametry and any number of the resulting structures and loading systems 
can be combined together to produce optimum structures in three dimensions. 
Two simple cases are illustrated in Fig. 4.10. The solution for (a) is obtained by 


---F--So<e 
> 
7 


Fig. 4.10 


placing one of the figures from Fig. 4.9 in each vertical plane through the 
diagonals of the square formed by the points of application of the forces F. The 
solution for (b) is obtained in the same way using one of the figures from 

Fig. 4.7. The result of using Fig. 4.7 (c) is shown. 


4.6 Equiangular spirals 


An interesting strain field may be obtained from (4.43) by taking a = b = | and 
making A = B. This gives 


A=B=/2.R exp (a+), (4.78) 
where R is a constant. It is convenient to take $9 = —z/4 and then (4.42) gives 
o=—n/4-atB. (4.79) 
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Equation (4.13) gives, with x9 =R, Yo =0, 
x + iy =r exp (10) =R exp{a +8 + i(8—a)}, (4.80) 
where polar coordinates r, 6 have been introduced. Equation (4.80) yields 
r=Rexp(a+$), @=8-—a(mod 2n). (4.81) 
The coordinate lines have the following equations: 
elines r=R exp (26-4), 
Glines r=R exp (2a+6). vey 


These are orthogonal equiangular spirals, with angles 7/4, encircling the origin 
(Fig. 4.11). 


Fig. 4.11 


The displacements take a particularly simple form when o7 = dg. This 
assumption will therefore be made for the rest of the section. Equation (4.51) 
gives, with wo = —€ 


w =—e—2e(a +B). (4.83) 
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Equation (4.19), with uo = Vo = 0, gives 
u + iv=(u, + iug) exp i(@ —¢) = V2. eR(1 —i)(a@ +8) exp (a+), (4.84) 


where the polar components u,, ug have been introduced. Equations (4.84) and 
(4.81) then yield 


u,=0, ug =—2er log (r/R), (4.85) 


which shows that the virtual deformation consists of the rotation of concentric 
circles, with a circle of radius r rotating through an angle —2e log (r/R). It is to 
be remarked that the circle r = R has zero rotation and can be taken as a fixed 

support. 

The present strain field was applied by Michell to the problem of the 
cantilever (Fig. 4.12). An optimum structure is required to transmit the force 
F to the fixed circular boundary r = R. A solution is obtained by using the 
strain field of Fig. 4.11. 


The force F requires concentrated members to carry it. These members lie 
along the coordinate lines a, 8 = (1/2) log (i/R), as is shown by (4.81). They 
carry forces of magnitude F/s/2. Equations (4.34, 35) then give, recalling 
(4.79), 


T, =—F/V2 at B= (1/2) log (YR), Ty = F/A/2 at «= (1/2) log (/R). (4.86) 


These are the boundary conditions for (4.52), which, with a = b = 1, can be 
written 


OT, /08=T,, dT, /da=T2. (4.87) 
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It is convenient to change to new coordinates a, 8 given by 


&= (1/2) log (/R)—a, B= (1/2) log (/R) — 8. (4.88) 
Equations (4.86, 87) then become 
aT, /a8=-T,, 9T,/da=—-To, 
7 - — (4.89) 
with T, =-F/x/2 at B=0 and T, =F//2 ata=0. 


Equations (4.89) agree with (4.43) with a = b =—1 and T, written for A and 
T, for B. The required solution is then given by (4.50), which yields 


T2&, 0) = -(F/V2) fo {2-VE 1)} + \, (-F/V2) Io [20 {£6 — n)} 148 


0 

= —(F//2) Io 2VEn)} — FV2) VQalE) hs 2VED)}, (4.90) 
where use has been made of the relation /,'(t) + (1/t)i (©) = Jo(t), between 
the modified Bessel functions of the zero and first orders. 

The value of T, follows by symmetry and is 
T(E, n) = -T2(n, 8). (4.91) 
It follows by (4.32) that t,, t2 >O and from (4.36) that 
(E, 7) = (F/oR) exp {-€E + 0)} Uo 2VEn)} + E+ A LVED}/2VEN)) \ (4.92) 

where E=(1/2)log(/R)—£, 1 =(1/2) log (/R)~ 2. 


The total volume of material is best obtained using (4.6) and (4.85). The result is 
Venin = (2FI/o) log (l/R). (4.93) 
An interesting application of the present solution to the problem of trans- 
mitting a pure bending moment is given in H. S. Y. Chan (1972). The structural 
layout and the given moments M, which are applied by rigid circles of radius R, 
are shown in Fig. 4.13. The virtual strain field is defined separately in each of 
the three regions: (a)R <r <I/\/2, 0<6 <a; (b) x+y <i, y 20, less the 
region (a); (c)x + >1, x <1. The deformation in the rest of the plane 
follows by symmetry. The displacement field is given by 


(a) up=0, ug =er {1 +2 log (I/2r)}, 
(b) u,=0, ug =er, (4.94) 
ic) uy =-e(I-x), uy =e(l-9). 

The field defined by (4.94) (a) differs from that defined by (4.85) by a rigid 


body rotation about the origin r = 0. It thus has the same strain field as Fig. 4.11, 
with 67 = 0¢ = 2. The field of (4.94) (b) is a rigid body rotation of amount € 
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7ex,, 


Seiecae pans os 
Fig. 4.13 


about the origin and so has zero strain everywhere. The field of (4.94) (c) has 
lines of principal strain with strain € parallel to Ox and lines of principal strain 
with strain —e parallel to Oy. 

The displacements are continuous at the boundary of regions (a) and (b), 
where r = I/,/2. At the boundary between (b) and (c), the displacements of (b) 
have the values u, =—ug sin 8 =—ey and uy = ug cos 6 = ex. On this boundary 
x +y =/and so continuity with the displacements of (c) is satisfied. The 
symmetry conditions are clearly satisfied by (a) and (b) at y = 0, since (4.94) 
(a) and (b) are valid for y < 0 as well. The symmetry condition at x = /, namely 
u, = 0, is also satisfied. 

The strain field defined by (4.94) is thus a continuous Michell strain field and 
so the layout of Fig. 4.13 satisfies the kinematic conditions for an optimum 
structure. It remains only to show that it can carry the moments M with stresses 
agreeing in sign with the corresponding strains. The loads in the tie and strut at 
the centre line are M/l. These require concentrated members, which lie along +e 
lines and enter the cantilevers of equiangular spirals in the directions of te 
strain lines respectively. The boundary conditions of (4.86) would require 
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modification, if detailed analysis was needed, but it is clear that a consistent 
solution for member sizes can be found.} The structure of Fig. 4.13 is thus an 
optimum structure, with a volume of material given by (4.6) and (4.94) as 


Veni = (2M/o) {1 + 2 In (I/\/2R)}. (4.95) 


4.7 Michell’s sphere 


The strain field of Fig. 4.11 may be used to construct the optimum structure 

for the transmission of a uniformly distributed torque at an inner circular 
boundary to a similarly uniformly distributed torque at a concentric outer 
circular boundary (Fig. 4.14). Michell extended this solution to three-dimensions 
and showed that the optimum structure, which balances two opposing torques 
applied to two rigid co-axial circular discs, consists of a framework of members 
lying along the rhumb linest on a spherical surface, which passes through the 
edges of the discs. Fig. 4.14 may be regarded as giving an end-view of this 
framework, showing one disc with its applied torque and the balancing forces 
from the second hemisphere. 


It is convenient to use spherical polar coordinates (r, 0, 6) with @ = 0, 7 on 
the axis of the discs. Let the spherical surface have radius R and let the edges 
of the discs be defined by r=R, 0 =@9 and 7-49. 


+ The integral in (4.90) is absent here and F/,/2 is replaced by M/L 
+ Lines which bisect the angles between the meridians and parallels in any spherical polar 
coordinate system. 
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A spherical shell, with middle surface at r = R, can be in equilibrium, when 
carrying only a shear stress resultant S, associated with the 8 and ¢ directions, 
and depending upon @ alone. The conditions of equilibrium? give 


dS/d0 + 2S cot @=0, (4.96) 

which integrates to give 

S=So sin? 99/sin? 6, (4.97) 
where So is the input of shear per unit length at the edges of the discs. The 
shear stress resultant S is equivalent to two principal stress resultants +S lying 
along the rhumb lines and requiring end load carrying material of equivalent 
thicknesses t, = t2 = S/o, where the simplifying assumption that op = 0¢ =o 
has been made. The required total thickness f, is thus given by 


t=t, +f, =2S/a=2So sin? 40/0 sin? 0.= (4.98) 
Consider now a virtual deformation defined by 
u,=0, ug =0, ug =2er sin @ log tan (8/2). (4.99) 


The corresponding components of strain can be calculated using the standard 
formulae.§ The only component which does not vanish is the shear strain Yoo: 
which has the value 2e. The corresponding principal strains are +€ and so the 
Strain system defined by (4.99) is limited in the way required by Michell’s 
conditions. These principal strains lie along the rhumb lines of all spheres 
centred at r = 0, including of course r = R. The signs of the principal strains 
agree with those of the principal stresses corresponding to S, since the sign 
convention for stresses and strains is chosen so this is true. Michell’s sphere is 
thus the optimum structure for the transmission of torque. 

The volume of material follows in the usual way from the virtual work. The 
axial rotation corresponding to ug is ug/r sin 8. If the torque is T, the work is 
2Tu,/r sin @ evaluated at 8 = 7 — 09. Hence 


Vain = (47/0) log cot (89/2). (4.100) 
The relation between T and So is 
T = 2nR? So sin? 8. (4.101) 
4.8 Cycloids 
A particular solution to (4.43) with a = —1, b = 1 is given by 
A=2hcos(a+$), B=2hsin (a+), (4.102) 


where / is a constant. The corresponding formula for ¢ is given by (4.42). Taking 
o = 0 yields 

o=atZ. (4.103) 
* Love (1927), Section 344. 


+ This is the same as would be required for a shell with allowable stress in shear = 0/2. 
§ Love (1927), Chap. 1, Section 22. 
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Equation (4.13) now gives, with xo = ¥o = 0, on substituting from (4.102, 103), 
the formulae: 


x= (h/2){sin (a +8) + 2(@-B)}, y= (h/2) {1 —cos 2(@ + ff, 
(4.104) 


which shows that the coordinate curves are cycloids (Fig. 4.15).+ The effective 
coverage of this coordinate system is the strip 0<y <h. A vanishes for 

+8 =n/2 and B for a+ B =O and the coordinate curves envelop these loci, 
which by (4.104) are seen to be the lines y = h and y = 0 respectively. 


Fig. 4.15 


The virtual displacements may now be calculated. Assuming for simplicity 
that op = dg = 0, (4.51), with wo = 7/2, gives 


w = 7e/2 + 2e(a— B), (4.105) 
and so (4.19), with uo = Yo = 0, integrates to give 
u,, t+ iuy =(u + ive =eh {((a—B + 7/4)e%#%* +a0+B—a/4 + 


+i@-pe-p+ary. 41° 


Equation (4.106) gives the Cartesian components (u,, uy) and (4.104) may 
then be used to express these in terms of (x, y). The result is 


u, =—ney/2+ eh((1/2) cos"! (1 — 2y/h) 
+ bx/h-V{Q/h)(1—y/W)}C — 2y/h)), } (4.107) 
ly = nex/2 + eh {(x/h)? — (y/h\(1 —y/h)}, 


where the inverse cosine is to be taken between 0 and 7. It is to be remarked 
that both u, and wy vanish at x = 0, y = 0 and at x =0, y =h. 
Equations (4.52) have a solution, when a =—1,5 = 1, of the form 


T, =2F cos(a+B), Tz =—2F sin (a+). (4.108) 


+ This is the Hencky net used by Prandtl to analyse the plastic deformation of a strip under 
transverse pressure. See Hill (1950), VIII S. 
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Equations (4.32), with op = 0¢ = a, then give 
t, =(F/oh) cot (a+B), — t, = (F/oh) tan (a + B), (4.109) 


which are positive 0 <a + 6 < n/2 and thus show that the signs of (4.108) agree 
with those of the virtual strains of (4.107). They become unbounded at, 
respectively, y = 0 and h. This is an indication that concentrated members are 
required along these lines and, for consistency with the strains, a tension 
member must lie along y = 0 and a compression member along y = h. 

The Cartesian stress resultants T, and S,, can be calculated from (4.108) 
in the usual way by considering an infinitesimal right-angled triangle, with 
hypotenuse of length dy parallel to Oy. This gives 


T, dy = T, dB cos $+ Tz dasin¢g, S,, dy=T, dB sin ¢—T, da cos ¢, 
(4.110) 


where da and df are related so that dx = 0. Substitution from (4.103, 108) 
gives, making use of (4.104), 


T, = 2F db/ {(x/h)(1 —y/h)f/dy, Sy = Fh. (4.111) 


The constant shear flow S,, requires balancing at y = 0, h. This will be assumed 
to be done by concentrated members lying along these lines. 

The resultant force components and couple corresponding to (4.111) are 
given by 


A 
{70-2 [Soar and (, teoyay= nen. (4.112) 
0 o o 


This indicates the loads that must be carried by the optimum structure, which 
carries the stresses (4.108) and ‘allows’ the strain field (4.107). 

The problem to which the present investigation will be applied is the 
optimum design of the ‘shear bracing’ of a long cantilever of length />h and 
uniform depth h under a tip shear F (Fig. 4.16). The unwanted bending 


Fig. 4.16 


moment of (4.112) is cancelled out by forces xF/4 applied to the flanges at 
the tip. External forces are provided at the root to balance the loads in the 
web. The ‘errors’ in the applied loading will only affect the volume of the 
structure by an amount of order Fh/o, which is small compared to the web 
volume of order Fi/o. 
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Michell’s conditions of (4.3) are satisfied everywhere except in the flanges 
up to a distance 7h/4 from the tip. This discrepancy will produce an error of 
order Feh in the virtual work or Fh/o in the volume, which can be neglected. 
The required volume is given by 

= = =h 
oe Main = (Telus) 24 +S ey [uy ]27"} dy + (F/A) [(ue), 222. 


(4.113) 
Substituting from (4.107, 111) in (4.113) then yields 


Venin = FI? /oh + 1Fl/20. (4.114) 


The first term of (4.114) is the amount of material required for the flanges of 
the girder irrespective of the nature of the shear bracing. The second term 
gives the material required for the shear bracing itself. In the present solution 
that shear bracing consists of a continuous array of orthogonal cycloids 
(Fig. 4.15). 

The optimum shear bracing, which consists of straight members joining flange 
to flange, may be shown to be given by the Warren truss, in which members are 
inclined at an angle of 7/4 to the flanges. The volume for the Warren truss is 


V=FI?/oh+2Fl/o (Warren truss), (4.115) 


and this has 4/7 times the material in the shear bracing as that of (4.114). 
A closer approximation to the layout of cycloids is shown in Fig. 4.17. This 
has a volume given by 
V=FI?/oh+V/3Fl/o (Truss of Fig. 4.17), (4.116) 


with a volume for the shear bracing intermediate between that of (4.114) and 
(4.115). 


439 Fields developed from circles 


It is possible to develop a Hencky net beginning with two orthogonal circles 
as starting curves 8 and a = 0 (Fig. 4.18).+ Taking a = b = 1 and $9 = 0 in (4.42) 


gives o=-a+t8. (4.117), 


+ This Hencky net has been much used in plastic flow investigations. See for example Hill 
(1950) Chap. VIL. 


x 
Fig. 4.18 
The boundary conditions for A and B are 


A=Ron$=0, B=Rona=0, (4.118) 
and with these (4.50) integrates, just as (4.90) did, to give 


A(@, 8) = B(B, &) = R [Io {2V(a8)} + VB/a) 11 2V(o8)}]. (4.119) 
Equation (4.13) may be integrated in terms of Lommel’s functions (Chan, 
1967), but the results are somewhat involved. Attention will be limited here to 
points for which x = y = L/\/2 and a= = y. For this (4.13), with xo = yo = 0, 
gives 


L= (RiV2) | {lot + h(O}dt = (R/V2) {Lo(2u) — 1 +2 5: (1)"lon +1 2w}, 
0 n=0 


(4.120)7 
from which values of L/R for various u can be obtained. 
The virtual strain field can now be determined. The assumption will again 
be made that op = 6¢ = o and so the principal strains are +€. The Hencky net 
generated by the circles 8 and a = 0 is extended backwards, for negative values 
of a and 8, by means of orthogonal circles and radii, to the centres of these 


+ The standard results /,(t) = /,(0), J, (0) + Jy, 4.1 (0 = 2/;,(t) (m7 = 1, 2...) have been used here 
t Hencky nets may be constructed graphically. See Hill (1950) Chap. V1, 5. 
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circles. The remaining triangular region in Fig. 4.18 is filled in with orthogonal 
straight lines. The boundary segment along x + y + R = 0 is taken as a fixed 
support. This is consistent with (4.25), which requires, in this case, the 
coordinate curves to meet it at an angle of 7/4. 

Equations (4.19, 21) can now be applied. These are valid for the whole 
region occupied by the Hencky net, including the fans and the triangle, 
because (4.18, 20) can be integrated along any path to give (4.19, 21) in spite 
of discontinuities in A and B, since ¢, (u, v) and w are continuous everywhere. 
Equation (4.25) gives w =—e on the fixed segment, since dx = —dy on the 
boundary, and (4.21) gives w = wo =—€ in the triangle. Equation (4.51) then 
gives 

w=—e-2e(a + B) (4.121) 


in the main field a, 8 > 0. The values of uo = €R, Vp — ER follow from the 
extension of the sides of the triangle. Equation (4.19) then gives, using (4.121), 


ut iv=eR[{l —i + 2(@-i8)} fo{2V(@8)} + 2(1 — 1) (aS) 1 2-V(a8)} I. 
(4.122) 


The simplest proof of (4.122) is by direct verification.+ 

The strain field of Fig. 4.18 may be applied to solve the problem of the 
optimum cantilever supported on a fixed segment for any system of loads which 
give tension members along a-lines and compression members along 6-lines. For 
the special case of a force F applied at a = 8 = p parallel to the rigid line of 
support, the required volume is, by (4.122), given by 


Vein = (/2FR/o) {(1 + 22) Lo(2u) + 2u 1 (2u)}- (4.123) 


The concentrated members required to carry F are shown in Fig. 4.18 and the 
required structure is enclosed within them. The circles of the fans and the 
orthogonal straight lines are not required in this case. 

The strain field of Fig. 4.18 may be extended beyond a = m/2 and 8 = n/2 
by just increasing the starting circles. However the layout so obtained will not 
be valid for the whole plane. It will break down for the same reason as does the 
complete array of concentric circles and radii, namely, discontinuity of displace- 
ment. An extension of Fig. 4.18 to the whole plane is possible, nevertheless, 
and a quadrant of this is shown in Fig. 4.19. The arc AB is a starting circle with 
centre C. The region ABC is a right-angled fan of circles and radii and the 
region ACD one of orthogonal straight lines. The arc AB and its ‘image’ 
generates the region ABE, just as in Fig. 4.18. The straight segment AD and the 
curve AE generate a Hencky net in ADFE consisting of straight lines and 
involutes. This can be obtained using (4.26, 27), since B and hence F(8) are 
known on the starting curve AE. The curve DF and its ‘image’ can now be used 


+ Equation (4.122) was first given in A. S. L. Chan (1960) where he applied it to the 
problem of the cantilever treated in this section, 
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Fig. 4.19 

to develop a Hencky net in a region, part of which, namely DFHI, is shown in 
Fig. 4.19.4 The region EFG must be filled in with orthogonal straight lines and 
a region FGJH of straight lines and involutes can then be generated. The curve GJ 
and its ‘image’ now generates a new part of the field. The way in which this 
process of extension of the strain field over increasing regions of the plane can 
be carried out is now clear, as is the fact that it has no limit. 

This extended field of Fig. 4.19 has been used by H. S. Y. Chan in his D.Phil. 
thesis to generate optimum structures for forces parallel or normal to a fixed 


+ This can be generated using (4.50), but the boundary conditions are now more 
complicated than (4.118). 
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Fig. 4.20 


Fig. 4.21 4 


line of support and applied at any point of the plane. He has also, in H. S. Y. Chan 
(1967), solved the same problem for structures restricted to lie in a half-plane. 
Examples of these important developments are shown in Figs. 4.20 and 4.21. It is 
to be remarked that, in Fig. 4.21, one set of layout lines are tangential to the 
boundary of the half-plane; this may be compared with the layout of cycloids 

in Fig. 4.15, where the structure is confined to a finite strip. 


5 | Plates loaded in their planes 


5.1 Formulation of the problem 


The present problem is to be interpreted as the design of plates of variable 
thickness t occupying a given region R of a plane referred to axes O(x, y). The 
plates are loaded in R by forces with components (X, Y) per unit area and on 
a part By of the boundary of R by forces with components (X,,, Y,,) per unit 
length. The rest of the boundary Bg is rigidly constrained. 

The state of stress in the plate is specified by stress resultants 7,, T, and S, 
with reference to the axes O(x, y) and by principal stress resultants T, and 7>, 
where 7, is inclined at an angle ¢ to Ox. The stress distribution is in general 
continuous, but discontinuities at a finite number of lines are possible. The two 
specifications of stress are related by 


T, =(T, + Tz)/2 + {(T, — Tz)/2} cos 26, 
T, = (T, + Tz)/2 — {(T, — T2)/2} cos 2¢, (5.1) 
S ={(T, — Tz)/2} sin 2¢. 


Equilibrium requires that 


8T,/dx + OS/ay+X=0,  aS/dx + aT,/ay + Y=0 (5.2) 
in R and, if (7, m) is the unit normal to the boundary, 
IT, +mS=X,, iS+mT,=Y,, (5.3) 


on the boundary, where (X,,, Y,,) are given forces on By and can be interpreted 
as forces of constraint on Bg. A final condition of equilibrium is that 


IT, +mS, IS+mT, are continuous at lines of discontinuity, (5.4) 
where (/, m) is now a unit normal to the line of discontinuity. 
Let it be assumed that the safety of the plate is governed by Tresca’s 
criterion, with a tensile yield stress 0. This can be expressed, in a similar manner 


to (3.4), as 
T, + py = ot, -T, +4, = 91, 


T, + p2 = ot, —T, + q2 = ot, (5.5) 


T,—-T,+p3=ot, —T,+T,+q3=01, 
where 


Pie Vi» Pre G2» Ps, Qs #0. (5.6) 
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Equations (5.5, 6) can also be represented by means of a hexagon in the manner 
of Fig. 3.2. The zeros of the slack functions p;, q; (i= 1, 2, 3), which correspond 
to the sides of the hexagon, define a set of points at which yielding is just beginning. 
It will be assumed, for the optimum structures studied here, that this set consists, 
for each slack variable, of a finite number of regions of R. This is the natural 
generalization to two-dimensions of the characteristic property of one-dimensional 
optimum structures found in Chapters 2 and 3. 

The optimum designs are chosen so as to make the volume V of material in 
the plate a minimum. This is expressed by 


min v= ff tdx dy. (5.7) 


The mathematical problem of optimization is thus to choose the functions T,, 
T,, S, T;, Ta, %, t, Pir 4; (i = 1, 2, 3) so as to make V of (5.7) a minimum, subject 
to the constraints (5.1, 2, 3, 4, 5, 6). This is a problem of the calculus of 
variations, but here the unknown functions are functions of the two variables 
x,y. 


5.2 Formal derivation of necessary conditions 
Astrict derivation of the necessary conditions for the minimum of (5.7), using 
the methods of Sections 2.1, 2.5, 2.6 and 3.1, runs into difficulties when trying 
to show that the equations for the variations and the equations for the 
Lagrangian multipliers can actually be solved. This is because the solution of 
boundary value problems for partial differential equations is a much more 
difficult matter than that for ordinary differential equations. The present 
section is devoted to purely formal processes, which are carried out as an act of 
faith. Fortunately they lead to a dual problem and to sufficient conditions for 
the minimum. 

The Lagrangian for the problem of Section 5.1 can be written 


V*= (Joc) {f {oer + (Ti tps —01) + KT +41 — ot) #Ro(Ts +Pa—o0)* 
R 


+ H2(-Tz + q2 — ot) + A3(T, — Tz + ps — of) + u3(-T, + Tr + Q3 — 08) + 
+ u(8T,/dx + AS/ay + X) + v(BS/ax + AT, /ay + Y)} dx dy — 


“(1/00 | {us (IT, + mS ~ Xj) + 04(1S + mT, ~ ¥q)} ds, (5.8) 
Br 


where T,, T,, and S are given by (5.1), € is the usual positive infinitesimal, 
Moxy), Oe, ¥) = 1, 2, 3), ue, »), UO»), 41 (5) and v;(s) are Lagrangian 
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multipliers and s is the arc length of the boundary. Application of Green’s 
theorem, recalling (5.4), gives 


{I {u(@T,,/ax + 8S/dy) + v(AS/ax + AT,/ay)} dx dy = 
R 
7 {u(IT, + mS) + v(IS + mT,)} ds + 
Br+Be 
+f {Au(IT, + mS) + Av(IS + mT,)} ds — 
D 
s } {| {T.du/ax + T, dv/ay + $(@v/ax + du/ay)}dxdy, (5.9) 
R 


where I is taken over the lines of discontinuity and Au, Av are the increments 
D 


of u, v on these lines. Substitution from (5.9) into (5.8) and the introduction of 
(5.1) gives 


v*=(1/0e) {j (ree + MCT + Ps 00) + wy(Ty +44 08) + 2a(Ts + Pa o1)+ 
R 


+ uo (-Ly + G2 — ot) + A3(T; — T, + p3 — ot) + w3(-T, + T, + q3 — ot) — 
— ((T, + Tz)/2 + {(T, — Tz)/2} cos 24] (8u/ax) — 
— (Ti + T2)/2 — {(1y — T)/2} cos 26] (Bv/ay) — 

[{(T; - T)/2} sin 26] — (@v/ax + au/ay) + uX + oY) dx dy + 


+ (a/oe) f {u—u,)(IT, + mS) + (v- 0% )US + mT,) + u1X,, + 01 ¥, }ds+ 
Br 

+ (I/oe) I {u(IT, + mS) + o(1S + mT, )} ds + 
Be 


+(1/oe) i {Au(IT,, + mS) + Av(IS + mT, }} ds. (5.10) 
D 


The conditions that V * of (5.10) should be a minimum can now be derived 
in the usual way. Remembering (5.6) and that ¢ > 0, then gives, in R, that 


3 3 
YAtu=e€ @>O0, YF Atud<e (¢=0), 
i=1 i=1 
Ay Ha tA3 Ha =€1, A2—H2—As tus =€2, (M1 -T.)7=9, (5.11) 
A =O0(p,> 0), uj, =O(qi, > 0), 4% O(p;=9), 4, > 0 (G; = 0) 
(= 1, 2,3), 
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where 
€1 =(1/2)(8u/dx + dv/dy) + (1/2)(Gu/dx — dv/dy) cos 26 + 
+ (1/2)(dv/4x + du/dy) sin 20, 
€) = (1/2)(du/ax + v/dy) — (1/2)(Gu/4x — Bv/dy) cos 26 — (5.12) 
—(1/2)(dv/dx + du/dy) sin 26, 
¥ =—(1/2)(8u/4x — dv/dy) sin 2g + (1/2)(dv/dx + du/ay) cos 29. 
Further conditions are, on the boundary, 
u,=u, v=vo0n Be, u=v=0 on Be, (5.13) 


and on the lines of discontinuity, u, v are continuous functions 
(5.14) 

The conditions of (5.11, 12, 13, 14) may be given a kinematic interpretation. 
The functions u, v can be regarded as Cartesian virtual displacement components. 
Equations (5.12) then show that €,, €2 and y are virtual strain components 
associated with the directions ¢ and ¢ + 7/2 of the principal stress resultants 
T, and T,. Equations (5.13) show that (u;, v, ) are the boundary values of 
(u, v) and that the displacements vanish on the constrained part of the 
boundary Be. Equation (5.14) implies that the displacements are continuous 
everywhere. Finally (5.11) imposes restrictions on the virtual strains. If 
T, # Tz then y = 0 and €;, €2 are principal strains with principal directions 
coinciding with those of the principal stress resultants. If 7, = Tz, then, by (5.1), 
¢ is not required for the specification of the stress resultants and so may be 
chosen to make yy = 0 under these conditions as well. The remaining equations 
of (5.11) limit the magnitude of these principal strains. 


5.3 The dual problem 


The problem of Section 5.1 may now be termed the primal problem. The dual 
is then obtained by imposing (5.11, 12, 13, 14) on to V * of (5.10) to obtain, 
as will be seen, a maximum problem 


max W = (I/ae) {J (Xu+ Yo)dxdy+ [ (Xqut Y,0)as} (5.15) 
R Br 


This depends upon the dual ‘variables’ u, v only and so, in formulating the dual 
constraints, the primal ‘variables’ ¢, @, T;, Tz, p; and q;(i = 1, 2, 3) should be 
eliminated from (5.11, 12, 13, 14). This gives, taking y = 0 everywhere, 


3 
y At) <e, Ay Ha tA3 M3 =€1*, X22 Az ts = €2*) 
(5.16) 
yw 20 (=1,2,3), u=v=Oo0nBe, u,v continuous in R, 


where e€,* and €,* are the quantities, which are obtained from (5.12) by writing 
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y = Oand eliminating ¢, i.e. the principal strains corresponding to the displace- 
ments u, v. The dual problem is thus to maximize W of (5.15), subject to the 
constraints of (5.16). 

Let T;, To, 9, T,, Ty, S, t, pz and q;(i = 1, 2, 3) be any solution of the primal 
constraints (5.1, 2, 3, 4, 5, 6) and let uw, v, A; and u;(i = 1, 2, 3)* be any solution 
of the dual constraints (5.16), then, beginning with (5.15), 


W= (ioe) {| cu + Yo)dxdy + I {u(IT, + mS) + v(ls + mt, has |, 
eke by (5.3, 16) 


=(1/ce) i} {a(uT, + uS)/ax + A(uS + vT,)/ay + Xu + Yu} dxdy, 
e by Green’s theorem and (5.4), 


=(I/oe) ii} {T, du/ax + T, dv]ay + S(dv/ax + du/ay)} dxdy, by (5.2), 
R 


=(1/oe) ff (T, * €,* + T,* €,*) dxdy, (5.17) 
R 


where 7,*, T:* and S* are the stress resultants referred to the principal directions 
of strain, for which (€,*, €2*, 0) are the corresponding components of strain. The 
invariance of T,, du/dx + T, dv/ay + S(dv/4x + du/dy) has been used in the last 
line of (5.17). 

Tresca’s conditions of (5.5, 6) imply that all shear stresses in the plate are 
not greater than 9/2. It thus follows that | 7,*|/2, | T>*|/2¢ and 
({(7,* —1,*)/2}° + S*?]"?/e cannot exceed this stress. Slack functions p;*, 
q;* (i= 1, 2, 3) can thus be found such that 


T,* + p,* = ot, -T,* + q,* =ot, 
T,* + p2* = ot, —T,* + q2* = ot, 
T,*—T,* + p3*=ot, -T,*+7,* +q3* = ot, 6.18) 
wits Pi*, 4°20 (f=1,2,3). 
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Returning to (5.17) and substituting from (5.16) gives 


W= (Joe) ff {Qu — 1 Xs — Hs) + Tea — a — Aa + Hs)} axed, 
R 


= aloe) {f Qu Ti" + (HTi*) + Xo Te* + ua Ta") + 3(T* - Ta") + 
- +us(-T,* + Tr*)} dxdy, 


3 3 
= (aloe {{ {or 5 stu) ¥ Oupi* tua} dxdy, by (5.18), 
R i=1 i=1 


< {j tdxdy, by (5.16, 18), 
R 


=V, by (5.7). (5.19) 


It follows by (5.19) that W is bounded above and that, if the primal has a 
solution Vin, then 


W<Vain <V, (5.20) 


and so feasible solutions to the primal and dual constraints give upper and lower 
bounds to Vain. If these solutions also satisfy (5.11), which means that 


3 
t > Qitus) =te, 
i=1 


€:=e,*, €2=€2*, (T1-T)y=9, (5.21) 
MPi=0, 4q;=0 (@=1,2,3), 


it follows that the principal directions of strain coincide with the principal 
directions of the stress resultants or, in the case of isotropic strain, can be taken 
to coincide and so T,* = T, and T,* = T;. Equations (5.5, 18) then give p;* = p; 
and q;* = q; (= 1, 2, 3). Reference to the third step of (5.19) and to (5.21) 
shows that W = V, which, by (5.20), must equal V,,;,- The dual and the primal 
both have solutions under the present assumptions and 


Wax = Venin- (5.22) 


The present considerations lead to sufficient conditions for an optimum 
design. They are expressed mathematically by (5.1—6), i.e. the constraints of the 
primal, and (5.11—14), with the understanding that y = 0 everywhere, even when 
T, = Tz. If the theorem is true that the necessary conditions for the Lagrangian 
V * of (5.8) to be a minimum also give the necessary conditions that V of (5.7) 
should be a minimum, when subjected to the constraints (5.1—6), then the 
conditions (5.1—6, 11—14) are necessary as well as sufficient for an optimum 
design. 
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5.4 Nature of the optimum solution 


Much progress can often be made using sufficient conditions, as Chapter 4 has 
shown, and the same is true here, since the present solutions have much in 
common with those for Michell continua. 

The optimum solution cannot have a finite region in which p; > 0, 

q; > 0 (i = 1, 2, 3). Equation (5.5) shows that ¢ > 0, while (5.11) gives 
A; = 4; = 0 (= 1, 2, 3) and hence r = 0. 

There are six possible cases for which one of p;, g; (i= 1, 2, 3) is equal to 
zero, while the rest are positive. Equation (5.11) then implies that all but one 
of d;, u; (i = 1, 2, 3) are zero and that the positive one is equal to €, since 
t > 0 for this case. The principal strains €,, €2 are also constant, and take the 
values +e or 0. In detail (5.5, 11) give: 


(a) T, =ot, O0<T,<ot, €; =€, €2=0, 

(b) T, =—ot, -ot <T, <0, €, =—€, €2 =0, 

(c) Tz, =at, O<T, <ot, €, =0, €2=€, 

(d) T, =—ot, —ot <<T, <0, €, =0, €, =—€, 

(e) 7; —T,=o0t, O<T, <ot, -ot <T, <0, €, =€, €2 =—€, 
(f) 7, —T,=~—ot, -ot <T, <0, O<T, <ot, €, =—€, €2 =€. 


(5.23) 


There are six further cases that can be obtained by setting two of 
Ps 4; (i= 1, 2, 3) equal to zero, while the rest are positive. These correspond 
to the corners of the hexagon of Fig. 3.2. Equations (5.5, 11) give for these: 


(g) T, =ot>0, T, =0, €, =€, —€ <2 <0, 
(h) T, =—-ot <0, T, =0, €, =—€, OS <e, 
(i) T, =0, T, =ot >0, —e<e, <0, €2 =€, 


(5.24) 
(j) T,=0, T; =-ot <0, 0<e, <e, € =—€, 
(k) T,; =T2 = ot >0, €, +2 =€, €, 20, €2 >0, 
() T, = Tz =—ot <0, €; t+€, =—€, & <0, e, <0. 
The degenerate case 
(m) T, = Tz =t=0, je1| + feo| + ler + €2| <2e (5.25) 


must be recorded for completeness. The stated restriction on strain can be 
derived from (5.11), remembering that p; =q; = 0 (i = 1, 2, 3) for this case. 
All regions of an optimum solution must satisfy 


T,€, + Tr€2 = o€t. (5.26) 


This follows from (5.23—25) and was proved, in passing, when (5.22) was 
established. If €;, €) are thought of as ‘rates of strain’, then since the states of 
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stress corresponding to 7, and T> in (5.23—25) all permit yielding, equation 
(5.26) means that the dissipation per unit volume is a constant (Drucker and 
Shield, 1956). The strains of (5.23, 24) satisfy the flow rules of plasticity theory, 
since the vector (€;, €2) is normal to the corresponding side of the ‘yield 
hexagon’ of Fig. 3.2 in the case of (5.23) and lies between the normals to 
adjacent sides in (5.24). 

A complete optimum solution will consist of a finite number of regions of 
types (a) to (m) completely filling the region R. Condition (5.3) must be 
satisfied on the boundary By, condition (5.4) at the boundaries between 
regions, condition (5.13) on Be and condition (5.14) everywhere. Suitable 
forms for T;, Tz, ¢, t, u and v must be obtained by integration of (5.1, 2) and 
(5.12), with y = 0, together with the relevant equations from (5.23, 24).+ The 
disposable functions thus obtained must be found from the boundary conditions 
listed above. In each region the inequalities of (5.23, 24) and those implied for 
(m) of (5.25) must be satisfied, before the solution can be regarded as complete. 
The difficulty of the present problem lies in fact that the various regions and the 
type of solution applicable to each cannot be determined in advance. The only 
method is to try appropriate solutions, which are possible at the various parts 
of the boundary and attempt to fit them together, with perhaps further solutions 
to completely fill the region R. A detailed knowledge of known solutions of the 
various types is an obvious prerequisite for this process. 


5.5 Virtual-strain fields 


It is convenient here, as in Section 4.2, to use the lines of principal virtual 
strain as coordinate lines for a system of curvilinear coordinates (a, 8). The 
purely geometric developments of equations (4.1 1—16) are thus available for 
application to the present problem. 

Equations (4.17) require modification, so as to deal with the slightly more 
general strain conditions of (5.23). It is sufficient to replace o€/oy by €, and 
~a¢/0¢ by €,, where €, and €2 are constants, which must be given the values 
required by (5.23) for regions of the plate satisfying one of the specifications 
(a) to (f). Equations (4.19, 21) then become 


_, p(B) 
utiv=e' uo tivo) + 2 | e® {4e, da + iBe, dp ie 
©. 4 iw (Adar iBag}, ©?” 
and 


‘a, 8) 
@ = Wo + (€1 — €2) {(2¢/8a) da — (46/98) 4B}, (5.28) 


«0, 0) 


+ This is a general statement and must not be taken too literally. The analysis will usually 
be carried out in appropriate curvilinear coordinates (see Chapter 4). 


110 Plates loaded in their planes 


where it should be noted, that here, as in Chapter 4, (u, v) are component dis- 
placements in the curvilinear coordinate system (along the coordinate lines) and 
not just Cartesian components, as in Sections 5.2—4 above. 

Equations (4.22—24) follow from (5.28) just as they did from (4.21). This 
means that, for the cases governed by (5.23), the lines of principal strain form 
Hencky nets. The considerations of Section 4.4 and the many examples of 
Sections 4.5—9 are thus applicable to cases (a) to (f) of (5.23). In the special 
cases (e) and (f) even the displacements are given by the same formulae, when 
O7 =8¢ =9. 

The boundary conditions on constrained parts of the boundary (5.13) give, 
by (5.27), a replacement for (4.25) in the form 


w=tV(-e1€2), A V(e1) da= +B-V(-e2) 8 on Be. (5.29) 


The cases governed by (g) to (j) of (5.24) as typified by (g), have been dealt 
with already in equations (4.26—31), under the special assumption, which is all 
that will be considered here, that the a-lines are straight. The only modifications 
required are to write a = 0 and to replace (4.31) by one of the rather more 
restrictive inequality systems of (5.24). 

The kinematic conditions of (k) and (1) of (5.24) impose no restrictions on 
the layout of the principal strain lines. Also since 7; = T2 no restriction arises 
from the statics either. This means that for simplicity it is permissible to take 
(a, 8) in (k) or (1) as rectangular Cartesian coordinates, with A = B = 1 and 
= o. The conditions for (k) then give 


du/da+ dv/aB=e,  dv/8a + Ou/dB = 0, (5.30) 
which integrate to give 
u=ea/2+F,(at+8)+F,(a-B), v=eB/2 — F,(a+ 8) + F(@—8), (5.31) 
where F, and F; are arbitrary functions. Equation (5.24) requires €,, €2 > O and 
so F, and F, must satisfy 
—e/2 <F,' (a+ 8) + Fy'(a— 8) <e/2. (5.32) 


The case (m) of (5.25) is even less determinate, since only inequalities are 
imposed. The strain field for this degenerate case can of course be taken as a 
strain field for a non-degenerate case, if this will give a continuous deformation. 


5.6 Conditions of equilibrium 


The equations of equilibrium (4.33) may be applied to the present problem, 
when body forces (X, Y) are absent, by writing 


T,=BT,, T=AT. (5.33) 
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Equations (4.34, 35) may be applied here as boundary conditions, with or 
without the concentrated members. Concentrated forces may be conveniently 
dealt with by introducing concentrated members, rather than singularities of 
plate thickness. 

The required plate thicknesses may be calculated from (5.23, 24). Introducing 
T, and T, from (5.33) gives 


cases (a), (b),(g)and(h) =| Ti|/oB, 
cases (c), (d), (i) and (j) t=|T,|/cA, (5.34) 
cases (e) and (f) t=|T,/B-T,/Al/o. 


It is to be remarked that case (e) agrees with (4.36), with a7 = 0c = 0, and that 
case (g) agrees with (4.37), with op =, since, by (4.26) B= a+ F(6). It follows, 
as should be expected, for cases (e), (f), (g), (h), (i) and (j), that the ‘equivalent 
thickness’ of a Michell continuum is equal to the thickness of the optimum plate, 
when the ‘loading’ 7, and T> is the same. Cases (a), (b). (c) and (d) do not arise 
for Michell continua. The plate has the advantage that, when stressed to the 

limit o in one sense, it can carry stresses of the same sign in a direction at right 
angles up to the limit o, without any additional material being required. 

The fact that, for cases (g), (h), (i) and (j), equilibrium, in the absence of 
body forces, requires one set of coordinate lines to be straight, follows from 
(4.33) as before. 

Application of (4.33), (5.33) and (4.14) to (k) and (1) of (5.24) gives 


cases (k) and (1) =|7;\/¢=|7,|/o = constant (5.35) 


which is to be expected, since a hydrostatic system of stress cannot vary in the 
absence of body forces. The magnitude of ¢ in (5.35) is determined by the 
boundary conditions. The Michell continuum is not to be recommended for 
cases (k) and (1). For 7; = 72, t) =t2 = |T\/or and so the equivalent thickness 
is 2\7;|/o7 or twice the t of (5.35), when oy = o. The advantage of plates for 
like signed principal stresses has been noticed above.+ 

Finally it may be remarked that integrals of the equilibrium equations (4.33), 
valid for special coordinate systems, like (4.41) and special forms of those 
equations, like (4.52), can be applied to plate problems. 


5.7 Examples 


It is now clear that any Michell continuum, which is based on a strain field 
€, = 0€/0-p, €2 = —0€/¢ and which has a stress distribution for which T, and 


+ This point was made in Richards and Chan (1966). 
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T> have corresponding signs, can be converted into an optimum plate for the 
same loading system by writing 07 = 0g = 0. Such examples give the case (e) 
of (5.23) or (g) of (5.24), with in the last case €2 at its lower limit —e. Cases 
(f) of (5.23) and (h), (i) and (j) of (5.24) can be realized by changing the sign 
of the loads and/or interchanging the role of a and 6. 

Fig. 4.5 may be interpreted as a plate bounded by a concentrated member, 
with the layout lines giving lines of principal stress. The thickness ¢ is equal to 
Ffor and by (4.56) Vain = 201F/o. 

The strain field of Fig. 4.6 is not valid for an optimum plate, since 
€, = €2 =—€ contradicts (5.24) (1) and indeed all the strain conditions of 
(5.23—25). However, if 89 = 7/4, the region of hydrostatic strain disappears 
and the resulting field of Fig. 4.8 is satisfactory. Fig. 4.7(a) gives an optimum 
design in which the fan may be replaced by a plate and the same may be stated 
for the designs of Fig. 4.9. The volume for cases (a) and (b) of Fig. 4.9 is by 
(4.70) given by Vain = (1 + 2)FI/o. 

The field of Fig. 4.11 is a valid strain field and the design of Fig. 4.12 may 
be interpreted as a plate with edge members. The required plate thickness for 
the optimum is given by (4.92) and the total volume by (4.93). The design of 
Fig. 4.13 for the transmission of moments can also be interpreted as a plate 
design. It is interesting to note that region (b) has €, = €2 = 0 and thus provides 
an example of (5.25). The volume of material for this last is given by (4.95). 
Fig. 4.14 also provides an example of a pure plate design based on the field of 
Fig. 4.11. 

The field of cycloids of Fig. 4.15 is also satisfactory for an optimum plate 
and the design of Fig. 4.16 with a web of thickness f = ft, + f2, as given by 
(4.109), is much more practical than a spider’s web of cycloids. Equation (4.114) 
is valid for the plate design. 

The designs of Figs. 4.18, 4.20, and 4.21 and of others like them can all be 
regarded as plate designs, with the same volumes of material as that given, in 
for example (4.123), for the corresponding Michell continua.+ 

The almost trivial example of a region R of any form, simply or multiply 
connected, loaded by a uniform normal pressure along its boundary, may also 
be noted. This is solved by (5.24) (I) with (say) €, = €2 =—e/2. Equation (5.35) 
shows that the plate thickness is uniform in the optimum design for this 
loading case. 

The examples given above are derived from a theory of optimum plastic 
design. They also provide examples of optimum elastic design if their virtual 
strain system can be their actual strain system. This is a severe requirement and 
is only met by two of the examples, namely that of Fig. 4.14 for the trans- 
mission of torque and the uniform plate under uniform pressure at its edges. 


+ Numerical examples of the designs of Figs. 4.18 and 4.12 are given in Richards and 
Chan (1966). 
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Sufficient conditions can be derived from (5.23, 24) in the form: 
(a) T,/t=9, T,/t=vo, €=o(1— YE, 
(e) 7,/t =—T2/t = 0/2, €=a(1 + v)/2E, 
(g) T,/t=0, T,/t=0, e=o/E, €, =—ve, 
(k) T,/t =T2/t =0, €; =€2 =a(1 —v)/E, €=20(1 —v)/E, 


(5.36) 


where E is Young’s modulus and p is Poisson’s ratio. The principal stresses are 
constant in each region. The structure for torque is an example of (e) and that 
for pressure of (k), with its signs reversed, i.e. of (I). 


5.8 Optimum design with the Mises—Hencky criterion 


If the Mises—Hencky criterion is adopted for defining the condition for the 
safety of the plate, then (5.5, 6) must be replaced by 


(T,2 + T,? — TT, + 387)? +p=ot, p2o. (5.37) 
The Lagrangian of (5.8) then becomes 
v*= (1/oe) | [oet + A{(T,? + T,? —T,T, + 3S7)'? + p—ot} + 
R 


+ u(OT,,/dx + AS/day + X) + (AS/ax + aT, /dy + Y)] dxdy — 
-(/oe) f {u, (IT, + mS — X,,) + o,(IS + mT, — Y,,)} ds. (5.38) 
Br 


The transformation of (5.9) can be applied once more to yield, instead of 
(5.10), the equation 


V* =(1/ce) ff [oet + A(T? +T,? -T,T, + 3521/2 +p-ot} = 
— T,€, —Tj€y —Sy + uX +vY] dxdy + 


(Ula) [ {ums )( + mS) + (0-04 )US + mT,) + m1 Xq +04 Yo} s+ 
Br 
+ (/0e) | {u(IT, + mS) + v(S +mT,)}ds + 

Bo 


+(I/oe) {Au(IT, + mS) + Av(IS + mT,)} ds, (5.39) 
D 
where, 
€, = du/dx, €, =duv/dy, y= du/dx + du/ay, (5.40) 


and if (u, v) is taken as a virtual displacement in the coordinate system (x, y) then 
(5.40) defines the usual strain components. 
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The conditions for V * to be a minimum can now be written as 


A=e (t>0), A<e (t=O), (5.41) 
6. = MT. —T [20 450,F — Ty +387) 2, 
€y = MT, ~ Tel2M(T? + T,? - Te, + 387)", ) 
¥ = 30S? + T,? - TT, +3587)", 

X=0 (p>0), ASO (P=), (5.43) 


u, =u, v,=vonBp, u=v=O0onBe andu, vcontinuous. (5.44) 
Equations (5.42) imply that 
€,? te? tee, 7/4 = 37/4 (5.45) 
and 
T.e,+T,e, +Sy=NT + Ty —T1, +357)". (5.46) 
Equations (5.42) are the flow rules for perfect plastic flow, if €,,€, and y are 
taken as rates of strain and so the virtual deformation can be thought of as a 
collapse mechanism. The dissipation per unit area is given by (5.46), with 
\ =e by (5.41) and the square root replaced by of from (5.37) with p = 0, 
which is required by (5.43). The dissipation per unit volume is thus a constant 
ge as is required by the condition of Drucker and Shield (1956). 


The dual problem can be obtained by imposing (5.41—44) on (5.39). This 
gives using (5.46) 


max W= (1/09) i} {. (Xu + Yo) dxdy + i (X,ut Y,2)4s}, (5.47) 
FF 


subject to 


(€,7 + €,? + €,€, + 77/4)? <V/3e/2, 


(5.48) 
u, v continuous and u = v= 0 on Bc, 


where (5.45) has been used to eliminate the stress resultants. 

The usual inequality between W and V, when calculated from solutions of 
the dual and primal constraints, can now be established. A beginning is made 
with 


W = (1/ae) i (T,€, + Ty €y + Sy)dxdy, (5.49) 
R 


which follows from (5.17, 40). The next step requires the inequality 
Treg + Tey + Sy <(2/V3){(Te? + Ty? ~ Te Ty + 38?) (6,2 + 6,2 + exe, + P/4}Y?, 
(5.50) 
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which is best established by introducing the ‘deviatoric components’ of the stress 
resultant and strains 


T,! = 2T,/3-T,|3, Ty’ =2T,/3—Tz/3, Tz =-T./3-T,/3, S'=S, 


' i re Fam 
Oy Fees Cy Flys Ce Fey Mey, YF = 7- 


Jesn 


Equations (5.51) imply 
T,'€,' + Ty'€,' + T,'€,' + 2S'(y'/2) = T,€, + Ty €y + Sy, 
TT, ty? +28 =O) tT — 27, t355)), (6:52) 
€,? + ey? +€,'? + 2(7'/2)? = 2(€,? +e,” + €,€, + 77/4). 


Application of Cauchy’s inequalityy to the left hand sides of the equations of 
(5.52), then leads by (5.52) to (5.50). Finally use of (5.50), (5.37) and (5.48) 
gives by (5.49) 


W <(l/ce) ite (2//3) ot (\/3e/2) dxdy = V. (5.53) 


In the special case when the solutions used to calculate W and V in (5.53) 
satisfy (5.41—44), it follows by (5.49), (5.46) and (5.37) that 


W =(1/ce) {J X(ot — p) dxdy = V, (5.54) 
R 


on using (5.41, 43). All the usual deductions can now be made. 

The solution of problems using the equations derived in this section is clearly 
not an easy matter. However, progress can be made by using (5.53) to establish 
bounds to an exact solution. For example let it be supposed that a solution to a 
problem using Tresca’s criterion is known and let the minimum volume be V4q. 
If a system of stress satisfies (5.5, 6) then it must satisfy (5.37). This follows 
from the graphical representation of (5.37), which gives as a boundary the 
circumscribed ellipse of the Tresca hexagon of Fig. 3.2. It follows that Va, 
is an upper bound to the minimum volume Wen for the Mises—Hencky 
criterion. Comparison of the constraints of (5.16), or better the equivalent 
result of (5.25), with (5.48) can also be made graphically. Equation (5.25) 
gives an hexagonal boundary, which can be inscribed in the ellipse, obtained 
from (5.48) by transforming to principal directions,t by replacing € in (5.25) 
by V3e/2. This reduces the virtual work of the external forces for the Tresca 
solution by a factor \/3/2. Dividing this reduced virtual work, whose correspond- 
ing strains satisfy (5.48), by ae gives the lower bound (/3/2) Vo . Collecting 


+ The scalar product of two vectors is less than or equal to the product of their magnitudes, 
which is true for any number of dimensions. 
t The strain condition of (5.48) is an invariant of the strain components. 
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these results gives 
(32) Voie < Va < Vins (5.55) 


which determines acti to within 7 per cent. 


5.9 The problem of maximum stiffness 


The general results obtained for other structures can be extended to plates 
loaded in their own plane, if the problem is formulated as that which seeks 
minimum stored strain energy for a given volume of structure. This can be 

written in the present case as 


min U= i [(1/2Et) {(T, + T,)? + 2 + v)(S?-7,T,)} | dxdy, 
R 


subject to, (5.56) 


{ t dxdy = V = constant, 
R 


where £ is Young’s modulus and v is Poisson’s ratio. 

A variation from f to ¢ + 5t induces corresponding changes 57,, etc., in the 
stress resultants due to the given external forces. These varied stress resultants 
are in equilibrium with the external forces and so the variations 57,, etc., form 
an internal system of stress. It follows therefore that, if 57,, etc., are imposed 
on U of (5.56), with t held unvaried, the resulting variation 5U will be zero 
(Castigliano’s principle). 

Varying t in (5.56) thus gives 


ii} [81/2EP){(T, + T,)? + 2(1 + »)(S? - T_T, } | axdy > 0, 
R 
when (5.57) 


i) 5tdxdy =0. 
R 


In regions where ¢ > 0, 5f can have any sign and the inequality of (5.57) must 
be replaced by an equality. It then follows that 


(1/2EP){(T, +7, + 2(1 + v)(S? —T,T,)} = constant, (t>0) (5.58) 


or that the density of strain energy is uniformly distributed throughout the 
plate. 

The result of (5.58) does not, of course, imply uniformity of stress distribu- 
tion. However it is implied by it and so any solution, of Section 5.7, which has a 
uniform stress distribution and is also an elastic solution, gives an example of a 
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plate of maximum stiffness. The only examples which meet these requirements 
are those for the transmission of torque (Fig. 4.14) and for the transmission of 
pressure, using a uniform thickness, cited at the end of Section 5.7. 

The optimum design of plates for given strength was found in Section 5.7 to 
be closely related to the corresponding problem of Michell continua. This is 
not the case for design for maximum stiffness. Michell continua can be treated 
like the frameworks of Section 1.12, but as has just been seen the problem for 
plates is strictly speaking a difficult non-linear problem in differential equations. 
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